Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,421 C→A 6.6% intergenic (+339/‑33) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,4210CA 6.6% 206.9 / 5.4 76intergenic (+339/‑33)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base C (43/28);  new base A (2/3);  total (45/31)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

AGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGA  >  NC_000913/1426272‑1426560
                                                                                                                                                     |                                                                                                                                           
aGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTTCAGAGATGTGATGTTGTCTTTTATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgc                                                                                                                                             >  1:142872/1‑150 (MQ=255)
aGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGATCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgc                                                                                                                                             >  2:666115/1‑150 (MQ=255)
 ggAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGAGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgcg                                                                                                                                            >  2:380802/1‑150 (MQ=255)
    gCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTg                                                                                                                                         <  2:588644/150‑1 (MQ=255)
     cAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTgg                                                                                                                                        <  2:627039/150‑1 (MQ=255)
       gtgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGt                                                                                                                                      <  2:305950/150‑1 (MQ=255)
         gtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGttt                                                                                                                                    >  1:10257/1‑150 (MQ=255)
            tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttttt                                                                                                                                 >  1:1147453/1‑150 (MQ=255)
            tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcag                                                                                                                                 <  1:984270/150‑4 (MQ=255)
              tAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgt                                                                                                                               <  2:703416/150‑1 (MQ=255)
              tAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcaggc                                                                                                                               <  2:816078/150‑6 (MQ=255)
                gttgttAATTACGTACACTCTTTTTCAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgttg                                                                                                                             <  2:237063/150‑1 (MQ=255)
                ggtgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTAAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATCTCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgttg                                                                                                                             <  2:1286801/148‑1 (MQ=255)
                     tAAATACGTACACTCTGTTACAGAGATGAGCTGGTGTCTTTAAATAAATCGATGATGCTCCTGGAGAAAAGCAATGCGTGGCCTCGTAATCGCTATATCTACTATTATGTAGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAgg                                                                                                                        <  2:289983/150‑1 (MQ=255)
                           cGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtttt                                                                                                                  >  1:1205497/1‑150 (MQ=255)
                           cGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGtttt                                                                                                                  <  1:890939/150‑1 (MQ=255)
                             tACACTCTGTTACAGAGATGTGATGGTGTCGTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaa                                                                                                                >  1:516585/1‑150 (MQ=255)
                              acacTCTGTTACAGAGATGTTATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaat                                                                                                               <  2:836755/150‑1 (MQ=255)
                                    tGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATgg                                                                                                         >  1:82830/1‑150 (MQ=255)
                                    tGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATgg                                                                                                         >  1:912496/1‑150 (MQ=255)
                                     gTTACAGAGCTGTGATGGGGTCTTTAATTAAATCGAGGATGCTCCAGGAGAAAAGCATTGCGTGGCCTAGTAATCGCTATCTCTACTAATATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGc                                                                                                        <  2:197213/150‑1 (MQ=255)
                                       tACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGATACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCag                                                                                                      <  2:768644/150‑1 (MQ=255)
                                       tACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTACTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATCATGTCGCCTGAAACCCACTTCTCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCag                                                                                                      >  1:689200/1‑150 (MQ=255)
                                        aCAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCaga                                                                                                     <  1:596491/150‑1 (MQ=255)
                                        aCAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGAAGTTTTAATAAATGGCaga                                                                                                     >  1:1086874/1‑150 (MQ=255)
                                              aTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAc                                                                                               <  1:846850/150‑1 (MQ=255)
                                               tgtgATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                              >  1:569470/1‑150 (MQ=255)
                                                gtgATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCt                                                                                             >  2:1123065/1‑150 (MQ=255)
                                                  gATGGTGTCTTTAATTAAATAGATGATGCTCCTGGAGAAAATCATTGAGTTGCCTCTTAATCGCTATATCTACTATTATGTCGCATGAAACCCACTTCGCGGTGCGTTTTTTGTTGTCAGGAGTTTTATTAAATGGCATAGCAAACCTca                                                                                           >  1:907937/1‑149 (MQ=255)
                                                    tGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTcgcg                                                                                         <  1:279275/150‑1 (MQ=255)
                                                     ggTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGt                                                                                        >  2:758738/1‑150 (MQ=255)
                                                     ggTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGt                                                                                        >  1:864191/1‑150 (MQ=255)
                                                     ggTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGAAGAGCAAACCTcgctt                                                                                        >  2:202256/1‑148 (MQ=255)
                                                         tCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCg                                                                                    <  1:1286898/150‑1 (MQ=255)
                                                          cTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                   >  1:797671/1‑150 (MQ=255)
                                                          cTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                   >  1:1276955/1‑150 (MQ=255)
                                                            ttaattaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaa                                                                                 <  1:1114185/150‑1 (MQ=255)
                                                               attaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaataa                                                                              >  2:1110815/1‑150 (MQ=255)
                                                               attaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaataa                                                                              >  2:895444/1‑150 (MQ=255)
                                                               attaaATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCaataa                                                                              >  2:199271/1‑150 (MQ=255)
                                                                      cGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACATCGCTTATCGAAATAattattca                                                                       >  1:416491/1‑148 (MQ=255)
                                                                       gatgatGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAACCCTCGCGTCTCGCAATAATTATTGAt                                                                      >  2:1037243/1‑150 (MQ=255)
                                                                       gatgatGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGAt                                                                      >  2:329393/1‑150 (MQ=255)
                                                                             gCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACt                                                                <  2:797671/150‑1 (MQ=255)
                                                                               tCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTgg                                                              <  1:199271/150‑1 (MQ=255)
                                                                                  tGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGc                                                           >  1:774077/1‑150 (MQ=255)
                                                                                  tGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGc                                                           <  2:912496/150‑1 (MQ=255)
                                                                                   ggAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCg                                                          >  2:1308755/1‑150 (MQ=255)
                                                                                   ggAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCg                                                          >  2:984892/1‑150 (MQ=255)
                                                                                     agaAAAGCATTGCGTTGCCTCGTAATCGCTATATCTACTATTATGTCTTCTGATTCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaa                                                        <  2:1003889/150‑1 (MQ=255)
                                                                                      taaaaGCAATGCGTTGCCTCGTTCTCGCTATCTCTACAATTATGTCGCCGGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGGTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                       <  2:122365/149‑1 (MQ=255)
                                                                                                gCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTg                                             <  1:1136278/150‑1 (MQ=255)
                                                                                                 cGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGa                                            >  1:178781/1‑150 (MQ=255)
                                                                                                     gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                        >  2:933316/1‑150 (MQ=255)
                                                                                                       cTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAATTAATGCTGACAATCt                                      >  2:616609/1‑150 (MQ=255)
                                                                                                             aTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcctc                                <  2:377099/150‑1 (MQ=255)
                                                                                                                gCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACctcctc                             >  2:1269509/1‑150 (MQ=255)
                                                                                                                    tatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTa                         <  2:1276955/150‑1 (MQ=255)
                                                                                                                     atCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAg                        >  1:1043406/1‑150 (MQ=255)
                                                                                                                     atCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTAGTCAGGAGTTTTAAAAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAg                        <  1:986982/150‑1 (MQ=255)
                                                                                                                     atCTACTATTATGTCGCCTGAAACCCACTTCGCGGGGTGTTTTTTGTTGTCAGTAGTTTTAATAAATGGCAGAGCAAACCTCGCGGCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAg                        <  2:967112/150‑1 (MQ=255)
                                                                                                                      tCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGt                       <  1:933316/150‑1 (MQ=255)
                                                                                                                        tactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCa                     >  2:954106/1‑150 (MQ=255)
                                                                                                                        tactaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCa                     >  1:292627/1‑150 (MQ=255)
                                                                                                                         actaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTTAAATCACTGGAGCGAAAAATTATGCTGACAATCTGACCTCCTCATTAGTCaa                    >  2:108321/1‑150 (MQ=255)
                                                                                                                          ctaTTATTTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGGCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaaa                   <  1:214434/150‑1 (MQ=255)
                                                                                                                           tattatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaaaa                  >  2:1111104/1‑150 (MQ=255)
                                                                                                                            attatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCTGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAAGTATTGATAGCAATGGAGCGATAAATAATGCTGACAATCTGACCTACTGATTAGTAAAAAt                 >  1:200996/1‑150 (MQ=255)
                                                                                                                                 gTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGc            >  1:834540/1‑150 (MQ=255)
                                                                                                                                  caacaaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACATCCTCATTAGTCAAAATGACGCa           >  2:319038/8‑150 (MQ=255)
                                                                                                                                   cGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGCCCTCCTCATTAGTCAAAATGACGCAg          >  2:822067/1‑150 (MQ=255)
                                                                                                                                   cGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAg          >  2:671158/1‑150 (MQ=255)
                                                                                                                                    gCCTGAAACCCACTTAGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTTACCTCCTCATTAGTCAAAATGACGCAgg         >  2:85408/1‑150 (MQ=255)
                                                                                                                                      aaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt       >  2:588288/4‑150 (MQ=255)
                                                                                                                                       tGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATTACGCAGTCTg      >  1:475523/1‑150 (MQ=255)
                                                                                                                                       tGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTg      >  2:18366/1‑150 (MQ=255)
                                                                                                                                         aaaCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTggg    <  1:594742/150‑1 (MQ=255)
                                                                                                                                           aCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGa  >  2:1200729/1‑150 (MQ=255)
                                                                                                                                           aCCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGATGACAAACTGAACTCCTCATTAGTCAAAATGACGAAGGCTGGGGa  >  2:584525/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
AGGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGA  >  NC_000913/1426272‑1426560

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