Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 867,917 T→C 6.4% L133P (CTG→CCG)  gsiA → glutathione ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,9170TC 6.4% 239.6 / 5.7 78L133P (CTG→CCG) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (33/40);  new base C (3/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

TCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGTGATGATTGCGATGG  >  NC_000913/867771‑868063
                                                                                                                                                  |                                                                                                                                                  
tCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTg                                                                                                                                                   >  7:97077/1‑148 (MQ=255)
  gcgAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCAt                                                                                                                                                >  3:174694/1‑149 (MQ=255)
  gcgAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCAt                                                                                                                                                >  4:74653/1‑149 (MQ=255)
    gAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCa                                                                                                                                              <  3:363882/149‑1 (MQ=255)
    gAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCa                                                                                                                                              <  3:64626/149‑1 (MQ=255)
      aGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGa                                                                                                                                            <  7:19256/149‑1 (MQ=255)
          aTTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGACCCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGc                                                                                                                                        >  6:193056/1‑149 (MQ=255)
          aTTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGc                                                                                                                                        >  5:200587/1‑149 (MQ=255)
          aTTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGc                                                                                                                                        <  4:227383/149‑1 (MQ=255)
           ttGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGc                                                                                                                                        >  3:352302/1‑148 (MQ=255)
             gAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCa                                                                                                                                      <  8:108284/148‑1 (MQ=255)
             gAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCa                                                                                                                                      <  8:213517/148‑1 (MQ=255)
              aaCTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGt                                                                                                                                    >  1:372496/1‑149 (MQ=255)
              aaCTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGt                                                                                                                                    >  5:377018/1‑149 (MQ=255)
               aCTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTc                                                                                                                                   >  4:57152/1‑149 (MQ=255)
                cTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCg                                                                                                                                  >  1:289458/1‑149 (MQ=255)
                 ttAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGt                                                                                                                                 >  1:4912/1‑149 (MQ=255)
                    gCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTgaa                                                                                                                              <  7:375941/149‑1 (MQ=255)
                    gCGAGCAGAACGCTGCACAAATGCGCCATGCTCGCGGTGCGGATATGGCGATGATATTACAGGAGCCGATGACCTCGCGGCACCCGGTAGTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTgaa                                                                                                                              <  2:277236/149‑1 (MQ=255)
                             aCGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATg                                                                                                                     <  3:347929/149‑1 (MQ=255)
                             aCGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATg                                                                                                                     <  2:4912/149‑1 (MQ=255)
                                cTGCACAAATGCGCCATGTTCGCGGGGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTc                                                                                                                  <  2:346863/149‑1 (MQ=255)
                                  gCACAAATGCGCCATGTTCGCGGTGCGGATATCGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTc                                                                                                                  >  6:12388/1‑147 (MQ=255)
                                   cacaAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAg                                                                                                               >  1:268672/1‑149 (MQ=255)
                                    acaAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACATATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAgg                                                                                                              >  8:53359/1‑149 (MQ=255)
                                      aaaTGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGc                                                                                                             >  3:362102/1‑148 (MQ=255)
                                        aTGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTTCCGAATCAAGTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCg                                                                                                            >  5:394898/1‑147 (MQ=255)
                                        aTGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGaa                                                                                                          <  5:193056/149‑1 (MQ=255)
                                           cgcCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGc                                                                                                        <  3:116115/148‑1 (MQ=255)
                                             ccATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGt                                                                                                                                 <  3:139614/121‑1 (MQ=255)
                                             ccATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGt                                                                                                                                 >  4:139614/1‑121 (MQ=255)
                                             ccATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGa                                                                                                     <  7:151396/149‑1 (MQ=255)
                                              cATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGAt                                                                                                    <  5:321452/149‑1 (MQ=255)
                                               aTGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATg                                                                                                   >  1:14983/1‑149 (MQ=255)
                                                 gTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCt                                                                                                 >  2:29255/1‑149 (MQ=255)
                                                 gTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCCGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCt                                                                                                 >  3:91533/1‑149 (MQ=255)
                                                       ggTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGTATGCTGGATCa                                                                                           <  3:57152/149‑1 (MQ=255)
                                                        gTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAg                                                                                          <  4:352302/149‑1 (MQ=255)
                                                         tGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAgg                                                                                         <  2:268672/149‑1 (MQ=255)
                                                           cGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGt                                                                                        <  6:200587/148‑1 (MQ=255)
                                                           cGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTa                                                                                       >  7:32841/1‑149 (MQ=255)
                                                            ggATATTGGGATGATATTTCAGGTGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTAc                                                                                      <  2:135250/149‑1 (MQ=255)
                                                                  tgCGATGATATTTCAGGAGCCGATGACATCGCTGAACACGGTATTTACTGTGGGTGAACAGATTGCCTAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTc                                                                                >  3:82078/2‑149 (MQ=255)
                                                                     gatgatATTTCAGGAGCCGATGACATCGCTGAAACCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGc                                                                                                        <  4:306314/122‑1 (MQ=255)
                                                                     gatgatATTTCAGGAGCCGATGACATCGCTGAAACCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGc                                                                                                        >  3:306314/1‑122 (MQ=255)
                                                                      atgatATTTCAGTAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCCGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAg                                                                                          <  1:266846/135‑1 (MQ=255)
                                                                      atgatATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGa                                                                            <  8:97077/149‑1 (MQ=255)
                                                                      atgatATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCCGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAg                                                                                          <  1:272649/135‑1 (MQ=255)
                                                                      atgatATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCCGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAg                                                                                          >  2:266846/1‑135 (MQ=255)
                                                                      atgatATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCCGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAg                                                                                          >  2:272649/1‑135 (MQ=255)
                                                                               cAGGAGCCGATGACATCGCTGAACCCGGTTTTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACaaa                                                                    <  1:107162/148‑1 (MQ=255)
                                                                                  gAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCAt                                                                >  4:171284/1‑149 (MQ=255)
                                                                                  gAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCAt                                                                <  6:377018/149‑1 (MQ=255)
                                                                                   aGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCAtt                                                               <  2:55394/149‑1 (MQ=255)
                                                                                        aTGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGAGAAGCGGATGCTGGATCAGGGACGCATTCCTGAGGCACAAACCAATCttt                                                           >  6:194948/1‑148 (MQ=255)
                                                                                           aCATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCATGTACGCATTCCTGAGGCACAAACCATTCTTTCAc                                                        <  4:362102/148‑1 (MQ=255)
                                                                                                      cccGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATc                                            <  3:171284/149‑1 (MQ=255)
                                                                                                       ccGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATc                                            >  2:411750/1‑148 (MQ=255)
                                                                                                        cGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCaa                                          <  1:29255/149‑1 (MQ=255)
                                                                                                                    tgGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCg                                <  7:53359/147‑1 (MQ=255)
                                                                                                                      ggTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGa                             <  3:234166/148‑1 (MQ=255)
                                                                                                                          aaCAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATgcgc                        >  4:155797/1‑149 (MQ=255)
                                                                                                                          aaCAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATgcgc                        <  1:411750/149‑1 (MQ=255)
                                                                                                                             aGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCa                                                                 <  1:199172/105‑1 (MQ=255)
                                                                                                                             aGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCa                                                                 >  2:199172/1‑105 (MQ=255)
                                                                                                                             aGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGGGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAg                     <  2:78755/149‑1 (MQ=255)
                                                                                                                                ttGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGa                  <  7:170362/149‑1 (MQ=255)
                                                                                                                                   ccGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGtg               <  8:32841/149‑1 (MQ=255)
                                                                                                                                     gAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGtgat             >  3:311019/1‑149 (MQ=255)
                                                                                                                                      aaTCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGtgat             >  4:145472/1‑148 (MQ=255)
                                                                                                                                       aTCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGtgatga           >  7:378847/1‑149 (MQ=255)
                                                                                                                                       aTCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGtgatga           <  5:12388/149‑1 (MQ=255)
                                                                                                                                         cAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGTGATGAtt         >  7:204687/1‑149 (MQ=255)
                                                                                                                                          aaTTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCg                                <  4:57164/125‑1 (MQ=255)
                                                                                                                                          aaTTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCg                                >  3:57164/1‑125 (MQ=255)
                                                                                                                                           aTTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAACCCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGTGATGATTGc       >  1:313188/1‑149 (MQ=255)
                                                                                                                                               gTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGTGATGATTGCGATg   <  3:74653/149‑1 (MQ=255)
                                                                                                                                                tCTGCATCAGAAAGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCTGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGTGATGATTGCGATgg  <  6:196831/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATATGGCGATGATATTTCAGGAGCCGATGACATCGCTGAACCCGGTATTTACTGTGGGTGAACAGATTGCCGAATCAATTCGTCTGCATCAGAACGCCAGTCGTGAAGAAGCGATGGTCGAGGCGAAGCGGATGCTGGATCAGGTACGCATTCCTGAGGCACAAACCATTCTTTCACGTTATCCGCATCAACTCTCTGGCGGGATGCGCCAGCGAGTGATGATTGCGATGG  >  NC_000913/867771‑868063

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: