Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 777,208 A→C 9.9% K289N (AAA→AAC)  tolA → membrane anchored protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913777,2080AC 9.9% 473.2 / 3.1 172K289N (AAA→AAC) tolAmembrane anchored protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base A (109/41);  new base C (3/19);  total (112/60)
Fisher's exact test for biased strand distribution p-value = 1.48e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 174 total aligned reads displayed.

GAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAATAA  >  NC_000913/777062‑777357
                                                                                                                                                  |                                                                                                                                                     
gagaAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCg                                                                                                                                                    >  2:3451943/1‑150 (MQ=255)
   aaaGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAATAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGcc                                                                                                                                                 >  2:1978689/1‑150 (MQ=255)
        tGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCaga                                                                                                                                            >  1:1004456/1‑150 (MQ=255)
        tGCAGACGACAACAAAGCAGCAGAAAATGCGGCTGCGGAAACGGAAGCAGCTGATAAGAAAGAAGCGGAATAAAAAGCCGCCGCAGACAAAAAAGCGGCTGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCaga                                                                                                                                            <  2:2147230/150‑1 (MQ=255)
               gaaaaaaaaCCAGCATAAAACGGGGCTGCTGAAACGGCCGCAGCTTAAAAGAACGCAGCGGCAGAAAATGCCGCCGCAGACAAAAAAGCGGCTGCGGCAAAAGCTGAAGCTGAAAAAGCCGCTGCAGCAAACGCGGCCGCAGAGGCAgat                                                                                                                                     <  2:1943018/147‑1 (MQ=255)
               gACAACAAAGCAGCATAAAAAGCGGCTGCTGAAACGGCAGCAGCTGATAAGACAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgat                                                                                                                                     <  2:2206664/150‑1 (MQ=255)
               gACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgat                                                                                                                                     >  2:1617587/1‑150 (MQ=255)
                aCAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatg                                                                                                                                    >  1:1979268/1‑150 (MQ=255)
                aCAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatg                                                                                                                                    >  1:2819935/1‑150 (MQ=255)
                  aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCCGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                  >  2:1969027/1‑150 (MQ=255)
                  aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAGAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                  >  1:3861098/1‑150 (MQ=255)
                  aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                  >  2:1580832/1‑150 (MQ=255)
                  aaaaaaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAgatgat                                                                                                                                  >  2:4116826/1‑150 (MQ=255)
                      aaGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCGGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCGGAAAAAGCCGCTGCCGCCCCCGCGGCCGCAGAGGCAGAGGATAttt                                                                                                                              >  2:3603281/1‑150 (MQ=37)
                        gcagcaGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTc                                                                                                                            <  2:3410512/150‑1 (MQ=255)
                         cagcagAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGAGGATATTTTCg                                                                                                                           >  2:3447871/1‑150 (MQ=255)
                           gcagAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGt                                                                                                                         >  2:871221/1‑150 (MQ=255)
                           gaagAAAACGCGGCTGCTGAAACGGCAGCAGCTGAAAAGAACGCAGCGGCAGAAACCGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAACGCTGCAGCTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGt                                                                                                                         <  2:58933/148‑1 (MQ=255)
                                aaaaGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCt                                                                                                                    <  2:3720911/150‑1 (MQ=255)
                                 aaaGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTa                                                                                                                   >  2:2436766/1‑150 (MQ=255)
                                 aaaGCGGCTGCTGAAAAGGCACCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGATGCAGCTGAAAAAGCCGCTGCAGCAAACGCGGCCGCAGAGGAAGATGATATTTTCGGTTAGCTa                                                                                                                   >  2:2436769/1‑150 (MQ=38)
                                    gCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGc                                                                                                                <  2:2312298/150‑1 (MQ=255)
                                    gCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGc                                                                                                                <  2:2243070/150‑1 (MQ=255)
                                    gCGGCTGCTGAAAAGGCAGAAGATTAAAACACTGCAGGGGATGAAAAAGCCGCCGCAGACAAAAACGCGTAAGCGGAAACAGCTGCAGCTGAAAAAGCCGCGGCAGAAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGc                                                                                                                <  2:1852473/150‑1 (MQ=255)
                                      ggCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGctc                                                                                                              >  2:4082255/1‑150 (MQ=255)
                                        ctgctgAAAAGGCAGCAGCTGATACGACAGCAGCGGCAGAAAACGCCGCCGCAGACAAAAAAGCGGCAGCGGCAACAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTg                                                                                                            <  2:3630999/150‑1 (MQ=255)
                                        ctgctgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTg                                                                                                            >  1:843954/1‑150 (MQ=255)
                                            tgAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTaa                                                                                                        >  1:1417917/1‑150 (MQ=255)
                                              aaaaGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCACAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGa                                                                                                      >  2:454798/1‑150 (MQ=255)
                                                   gctgctgcTTATAATAACGCCGCGCAATAAAACGCCGCCCCACACAAAAACGCGGCAGCGGAAAACGCTGCAGTTTAAAAAGCGGCTGCAGAAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCa                                                                                                 <  2:1741663/144‑1 (MQ=255)
                                                   gcagcaGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCa                                                                                                 <  2:3982297/150‑1 (MQ=255)
                                                   gcagcaGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCa                                                                                                 >  1:1555528/1‑150 (MQ=255)
                                                   gcagaaGCTCATAACAAACATGCGCAAAAAAACGCCGCCGCAGACAAAAAAGCGGCAGCGGAAAAAGCTGCAGCTTAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCa                                                                                                 <  2:1460775/150‑1 (MQ=255)
                                                             aTAAGAACGCAGCGGCATAAAACGCCGCCGAAGACAAACAAGCGGCAGCGGCAAACGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTACGAATGCACCGAAAACgg                                                                                       <  2:2442337/150‑1 (MQ=255)
                                                              tAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggg                                                                                      >  1:4077092/1‑150 (MQ=255)
                                                               aaGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCTGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggg                                                                                     <  2:1574305/150‑1 (MQ=255)
                                                               aaGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTG‑AAAAGCCGCATCATCAAAACCCGCCGCAGGGGCAGATGAAATTTTAGGTTATCAAATCAAATGTAAGAATTCACCCGAAAAggggg                                                                                    >  2:696375/1‑150 (MQ=255)
                                                                aGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAACGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACggggg                                                                                    <  2:1345409/150‑1 (MQ=255)
                                                                 gAAAGCAGCGGCAGAAAAAGCCGCCGCAGCCCAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCACAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACgggggg                                                                                   >  2:2674160/1‑150 (MQ=255)
                                                                  aaaGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGa                                                                                  >  2:1373021/1‑150 (MQ=255)
                                                                  aaaGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGa                                                                                  >  1:3493100/1‑150 (MQ=255)
                                                                     gcagcGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCATAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAggg                                                                               >  2:1500500/1‑150 (MQ=255)
                                                                       agcGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGc                                                                             >  1:2148806/1‑150 (MQ=255)
                                                                         cGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAACAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGa                                                                           >  2:3592610/1‑150 (MQ=255)
                                                                         cGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGa                                                                           >  1:888006/1‑150 (MQ=255)
                                                                          ggCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGaa                                                                          >  1:417161/1‑150 (MQ=255)
                                                                           gCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAATAGCTGCAGCTGAAAATGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGaaa                                                                         >  2:1506531/1‑150 (MQ=255)
                                                                             aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGACAgg                                                                       >  2:3109949/1‑150 (MQ=255)
                                                                             aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                       >  1:2917447/1‑150 (MQ=255)
                                                                             aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                       >  1:6956/1‑150 (MQ=255)
                                                                             aGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAgg                                                                       >  2:2566439/1‑150 (MQ=255)
                                                                              gAAAACGCCGCCGCAGACAAAAAAGCGGCCGCGCCAAACGCTGCAGCTGAAAACGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAggg                                                                      <  2:1564836/150‑1 (MQ=255)
                                                                              gAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAggg                                                                      <  2:3770392/150‑1 (MQ=255)
                                                                               aaaaaGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGa                                                                     >  1:94711/1‑150 (MQ=255)
                                                                                  aaGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACa                                                                  >  1:3968163/1‑150 (MQ=255)
                                                                                     ccgccgCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATg                                                               >  2:940047/1‑150 (MQ=255)
                                                                                      cgccgcAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGc                                                              >  2:3728756/1‑150 (MQ=255)
                                                                                           cAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGc                                                         >  2:2232866/1‑150 (MQ=255)
                                                                                             gaaaaaaaaGGGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGAAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGACGAAAGGGAACAATGCTTCGCCt                                                       <  2:3285957/147‑1 (MQ=255)
                                                                                               cAAAAACGCGCCTGCGCAAAACGCTCCACCTAAAAACGCCCGTGCACCAAACGCGGCCGCAGTGGCAGTTGTTTTTTTCGGTGAGCTAAGCTCTGGTATGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGc                                                     <  1:3948956/150‑1 (MQ=255)
                                                                                                aaaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGcc                                                    >  2:2957135/1‑150 (MQ=255)
                                                                                                aaaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGcc                                                    >  1:489179/1‑150 (MQ=255)
                                                                                                aaaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGcc                                                    >  2:2555660/1‑150 (MQ=255)
                                                                                                 aaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCg                                                   >  2:818046/1‑150 (MQ=255)
                                                                                                 aaaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCg                                                   >  2:1089002/1‑150 (MQ=255)
                                                                                                  aaaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCgg                                                  >  1:1750216/1‑150 (MQ=255)
                                                                                                   aaaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCggg                                                 >  1:2519130/1‑150 (MQ=255)
                                                                                                   aaaGCGGCAGCGGCAAAAGCTGCAGCTG‑AAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGa                                                >  1:843031/1‑150 (MQ=255)
                                                                                                    aaGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGa                                                >  2:1724464/1‑150 (MQ=255)
                                                                                                       cggcagcgtcaAAAGCTGCAGCTGAAAACGCCGCTGCAGAAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGATTGCACCGAAAACGGGGGGAGGGGCTAAAGGGAACAATGCTTCGCCTGCCGGGAGTg                                             <  2:1431136/150‑1 (MQ=255)
                                                                                                       cggcagcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTg                                             >  2:3890247/1‑150 (MQ=255)
                                                                                                           agcggcaAAAGCTCCAGCTCAAAAACCCGCTGCACCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCCGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTaa                                         <  1:2710531/150‑1 (MQ=255)
                                                                                                            gcggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAAt                                        >  1:4020045/1‑150 (MQ=255)
                                                                                                             cggcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATa                                       <  2:2407307/150‑1 (MQ=255)
                                                                                                              ggcaACCGCTGCAGCTCAAACCGCCGCTGCACCAACCGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAAGGCACCGAACACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATAc                                      <  1:3368479/150‑1 (MQ=255)
                                                                                                               gcaAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACt                                     >  1:3319853/1‑150 (MQ=255)
                                                                                                               gcaAAAGCTGCAGATGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACt                                     >  1:3828883/1‑150 (MQ=255)
                                                                                                                    aGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaaaa                                <  2:3094499/150‑1 (MQ=255)
                                                                                                                    aGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTaaaaa                                >  1:2447805/1‑150 (MQ=255)
                                                                                                                       tGCCGCTGAAAATGCCGGGGCAGAAAACGCGGCCGCTGCGGCTGTTGATATTTTCGGTGAGCTAAGCTGGGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACaa                             <  1:3273761/150‑1 (MQ=255)
                                                                                                                       tGCAGCTGAAAAAGCCGCTGCAGCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACaa                             <  2:2554258/150‑1 (MQ=255)
                                                                                                                        gCAGCTGAAAACGCGGCTGCAGCAACCGCGGCCGCTGAGGCAGATGATATTTTCGGTGAGCTATGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAAt                            <  1:242067/150‑1 (MQ=255)
                                                                                                                         cAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATg                           <  2:1360634/150‑1 (MQ=255)
                                                                                                                         cAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATg                           >  2:543212/1‑150 (MQ=255)
                                                                                                                           gCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgc                         <  2:1008465/150‑1 (MQ=255)
                                                                                                                            cTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgcg                        >  1:3131218/1‑150 (MQ=255)
                                                                                                                            cTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgcg                        >  1:3921268/1‑150 (MQ=255)
                                                                                                                             tGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGgcgc                       >  2:1758209/1‑150 (MQ=255)
                                                                                                                               aaaaaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCAt                     >  2:1233249/1‑150 (MQ=255)
                                                                                                                                  aaGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAg                  >  1:3705129/1‑150 (MQ=255)
                                                                                                                                   cgCCGCTGCAGCAAACGCGGCCGCAGCGGCAGATGATATTTTCGGTGAGCTAAGCTCTGTTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgg                 <  1:1114724/149‑1 (MQ=255)
                                                                                                                                     ccGCTGCAGCAAAAGCGGCCGCAGAGGCTGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               <  2:798253/150‑1 (MQ=255)
                                                                                                                                     ccGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               >  2:2795753/1‑150 (MQ=255)
                                                                                                                                     ccGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAgggg               >  1:1604742/1‑150 (MQ=255)
                                                                                                                                         tGCAGCAAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGTTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGa           <  1:2526741/150‑1 (MQ=255)
                                                                                                                                         tGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGa           <  2:4077092/150‑1 (MQ=255)
                                                                                                                                            agcaAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat        >  1:31545/1‑150 (MQ=255)
                                                                                                                                            accaAAAGCGGCCGCAGAGGCAGATTATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGatat        <  1:2136646/148‑1 (MQ=255)
                                                                                                                                              caAACGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCa      <  1:2643365/150‑1 (MQ=255)
                                                                                                                                                  aGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAATaa  <  2:538676/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GAGAAAGCTGCAGCCGACAAAAAAGCAGCAGAAAAAGCGGCTGCTGAAAAGGCAGCAGCTGATAAGAAAGCAGCGGCAGAAAAAGCCGCCGCAGACAAAAAAGCGGCAGCGGCAAAAGCTGCAGCTGAAAAAGCCGCTGCAGCAAAAGCGGCCGCAGAGGCAGATGATATTTTCGGTGAGCTAAGCTCTGGTAAGAATGCACCGAAAACGGGGGGAGGGGCGAAAGGGAACAATGCTTCGCCTGCCGGGAGTGGTAATACTAAAAACAATGGCGCATCAGGGGCCGATATCAATAA  >  NC_000913/777062‑777357

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