Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,309,009 C→T 5.1% intergenic (‑85/+24) nlpI ← / ← pnp lipoprotein NlpI/polynucleotide phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,309,0090CT 5.1% 505.5 / 9.6 156intergenic (‑85/+24)nlpI/pnplipoprotein NlpI/polynucleotide phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (77/71);  new base T (5/3);  total (82/74)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 160 total aligned reads displayed.

TATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTT  >  NC_000913/3308861‑3309157
                                                                                                                                                    |                                                                                                                                                    
tATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCccgccg                                                                                                                                                     >  2:1686861/1‑150 (MQ=255)
tATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCccgccg                                                                                                                                                     <  1:838997/150‑1 (MQ=255)
    aCTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGc                                                                                                                                                 <  2:1477074/150‑1 (MQ=255)
        cATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAg                                                                                                                                             >  2:1126700/1‑150 (MQ=255)
          tCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAggg                                                                                                                                           <  1:518071/150‑1 (MQ=255)
           ccTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAATAGGCTTAATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCTCCGCAGCGGAGGGc                                                                                                                                          <  2:258208/150‑1 (MQ=255)
           ccTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGc                                                                                                                                          >  1:1205392/1‑150 (MQ=255)
            cTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGct                                                                                                                                         <  2:1439199/150‑2 (MQ=255)
            cTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCGGCGGAGGGCa                                                                                                                                         <  1:154498/150‑1 (MQ=255)
            cTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCa                                                                                                                                         >  2:1388989/1‑150 (MQ=255)
            cTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCa                                                                                                                                         >  2:2391653/1‑150 (MQ=255)
             tGCAAGCGTAAGTACTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCaa                                                                                                                                        >  1:2368947/1‑150 (MQ=255)
                aaGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCTCCGCTGCGGCGGGGGGctttta                                                                                                                                     <  2:325314/150‑6 (MQ=255)
                aaGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATg                                                                                                                                     <  2:930583/150‑1 (MQ=255)
                 aGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCTCCGCTGCGGCggggggcttt                                                                                                                                    <  2:2161470/150‑8 (MQ=255)
                      aaGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAAcc                                                                                                                               <  1:1498661/150‑1 (MQ=255)
                       aGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCt                                                                                                                              >  1:2080277/1‑150 (MQ=255)
                           cTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACt                                                                                                                          >  1:1055745/1‑150 (MQ=255)
                             gTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCg                                                                                                                        <  2:1238125/150‑1 (MQ=255)
                              tCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGc                                                                                                                       <  2:882608/150‑1 (MQ=255)
                               cGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGcc                                                                                                                      >  2:1806686/1‑150 (MQ=255)
                                gCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGccc                                                                                                                     >  1:676595/1‑150 (MQ=255)
                                gCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGccc                                                                                                                     >  2:683849/1‑150 (MQ=255)
                                gCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGccc                                                                                                                     >  1:1054695/1‑150 (MQ=255)
                                gCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGccc                                                                                                                     >  1:1963107/1‑150 (MQ=255)
                                 cAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTc                                                                                                                         <  2:1265164/145‑1 (MQ=255)
                                 cAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTc                                                                                                                         >  1:1265164/1‑145 (MQ=255)
                                 cAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCt                                                                                                                    >  2:2116905/1‑150 (MQ=255)
                                    cGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGtt                                                                                                                 <  1:727719/150‑1 (MQ=255)
                                     gAAACACCAGAGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATTGCAACCTTACTCGCCCTGTTc                                                                                                                <  1:1011723/150‑1 (MQ=255)
                                      aaaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTca                                                                                                               >  2:1546852/1‑150 (MQ=255)
                                      aaaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTca                                                                                                               <  1:430089/150‑1 (MQ=255)
                                       aaCACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTcag                                                                                                              <  1:2331158/150‑1 (MQ=255)
                                         cacCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTcagca                                                                                                            >  1:57186/1‑150 (MQ=255)
                                           ccAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCagcagc                                                                                                          >  1:1728063/1‑150 (MQ=255)
                                                 cAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGc                                                                                                    >  2:1394592/1‑150 (MQ=255)
                                                   aaaaaGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCtt                                                                                                  >  1:88643/1‑150 (MQ=255)
                                                            tCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGtgctg                                                                                         <  1:261644/150‑1 (MQ=255)
                                                             cATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCACCGGAAGCTCCGGctgctga                                                                                        <  1:2078894/150‑2 (MQ=255)
                                                              aTTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGCTGCTgca                                                                                       >  1:1168874/1‑150 (MQ=255)
                                                                  cccACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGcagcag                                                                                   >  2:812464/1‑150 (MQ=255)
                                                                   ccACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAgg                                                                                  >  1:2366231/1‑150 (MQ=255)
                                                                      cTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTg                                                                               >  2:432237/1‑150 (MQ=255)
                                                                       tCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                              <  2:1410670/150‑1 (MQ=255)
                                                                       tCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                              <  2:846636/150‑1 (MQ=255)
                                                                       ccccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGTGTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTga                                                                              <  2:1922390/149‑1 (MQ=255)
                                                                        cccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCACCCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTgag                                                                             >  2:1546234/1‑150 (MQ=255)
                                                                        cccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTgag                                                                             >  1:304164/1‑150 (MQ=255)
                                                                        cccGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTgag                                                                             <  2:894412/150‑1 (MQ=255)
                                                                         ccGAAGACCACGGTTGAATGAACGTCCTGTTCCGGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTgaga                                                                            <  1:1214942/150‑1 (MQ=255)
                                                                          cGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGAc                                                                           <  1:1394592/150‑1 (MQ=255)
                                                                             aGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGc                                                                        <  1:194696/150‑1 (MQ=255)
                                                                               aCCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTc                                                                      <  2:30926/150‑1 (MQ=255)
                                                                                   cGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTc                                                                  <  1:1894146/150‑1 (MQ=255)
                                                                                   cGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTc                                                                  <  1:195726/150‑1 (MQ=255)
                                                                                    ggTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCg                                                                 <  2:1282270/150‑1 (MQ=255)
                                                                                         aatgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACATTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGcttct                                                            >  2:1966082/1‑150 (MQ=255)
                                                                                          atgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGcttctt                                                           <  1:2211984/150‑1 (MQ=255)
                                                                                          atgaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGCTGCTGCAGCAGGTTGAGACTGCTCAGGCGctgatt                                                           >  2:1086742/1‑146 (MQ=255)
                                                                                            gaaCGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTa                                                         <  2:725377/150‑1 (MQ=255)
                                                                                               cGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATg                                                      <  1:1003898/150‑1 (MQ=255)
                                                                                                 tCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCt                                                    >  1:2402793/1‑150 (MQ=255)
                                                                                                 tCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCt                                                    <  1:1126700/150‑1 (MQ=255)
                                                                                                  ccTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTc                                                   <  1:2347829/150‑1 (MQ=255)
                                                                                                  ccTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTc                                                   <  2:1814233/150‑1 (MQ=255)
                                                                                                   cTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGGGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCa                                                  >  1:1173566/1‑150 (MQ=255)
                                                                                                   cTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCa                                                  <  1:1806686/150‑1 (MQ=255)
                                                                                                       tCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGAct                                              <  2:1841005/150‑2 (MQ=255)
                                                                                                          cGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGAt                                           <  1:1546852/150‑1 (MQ=255)
                                                                                                            gTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATAc                                         >  2:207613/1‑150 (MQ=255)
                                                                                                            gTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATAc                                         >  2:2262046/1‑150 (MQ=255)
                                                                                                              tGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACgg                                       <  2:828906/150‑1 (MQ=255)
                                                                                                                cTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGcc                                     >  2:2137335/1‑150 (MQ=255)
                                                                                                                  aaCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCt                                   >  1:1461062/1‑150 (MQ=255)
                                                                                                                   aCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATa                                          <  1:173540/142‑1 (MQ=255)
                                                                                                                   aCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTg                                  <  1:1993968/150‑1 (MQ=255)
                                                                                                                   aCAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGAAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATa                                          >  2:173540/1‑142 (MQ=255)
                                                                                                                    cAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGGAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTgg                                 <  1:708360/150‑1 (MQ=255)
                                                                                                                      aGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCg                               <  2:653647/150‑1 (MQ=255)
                                                                                                                       ggCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGa                              >  2:2293421/1‑150 (MQ=255)
                                                                                                                         cGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATc                            <  1:326451/150‑1 (MQ=255)
                                                                                                                         cGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATc                            >  2:598113/1‑150 (MQ=255)
                                                                                                                          gTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCa                           >  1:1588653/1‑150 (MQ=255)
                                                                                                                                 ccGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCa                    >  1:159967/1‑150 (MQ=255)
                                                                                                                                  cGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAg                   >  1:367159/1‑150 (MQ=255)
                                                                                                                                   ggTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGa                  >  1:766959/1‑150 (MQ=255)
                                                                                                                                    gTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGaa                 >  1:1510745/1‑150 (MQ=255)
                                                                                                                                    gTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGaa                 <  2:2017699/150‑1 (MQ=255)
                                                                                                                                    gTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGaa                 <  1:1686861/150‑1 (MQ=255)
                                                                                                                                    gTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGaa                 >  1:1381478/1‑150 (MQ=255)
                                                                                                                                      tAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt               <  1:2145925/150‑1 (MQ=255)
                                                                                                                                      tAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACt               <  1:2137335/150‑1 (MQ=255)
                                                                                                                                       aaaaGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGc                                                                        <  2:1603096/92‑1 (MQ=255)
                                                                                                                                           gCCCTCCGCTGCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAc          >  1:1781778/1‑150 (MQ=255)
                                                                                                                                            ccctccGCTGCGGCGGGGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAcc         >  1:950020/1‑150 (MQ=255)
                                                                                                                                            ccccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCAcc         <  2:2368947/150‑1 (MQ=255)
                                                                                                                                              ccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGTTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCgg       <  2:74430/150‑1 (MQ=255)
                                                                                                                                              ccccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACAGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCgg       >  1:1946470/1‑150 (MQ=255)
                                                                                                                                               cccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGt      <  1:2391653/150‑1 (MQ=255)
                                                                                                                                               cccGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGt      >  2:514723/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TATTACTGCATCCTGCAAGCGTAAGTGCTGTCGCAACGAAACACCAGCGCAAAAAAGGCTTCATTTCCCACTCCCGAAGACCACGGTTGAATGAACGTCCTGTTCCCGGTTGCTAACAAGGCGTCCTGCCCGGTTAAAAGCCCCCCGCCGCAGCGGAGGGCAAATGGCAACCTTACTCGCCCTGTTCAGCAGCCGGAGCTTCCGGTGCTGCAGCAGGTTGAGACTGCTCAGTCGCTTCTTTAATGCTCAGACGGATACGGCCCTGGCGATCAACTTCCAGAACTTTCACCGGTACTT  >  NC_000913/3308861‑3309157

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: