Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,421 C→A 12.1% 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,4210CA12.1% 316.3 / 28.3 121intergenic (+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 (51/55);  new base A (4/11);  total (55/66)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 129 total aligned reads displayed.

GGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCA  >  NC_000913/1426273‑1426570
                                                                                                                                                    |                                                                                                                                                     
ggAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTcgcg                                                                                                                                                      >  2:1370569/1‑150 (MQ=255)
 gAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCgg                                                                                                                                                     <  2:1019756/150‑1 (MQ=255)
     aGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTggg                                                                                                                                                 <  1:1037848/150‑1 (MQ=255)
     aGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTggg                                                                                                                                                 <  1:1478049/150‑1 (MQ=255)
      gtgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGt                                                                                                                                                <  1:236814/150‑1 (MQ=255)
       tgtgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtt                                                                                                                                               <  1:1116347/150‑1 (MQ=255)
        gtgtgtTGTTGTTTATTACGTACACTCGGTTACAGGGGTGTGTTGGGGTCTTTTATTAAATAGATGATGCTCCTGGAGAAAAGCATTGCGGGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGttt                                                                                                                                              >  2:1351767/1‑150 (MQ=255)
         tgtgtAGTTGTTATTTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttc                                                                                                                                             <  2:472849/150‑2 (MQ=255)
         tgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttt                                                                                                                                             <  1:778155/150‑1 (MQ=255)
         tgtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttt                                                                                                                                             >  2:293945/1‑150 (MQ=255)
          gtgtAGTTTTTAATTACGTACACTCTGTTACATAGATGAGATTGTGTCTTTAATTAAATCTATGATGCTCCTGGAGAAAAGCATTGCGTGTCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttca                                                                                                                                            <  2:235081/150‑3 (MQ=255)
          gtgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttca                                                                                                                                            <  2:1658417/150‑3 (MQ=255)
           tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcag                                                                                                                                           <  1:642113/150‑4 (MQ=255)
           tgtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcag                                                                                                                                           <  2:1087709/150‑4 (MQ=255)
            ttAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                          <  1:1611683/149‑5 (MQ=255)
            gtAGTTTTTAATACCGAACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                          <  1:273250/150‑5 (MQ=255)
            gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGtttcagg                                                                                                                                          <  1:628068/150‑5 (MQ=255)
            gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttg                                                                                                                                          <  1:529589/150‑1 (MQ=255)
            gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                          <  2:1595045/150‑5 (MQ=255)
            gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcagg                                                                                                                                          <  1:1717214/150‑5 (MQ=255)
            gtAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcaggc                                                                                                                                         <  1:720291/150‑6 (MQ=255)
             tAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgt                                                                                                                                         >  2:252587/1‑150 (MQ=255)
             tAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACCGCGAAGTGGGtttcaggc                                                                                                                                         <  1:1180930/150‑6 (MQ=255)
             tAGATGTTAATTACGTACACTCTGTTACAGAGATGTGATTGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgt                                                                                                                                         <  1:867812/150‑1 (MQ=255)
              aGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgtt                                                                                                                                        <  2:51411/150‑1 (MQ=255)
               gttgttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTttgttg                                                                                                                                       <  1:658319/150‑1 (MQ=255)
                  gttAATTACGTACTCTCTGTTACAGAGATGTGATGGTGTCTTTATTTTAATCGATGATGCTCCTTGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAATCCCACTTCGCGTTGGGTTTTTTGTTGTCa                                                                                                                                    <  1:1176714/150‑1 (MQ=255)
                   ttAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAg                                                                                                                                   <  2:1537794/150‑1 (MQ=255)
                    tAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAgg                                                                                                                                  >  2:314166/1‑150 (MQ=255)
                    tAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAgg                                                                                                                                  >  1:520726/1‑150 (MQ=255)
                        tACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCAGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGACTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGGGtt                                                                                                                              >  1:960323/1‑150 (MQ=255)
                               acTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTtaataa                                                                                                                       >  2:65467/1‑150 (MQ=255)
                                    gTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGc                                                                                                                  <  1:67988/150‑1 (MQ=255)
                                        catagaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCAGGACAAAAGCATTGCGTGGCCTCGTAATCGCTAAATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGCGTTTTTAGTTGTCATGAGTTTTAATAAATGGCatag                                                                                                              >  2:284838/4‑150 (MQ=255)
                                         agacATGTTATGTTGTGTTTAATTAAATCGATGATGCTCCTGGAGAAACGCATGGCGTGGCCGCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGc                                                                                                             <  1:1263346/150‑1 (MQ=255)
                                            gaTGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCaaa                                                                                                          >  2:1021679/1‑150 (MQ=255)
                                              tgtgTTGTTGTCTTTAATTAAATCGATGATGCTACTGGAGAAAAGCATTGCGTGGCCTCGTAATCTCTATATCTTCTATTATGTCGCATGAAACCCACATCGCGGTGGTTTTTTTGTTGGCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                        <  1:252587/150‑1 (MQ=255)
                                              tgtgATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAAcc                                                                                                        <  2:922702/150‑1 (MQ=255)
                                                  aTGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTcgc                                                                                                    <  1:1672333/150‑1 (MQ=255)
                                                     gtgtCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGtc                                                                                                 >  1:857701/1‑150 (MQ=255)
                                                      tgtCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGtct                                                                                                >  1:1374465/1‑150 (MQ=255)
                                                      tgtCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGtct                                                                                                >  1:783884/1‑150 (MQ=255)
                                                        tCTTTAACTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCg                                                                                              >  2:205156/1‑150 (MQ=255)
                                                        gcTTTAATTAAATAGATGATGCTCCTGGAGAAAAGCATTGCGTGGCGTCGTAATCGCTATATCTACTATTATGTCGCCTGACACCCACTTCGCGGTGGGGTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCg                                                                                              <  1:296245/149‑1 (MQ=255)
                                                         cTTTTATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTTTATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGc                                                                                             <  1:65467/150‑1 (MQ=255)
                                                          tttAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCa                                                                                            <  2:644698/150‑1 (MQ=255)
                                                            taattaaATCGATGATGCTCCTGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATATACTATTATGTCGCCTGAAACCCAATTCGCGGAGGGTTTTTTGTTGTAAGGAGTTTTAATAAATGGCAGAGAAAACCTCGAGTCTCACaata                                                                                         >  2:1592175/1‑150 (MQ=255)
                                                                   aTCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAattatt                                                                                   <  1:1001660/150‑1 (MQ=255)
                                                                     cGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGa                                                                                 >  2:582913/1‑150 (MQ=255)
                                                                           tgCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCAc                                                                           >  2:467875/1‑150 (MQ=255)
                                                                            gCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCTCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACt                                                                          <  1:647421/150‑1 (MQ=255)
                                                                                cTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAg                                                                      <  2:789061/150‑1 (MQ=255)
                                                                                cTGGAGAAAAGCATTGCGTGGCCTCGGAATCGCTTTATCTACTAGTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGGCTGGAGTTTTTATAAATGGCAGAGCACCCCTCGCGTCTCGCAAAAATTATTGATAGCACTggcg                                                                      >  2:332758/1‑148 (MQ=255)
                                                                                   gagaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCTCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGa                                                                   <  1:1547862/150‑1 (MQ=255)
                                                                                    agaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATATTTATTGATAGCACTGGAGCGaa                                                                  >  1:385307/1‑150 (MQ=255)
                                                                                     gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                                 >  1:1321372/1‑150 (MQ=255)
                                                                                     gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                                 <  1:1686208/150‑1 (MQ=255)
                                                                                     gaAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaa                                                                 <  2:519815/150‑1 (MQ=255)
                                                                                     gaAAAGAATTGCGTGGACTCGTAATCTCTATATCTACTATTATGTCGCCTGAAACCGACTTCGCGTTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTTATAGCACTGGAGCGaaa                                                                 <  1:1592175/150‑1 (MQ=255)
                                                                                       aaaGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGaaaaa                                                               >  2:1464726/1‑150 (MQ=255)
                                                                                           cATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                                  >  1:592632/1‑127 (MQ=255)
                                                                                           cATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTg                                                                                  <  2:592632/127‑1 (MQ=255)
                                                                                           cATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAaataat                                                           <  2:520726/150‑1 (MQ=255)
                                                                                            aTTGCGTGGTCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATg                                                          <  1:1592427/150‑1 (MQ=255)
                                                                                               gCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTg                                                       >  1:988447/1‑150 (MQ=255)
                                                                                                cGTGTCCTCGTAATCGCTATATCTACTATTATGTCTCATGAAACCCACTTCGCGGTTGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGa                                                      <  1:60840/150‑1 (MQ=255)
                                                                                                cGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGa                                                      >  1:925188/1‑150 (MQ=255)
                                                                                                  tGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACa                                                    >  1:1067219/1‑150 (MQ=255)
                                                                                                    gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                                  >  1:1145106/1‑150 (MQ=255)
                                                                                                    gCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAAt                                                  >  1:1567717/1‑150 (MQ=255)
                                                                                                     ccTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATc                                                 >  1:317323/1‑150 (MQ=255)
                                                                                                          tAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcc                                            <  1:293945/150‑1 (MQ=255)
                                                                                                           aaTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcct                                           >  2:724789/1‑150 (MQ=255)
                                                                                                           aaTCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGAcct                                           <  1:1351767/150‑1 (MQ=255)
                                                                                                               gCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACctcctc                                       >  1:1628034/1‑150 (MQ=255)
                                                                                                                cTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCa                                      >  1:1320864/1‑150 (MQ=255)
                                                                                                                  atatCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTGTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTAGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTTACATCATCAtt                                    >  1:524569/1‑150 (MQ=255)
                                                                                                                     tCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGt                                 >  1:1417292/1‑150 (MQ=255)
                                                                                                                        actaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaa                              >  1:52743/1‑150 (MQ=255)
                                                                                                                         ctaTTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCaaa                             <  2:1242040/150‑1 (MQ=255)
                                                                                                                           attatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAAt                           <  1:467875/150‑1 (MQ=255)
                                                                                                                           attatGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAAt                           <  2:1067219/150‑1 (MQ=255)
                                                                                                                              atGTCTCATGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGTAGTTTTAATAAATGTCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAACATGAc                        <  2:23484/150‑1 (MQ=255)
                                                                                                                              atGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGAc                        >  1:1590510/1‑150 (MQ=255)
                                                                                                                               tGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACg                       <  2:1354163/150‑1 (MQ=255)
                                                                                                                                gTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATATTTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGc                      >  1:1385732/1‑150 (MQ=255)
                                                                                                                                gTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGc                      <  1:560437/150‑1 (MQ=255)
                                                                                                                                   gCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgg                   <  2:1280516/150‑1 (MQ=255)
                                                                                                                                   gCCTGAAACCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGATGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgga                  >  1:1623145/1‑149 (MQ=255)
                                                                                                                                   acaaaaaaCCCACCGCGAAGTGGGTTTCTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAgg                   >  1:1700356/6‑150 (MQ=255)
                                                                                                                                    caaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTACTCATTAGTCAAAATGACGCAGGc                  >  2:845169/5‑150 (MQ=255)
                                                                                                                                     cTGAAACCCCCTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATACTTATTGATAGCTCTGGAGCGAAAAATAATGCTGACAATCGGCCCTCCTCATTAGTGACAATGAGGCAGGCt                 >  1:267054/1‑150 (MQ=255)
                                                                                                                                     aaaaaaCCCACCGCGAAGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCTCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCt                 >  2:710948/4‑150 (MQ=255)
                                                                                                                                      tGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTg                >  1:1496643/1‑150 (MQ=255)
                                                                                                                                                   gcgGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCATAGCAAACCTAGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAATAATAATGCTGACAATGTGAGCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAActgatg  <  1:284838/150‑2 (MQ=255)
                                                                                                                                                   gcgGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACtgctgc   <  2:1145106/150‑1 (MQ=255)
                                                                                                                                                    cgGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCa  <  2:857701/150‑1 (MQ=255)
                                                                                                                                                    cgGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCa  <  2:482490/150‑1 (MQ=255)
                                                                                                                                                    cgGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCa  >  2:1030745/1‑150 (MQ=255)
                                                                                                                                                    cgGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCa  >  2:1501168/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAGCAGTGTGTGTAGTTGTTAATTACGTACACTCTGTTACAGAGATGTGATGGTGTCTTTAATTAAATCGATGATGCTCCTGGAGAAAAGCATTGCGTGGCCTCGTAATCGCTATATCTACTATTATGTCGCCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTGGAGCGAAAAATAATGCTGACAATCTGACCTCCTCATTAGTCAAAATGACGCAGGCTGGGGAAACTGCTGCA  >  NC_000913/1426273‑1426570

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