Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,488,426 A→T 7.8% intergenic (+26/+135) lptG → / ← yjgR lipopolysaccharide transport system protein LptG/DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,488,4260AT 7.8% 559.6 / 27.1 180intergenic (+26/+135)lptG/yjgRlipopolysaccharide transport system protein LptG/DUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base A (85/81);  new base T (5/9);  total (90/90)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

Only 100 of 184 total aligned reads displayed.

CGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCC  >  NC_000913/4488277‑4488574
                                                                                                                                                     |                                                                                                                                                    
 gTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGaaa                                                                                                                                                     >  2:1199882/1‑150 (MQ=255)
  tACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAg                                                                                                                                                    >  2:108912/1‑150 (MQ=255)
    cTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGGGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCt                                                                                                                                                  <  2:1271328/150‑1 (MQ=255)
    cTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCt                                                                                                                                                  >  2:633113/1‑150 (MQ=255)
     tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTc                                                                                                                                                 <  2:92575/150‑1 (MQ=255)
     tGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTc                                                                                                                                                 >  1:159842/1‑150 (MQ=255)
      ggACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                                                                                                                                                <  2:626072/150‑1 (MQ=255)
      ggACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTcc                                                                                                                                                >  2:581038/1‑150 (MQ=255)
       gACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTccc                                                                                                                                               >  2:130564/1‑150 (MQ=255)
         ccAGATCTTCGGCCCGCTGACGTCGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCtt                                                                                                                                             <  2:159024/150‑1 (MQ=255)
           aGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTtta                                                                                                                                           <  2:742172/150‑1 (MQ=255)
             aTCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttatt                                                                                                                                         <  2:496456/150‑1 (MQ=255)
               cTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTtattat                                                                                                                                       >  1:1004380/1‑150 (MQ=255)
                  cGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCTTTCGCTCCCTttattcgttc                                                                                                                                    <  1:113769/150‑6 (MQ=255)
                   ggCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGtt                                                                                                                                   <  2:314527/150‑1 (MQ=255)
                    gCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTa                                                                                                                                  <  2:124257/150‑1 (MQ=255)
                     cccGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAg                                                                                                                                 <  2:429787/150‑1 (MQ=255)
                      ccGCTGACGTTGGTTTATTGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATCAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGc                                                                                                                                <  2:394537/150‑1 (MQ=255)
                        gCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCaa                                                                                                                              <  2:1394465/150‑1 (MQ=255)
                           gACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTtattat                                                                                                                                       <  1:1392691/138‑1 (MQ=255)
                           gACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTtattat                                                                                                                                       >  2:1392691/1‑138 (MQ=255)
                             cGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgtgt                                                                                                                         >  2:331620/1‑150 (MQ=255)
                              gTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGtgtg                                                                                                                        >  1:848704/1‑150 (MQ=255)
                               ttGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGc                                                                                                                       <  2:963918/150‑1 (MQ=255)
                                   tttATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgtgt                                                                                                                         <  1:294036/144‑1 (MQ=255)
                                   tttATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAgtgt                                                                                                                         >  2:294036/1‑144 (MQ=255)
                                     tATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt                                                                                                                 >  1:576061/1‑150 (MQ=255)
                                       tGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCACACTTTGCtaacaata                                                                                                               <  2:47334/150‑8 (MQ=255)
                                       tGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCa                                                                                                               <  2:979176/150‑1 (MQ=255)
                                         gCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAtt                                                                                                             >  1:1197683/1‑150 (MQ=255)
                                         acagtgcgccgATGATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCAtt                                                                                                             >  1:372631/7‑150 (MQ=255)
                                             cccGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTg                                                                                                         <  2:230258/150‑1 (MQ=255)
                                             cccGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTg                                                                                                         >  1:1137292/1‑150 (MQ=255)
                                                gccgATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGaaa                                                                                                      <  2:447413/150‑1 (MQ=255)
                                                      catcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataat                                                                                                >  1:324623/1‑150 (MQ=255)
                                                      catcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataat                                                                                                >  1:346828/1‑150 (MQ=255)
                                                      catcGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAaataat                                                                                                >  1:804649/1‑150 (MQ=255)
                                                           gcgcACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTaaaa                                                                                           <  2:671334/150‑1 (MQ=255)
                                                                cTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCaa                                                                                      <  2:1291692/150‑1 (MQ=255)
                                                                  gTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAata                                                                                    >  2:401723/1‑150 (MQ=255)
                                                                  gTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCATGAAAAATAATTAAAATTCAata                                                                                    >  2:576044/1‑150 (MQ=255)
                                                                      acAAGCGACAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATAtt                                                                                >  2:1030307/2‑150 (MQ=255)
                                                                        aaGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGc                                                                              <  1:1400487/150‑1 (MQ=255)
                                                                             ccAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtgg                                                                         <  2:94044/150‑1 (MQ=255)
                                                                               aGCTTCTTCTTAATCAGCCTGTGGCTGTTAACGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtgggg                                                                       <  1:512703/150‑1 (MQ=255)
                                                                                  ttcttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAgtggggtgg                                                                    <  1:657521/150‑1 (MQ=255)
                                                                                    cttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGc                                                                  <  2:236259/150‑1 (MQ=255)
                                                                                    cttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGc                                                                  >  2:1326790/1‑150 (MQ=255)
                                                                                     ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTc                                                                                                                <  1:1258494/103‑1 (MQ=255)
                                                                                     ttcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTCGGCa                                                                 >  2:1057536/1‑150 (MQ=255)
                                                                                      tcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCAt                                                                >  2:1206073/1‑150 (MQ=255)
                                                                                      tcttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCAt                                                                >  1:758498/1‑150 (MQ=255)
                                                                                       gttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATg                                                               <  1:711682/149‑1 (MQ=255)
                                                                                       cttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATg                                                               <  2:651190/150‑1 (MQ=255)
                                                                                       cttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATg                                                               >  1:338671/1‑150 (MQ=255)
                                                                                        ttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCa                                                                             <  2:796868/135‑1 (MQ=255)
                                                                                        ttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCa                                                                             >  1:796868/1‑135 (MQ=255)
                                                                                        ttAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGa                                                              <  2:348889/150‑1 (MQ=255)
                                                                                         tAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGAt                                                             <  1:180330/150‑1 (MQ=255)
                                                                                             cAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGa                                                         >  2:1342028/1‑150 (MQ=255)
                                                                                             cAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGa                                                         >  2:1121289/1‑150 (MQ=255)
                                                                                              aGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAc                                                        <  2:368167/150‑1 (MQ=255)
                                                                                               gCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcg                                                       >  1:499266/1‑150 (MQ=255)
                                                                                               gCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcg                                                       >  1:111425/1‑150 (MQ=255)
                                                                                                ccTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGAcgc                                                      <  1:13539/150‑1 (MQ=255)
                                                                                                  tgtgGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACgcgc                                                    <  1:1328515/150‑1 (MQ=255)
                                                                                                   gtgGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGcc                                                   <  2:159842/150‑1 (MQ=255)
                                                                                                   gtgGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGcc                                                   <  1:108912/150‑1 (MQ=255)
                                                                                                     ggCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGtt                                                                                                                 >  1:98225/1‑86 (MQ=255)
                                                                                                       cTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCg                                               <  1:1206073/150‑1 (MQ=255)
                                                                                                        tGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGAt                                                             <  1:225135/135‑1 (MQ=255)
                                                                                                        tGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGAt                                                             >  2:225135/1‑135 (MQ=255)
                                                                                                         gTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTc                                             >  2:800273/1‑150 (MQ=255)
                                                                                                           tAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGc                                           >  2:435138/1‑150 (MQ=255)
                                                                                                            aaTGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATATTTAAAATTCAATATATTGCAGTGGGGTGGGCAGGATAAGACGCGCCAGCGTCGCa                                          >  2:353392/1‑150 (MQ=255)
                                                                                                                agaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAg                                      >  2:1039322/1‑150 (MQ=255)
                                                                                                                 gaAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAgg                                     >  2:279834/1‑150 (MQ=255)
                                                                                                                     aTCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCaaa                                 >  2:448095/1‑150 (MQ=255)
                                                                                                                      tCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAg                                <  2:1197683/150‑1 (MQ=255)
                                                                                                                        gTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGAc                              <  1:879498/150‑1 (MQ=255)
                                                                                                                            ccGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAAt                          <  1:943998/150‑1 (MQ=255)
                                                                                                                                 tGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGc                     >  1:136936/1‑150 (MQ=255)
                                                                                                                                  gAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGcc                    <  1:279834/150‑1 (MQ=255)
                                                                                                                                  gAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGcc                    >  2:901984/1‑150 (MQ=255)
                                                                                                                                  gAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGcc                    >  2:1284766/1‑150 (MQ=255)
                                                                                                                                   caataaTAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCt                   >  1:1056180/5‑150 (MQ=255)
                                                                                                                                   aaCGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCt                   <  2:298907/150‑1 (MQ=255)
                                                                                                                                     cGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGa                 >  1:1069006/1‑150 (MQ=255)
                                                                                                                                       aaTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCa                                                                 <  2:60194/100‑1 (MQ=255)
                                                                                                                                       aaTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCa                                                                 >  1:60194/1‑100 (MQ=255)
                                                                                                                                       aaTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACt               <  1:1002675/150‑1 (MQ=255)
                                                                                                                                        aTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACtt              <  2:400700/150‑1 (MQ=255)
                                                                                                                                        aTAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACtt              <  2:1004380/150‑1 (MQ=255)
                                                                                                                                         tAAAGGGAGCTTTCGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCg                                            <  2:343909/119‑1 (MQ=255)
                                                                                                                                              ggAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttc        <  1:1236945/150‑1 (MQ=255)
                                                                                                                                               gAGCTAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGAATTATAATACTCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttct       >  1:528511/1‑150 (MQ=255)
                                                                                                                                               gAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACcttct       <  1:435138/150‑1 (MQ=255)
                                                                                                                                                   gAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTtctc   >  1:621494/1‑150 (MQ=255)
                                                                                                                                                    aaaGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTcc  <  2:1224060/150‑1 (MQ=255)
                                                                                                                                                    aaaGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTcc  <  2:1167926/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGTACTGGACCAGATCTTCGGCCCGCTGACGTTGGTTTATGGCATCCCGCCGATCATCGGCGCACTGTTGCCAAGCGCCAGCTTCTTCTTAATCAGCCTGTGGCTGTTAATGAGAAAATCGTAACCGAATGAACGAATAAAGGGAGCGAAAGCTCCCTTTATTATTGTTAGCAAAGTGTGCTTCGTTCATTCCTGAAAAATAATTAAAATTCAATATATTGCAGTGGGGTGGGCATGATAAGACGCGCCAGCGTCGCATCAGGCAAAGACAAATCTCGCCTGACTTACCTTCTTCTCC  >  NC_000913/4488277‑4488574

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: