Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 928,302 C→A 5.5% T298T (ACG→ACT)  cydC ← glutathione/L‑cysteine ABC exporter subunit CydC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913928,3020CA 5.5% 243.7 / 2.9 73T298T (ACG→ACT) cydCglutathione/L‑cysteine ABC exporter subunit CydC
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base A (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

AGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATC  >  NC_000913/928155‑928449
                                                                                                                                                   |                                                                                                                                                   
aGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                                                                                                                   >  2:1510422/1‑150 (MQ=255)
aGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGt                                                                                                                                                   >  1:1402069/1‑150 (MQ=255)
      ctgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTgg                                                                                                                                             >  1:847203/1‑150 (MQ=255)
      ctgaacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTgg                                                                                                                                             <  2:651979/150‑1 (MQ=255)
         aacATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGc                                                                                                                                          <  2:751333/150‑1 (MQ=255)
          acATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGcc                                                                                                                                         >  2:1458522/1‑150 (MQ=255)
              cccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTg                                                                                                                                     >  2:1316166/1‑150 (MQ=255)
               ccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGCTACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                                                                                                                    >  1:775794/1‑150 (MQ=255)
               ccGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGc                                                                                                                                    >  1:416351/1‑150 (MQ=255)
                cgaaaCGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                   >  1:875286/4‑150 (MQ=255)
                cGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                   <  2:1567117/150‑1 (MQ=255)
                cGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                   <  1:1418918/150‑1 (MQ=255)
                cGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCt                                                                                                                                   >  1:875287/1‑150 (MQ=255)
                      gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                             >  1:335177/1‑150 (MQ=255)
                      gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                             >  1:341098/1‑150 (MQ=255)
                      gTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                             <  2:847203/150‑1 (MQ=255)
                      gTCAGCGAAACGCGATCGGCAACACGAGTTTCGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAc                                                                                                                             <  2:1460786/150‑1 (MQ=255)
                         aGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAAcgcg                                                                                                                          <  2:1402069/150‑1 (MQ=255)
                           cGAAACGCGATCGGCAACACGAGTTTGGGTATCAGCAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGCTGCCAGTGCTTCAAACGCGGc                                                                                                                        >  1:1390420/1‑150 (MQ=255)
                            gAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCt                                                                                                                       <  2:1341580/150‑1 (MQ=255)
                              aaCGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTaa                                                                                                                     >  1:1367162/1‑150 (MQ=255)
                                cgcgATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAAcg                                                                                                                   <  2:1213607/150‑1 (MQ=255)
                                   gATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACgcgc                                                                                                                <  1:1455350/150‑1 (MQ=255)
                                     tCGGCAACACGAGTCTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGGCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATCCACCCGCTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCag                                                                                                              <  2:1239662/150‑1 (MQ=255)
                                      cGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCaga                                                                                                             >  2:978548/1‑150 (MQ=255)
                                       ggCAACACGAGTTGGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCagaa                                                                                                            <  1:541792/150‑1 (MQ=255)
                                        gCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCagaag                                                                                                           <  1:284578/150‑1 (MQ=255)
                                            cacGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACaa                                                                                                       >  1:622520/1‑150 (MQ=255)
                                                aGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCGGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAg                                                                                                   >  2:1553601/1‑150 (MQ=255)
                                                     gggTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCaa                                                                                              <  2:923979/150‑1 (MQ=255)
                                                        tATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTa                                                                                           >  2:373751/1‑150 (MQ=255)
                                                        tATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTa                                                                                           >  2:1020492/1‑150 (MQ=255)
                                                         aTCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTaa                                                                                          <  1:250357/150‑1 (MQ=255)
                                                            aGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAAcgc                                                                                       >  2:304971/1‑150 (MQ=255)
                                                             ggAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAAcgcg                                                                                      <  1:616097/150‑1 (MQ=255)
                                                              gAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACgcgc                                                                                     >  1:1560179/1‑150 (MQ=255)
                                                                aaGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCg                                                                                   <  2:696526/150‑1 (MQ=255)
                                                                aaGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACACGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCg                                                                                   >  2:231979/1‑150 (MQ=255)
                                                                  ggTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCggg                                                                                 >  1:192266/1‑150 (MQ=255)
                                                                     gACCTCGGGTTTTTGATCCGTTAAGTCAGCGATACGACCGGCAGAGGCAATGACTTGCCCCAGATTCTGAAATGCACCCGTTTCTTTTGCCAGTGCTTCAAACGCGGCTAACGCGCTGAAGACAAACAGGGCTATTAACGCGCCGGGTTg                                                                              <  2:1449763/150‑1 (MQ=255)
                                                                     gACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTg                                                                              >  2:252960/1‑150 (MQ=255)
                                                                      aCCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGa                                                                             <  2:579843/150‑1 (MQ=255)
                                                                       ccTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAg                                                                            <  1:1316166/150‑1 (MQ=255)
                                                                       ccTCCGG‑TTTTGATCCGTTAAGTCATAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGc                                                                           <  2:1060109/150‑1 (MQ=255)
                                                                        cTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGc                                                                           >  2:379228/1‑150 (MQ=255)
                                                                               ttttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTgccg                                                                    >  2:865405/1‑150 (MQ=255)
                                                                                 ttGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTgccgcc                                                                  <  1:1510422/150‑1 (MQ=255)
                                                                                  tGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCa                                                                 >  1:672215/1‑150 (MQ=255)
                                                                                        gTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAAcgccg                                                           <  2:775794/150‑1 (MQ=255)
                                                                                         ttAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAAcgccgc                                                          <  1:304971/150‑1 (MQ=255)
                                                                                          tAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACgccgcc                                                         >  1:878374/1‑150 (MQ=255)
                                                                                          tAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACgccgcc                                                         <  2:54873/150‑1 (MQ=255)
                                                                                          tAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACgccgcc                                                         >  2:651530/1‑150 (MQ=255)
                                                                                           aaGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCa                                                        >  1:806881/1‑150 (MQ=255)
                                                                                                 gagaTACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGcc                                                  >  1:86315/1‑150 (MQ=255)
                                                                                                     tacgtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATcc                                              <  2:806881/150‑1 (MQ=255)
                                                                                                     tacgtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATcc                                              <  1:745941/150‑1 (MQ=255)
                                                                                                       cgtacgGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCcac                                            >  2:580100/1‑150 (MQ=255)
                                                                                                                gagGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAgg                                   <  2:1404819/150‑1 (MQ=255)
                                                                                                                  ggCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGAt                                 >  1:1125683/1‑150 (MQ=255)
                                                                                                                    cAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGTATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCa                               >  1:1039767/1‑150 (MQ=255)
                                                                                                                      aTGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCAcc                             >  1:188735/1‑150 (MQ=255)
                                                                                                                          cTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTACCGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTa                         <  2:878374/150‑1 (MQ=255)
                                                                                                                             gCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAAcg                      >  2:874229/1‑150 (MQ=255)
                                                                                                                                 cAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGcc                  <  2:416351/150‑1 (MQ=255)
                                                                                                                                       cTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag            <  2:763345/150‑1 (MQ=255)
                                                                                                                                       cTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag            >  2:96312/1‑150 (MQ=255)
                                                                                                                                       agcactggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGag            >  2:1275841/8‑150 (MQ=255)
                                                                                                                                         gAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGagca          <  1:379228/150‑1 (MQ=255)
                                                                                                                                          actggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATg              >  1:1174956/5‑145 (MQ=255)
                                                                                                                                          actggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATg              <  2:1174956/141‑1 (MQ=255)
                                                                                                                                          actggCACCAGTAACGGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGagcag         <  2:954385/146‑1 (MQ=255)
                                                                                                                                               cACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTa    <  1:978548/150‑1 (MQ=255)
                                                                                                                                                 cccGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATc  >  1:1327696/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGTGAACTGAACATCCCGTAACGTCAGCGAAACGCGATCGGCAACACGAGTTTGGGTATCAGGAAAGGTGACCTCCGGTTTTTGATCCGTTAAGTCAGAGATACGTACGGCAGAGGCAATGACTTGCCCCAGATGCTGAAATGCACCCGTTACTGGTGCCAGTGCTTCAAACGCGGCTAACGCGCAGAAGACAAACAGGGCAATTAACGCGCCGGGTTGAGCATTGCCGCCAACGCCGCCAGACGCCATCCACAGCATCAGGATCACCGCTAACGCGCCAATGAGCAGCATTATC  >  NC_000913/928155‑928449

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