Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 164,101 T→A 14.3% L666Q (CTG→CAG)  hrpB → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913164,1010TA14.3% 418.5 / 55.8 140L666Q (CTG→CAG) hrpBputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base T (63/57);  new base A (9/11);  total (72/68)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 153 total aligned reads displayed.

TGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATC  >  NC_000913/163955‑164249
                                                                                                                                                  |                                                                                                                                                    
tGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCtgctgc                                                                                                                                                    <  4:836277/149‑1 (MQ=255)
  aTAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGACACGTGgctgct                                                                                                                                                      <  4:100516/145‑1 (MQ=255)
  aTAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCca                                                                                                                                                  >  1:237792/1‑149 (MQ=255)
      aGGTTTAAGCGTGCTCAACTGGACGGCGGAAGTGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcacgtt                                                                                                                                              <  3:449209/149‑4 (MQ=255)
      aGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACata                                                                                                                                              >  2:28176/1‑149 (MQ=255)
      aGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcacgtt                                                                                                                                              <  2:776855/149‑4 (MQ=255)
       ggTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACatat                                                                                                                                             <  3:828051/149‑1 (MQ=255)
       ggTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGCTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACatat                                                                                                                                             >  2:1275812/1‑149 (MQ=255)
         tttAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGa                                                                                                                                           <  3:1099343/149‑1 (MQ=255)
          ttAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGAc                                                                                                                                          >  1:80688/1‑149 (MQ=255)
              gCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGc                                                                                                                                      >  1:1238699/1‑149 (MQ=255)
                    tCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAt                                                                                                                                >  1:373039/1‑149 (MQ=255)
                    tCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTAGCAGCGCTGGAAACGTGGCAGCAGCcacgtt                                                                                                                                              >  2:208184/1‑132 (MQ=255)
                      aaCTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTc                                                                                                                              >  3:170200/1‑149 (MQ=255)
                      aaCTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTc                                                                                                                              >  3:70472/1‑149 (MQ=255)
                      aaCTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTc                                                                                                                              >  1:1061542/1‑149 (MQ=255)
                      aaCTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTc                                                                                                                              >  1:914907/1‑149 (MQ=255)
                       aCTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCa                                                                                                                             <  2:80688/149‑1 (MQ=255)
                           gACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTAc                                                                                                                         <  4:170200/149‑1 (MQ=255)
                               gcggaagcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCgg                                                                                                                     >  1:592060/1‑149 (MQ=255)
                                cggaagcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGc                                                                                                                    <  4:70472/149‑1 (MQ=255)
                                  gaagcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCt                                                                                                                  <  1:397768/149‑1 (MQ=255)
                                     gcggaaCAGCTATGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCca                                                                                                                                                  <  2:631642/114‑1 (MQ=255)
                                     gcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGcc                                                                                                                                                   <  2:504239/113‑1 (MQ=255)
                                     gcggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGcc                                                                                                                                                   >  1:504239/1‑113 (MQ=255)
                                      cggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGaaa                                                                                                              <  3:1145542/149‑1 (MQ=255)
                                       ggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGaaa                                                                                                              <  2:590684/148‑1 (MQ=255)
                                       ggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcac                                                                                                                                                 >  1:1261674/1‑113 (MQ=255)
                                       ggaaCAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcac                                                                                                                                                 <  2:1261674/113‑1 (MQ=255)
                                                   cTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACAttt                                                                                                 >  3:1036807/1‑149 (MQ=255)
                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAtt                                                                                                                               <  2:559357/118‑1 (MQ=255)
                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAtt                                                                                                                               >  1:559357/1‑118 (MQ=255)
                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCa                                                                                                                             <  4:84739/120‑1 (MQ=255)
                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcacgtt                                                                                                                                              >  2:434173/1‑100 (MQ=255)
                                                    ttGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcacgtt                                                                                                                                              <  1:434173/103‑4 (MQ=255)
                                                      gCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAtt                                                                                                                               >  1:237294/1‑116 (MQ=255)
                                                      gCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACAtt                                                                                                                               <  2:237294/116‑1 (MQ=255)
                                                      gCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATc                                                                                              <  2:485103/149‑1 (MQ=255)
                                                        gTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTa                                                                                                                                   <  2:1206475/110‑1 (MQ=255)
                                                        gTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTa                                                                                                                                   >  1:1206475/1‑110 (MQ=255)
                                                        gTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAg                                                                                            >  2:1397878/1‑149 (MQ=255)
                                                         tttGTTATGCGCCGCAAAGTGGTTGCCGTAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGAc                                                                                                     >  3:95020/1‑139 (MQ=255)
                                                         tttGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGAc                                                                                                     <  4:95020/139‑1 (MQ=255)
                                                             ttATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCt                                                                                                                  >  1:654721/1‑122 (MQ=255)
                                                                 gcgcCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcacgtt                                                                                                                                              >  4:1114819/1‑87 (MQ=255)
                                                                          aGTGGTTGCCGTAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTTGCTGCTGCCACATATGACTGGCGTACATTCACTACGCCGCCTGAAATCACTCGATATTTATCAGGCACTACGCGGATTACtt                                                                          <  2:373039/149‑1 (MQ=255)
                                                                            tGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGa                                                                        <  1:830349/149‑1 (MQ=255)
                                                                             ggTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGAt                                                                       <  3:651156/149‑1 (MQ=255)
                                                                               ttGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCAGCAGCcacgt                                                                                                                                               >  2:835843/1‑73 (MQ=255)
                                                                                     gAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGc                                                                                                                                      <  3:788932/78‑1 (MQ=255)
                                                                                     gAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAAt                                                               <  2:914907/149‑1 (MQ=255)
                                                                                        tatGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                            >  4:512278/1‑149 (MQ=255)
                                                                                        tatGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                            <  1:1397878/149‑1 (MQ=255)
                                                                                         atGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcag                                                           <  3:512278/149‑1 (MQ=255)
                                                                                         atGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgca                                                            <  1:1413830/148‑1 (MQ=255)
                                                                                         atGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgaa                                                            >  2:1413830/1‑146 (MQ=255)
                                                                                          tGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcagc                                                          >  3:1415268/1‑149 (MQ=255)
                                                                                           gACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcagca                                                         >  2:1209246/1‑149 (MQ=255)
                                                                                           gACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcagca                                                         >  3:1099186/1‑149 (MQ=255)
                                                                                            aCTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcagca                                                         <  1:646704/148‑1 (MQ=255)
                                                                                            aCTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATgcagca                                                         >  2:646704/1‑148 (MQ=255)
                                                                                             cTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGAt                                                                               <  2:28585/125‑1 (MQ=255)
                                                                                             cTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGAt                                                                               >  1:28585/1‑125 (MQ=255)
                                                                                             cTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAAc                                                       <  4:1099186/149‑1 (MQ=255)
                                                                                               ggCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGa                                                                                >  2:1464822/1‑122 (MQ=255)
                                                                                                  cAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACt                                                                           <  1:170503/124‑1 (MQ=255)
                                                                                                  cAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACt                                                                           >  2:170503/1‑124 (MQ=255)
                                                                                                  cAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTg                                                  >  4:1232855/1‑149 (MQ=255)
                                                                                                   aGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTgg                                                 >  1:1447985/1‑149 (MQ=255)
                                                                                                    gCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGa                                                <  3:1232855/149‑1 (MQ=255)
                                                                                                        ttGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGt                                            >  3:302191/1‑149 (MQ=255)
                                                                                                          gatgatgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCa                                                                                             >  1:347785/1‑98 (MQ=255)
                                                                                                             gatgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAAtt                                       <  4:1415268/149‑1 (MQ=255)
                                                                                                               tgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTggg                                                                   >  4:726509/1‑119 (MQ=255)
                                                                                                               tgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTggg                                                                   <  3:726509/119‑1 (MQ=255)
                                                                                                               tgaAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGc                                     >  2:365035/1‑149 (MQ=255)
                                                                                                                  aaGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGaa                                                                <  3:281234/119‑1 (MQ=255)
                                                                                                                  aaGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGaa                                                                >  4:281234/1‑119 (MQ=255)
                                                                                                                    gTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTggg                                                                   <  2:1229454/114‑1 (MQ=255)
                                                                                                                     tttATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTgcgc                               >  3:719360/1‑149 (MQ=255)
                                                                                                                         ttGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTa                           >  3:1043565/1‑149 (MQ=255)
                                                                                                                           ggCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAAt                                                               >  3:942509/1‑111 (MQ=255)
                                                                                                                           ggCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTg                                  >  1:587929/1‑140 (MQ=255)
                                                                                                                              aGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACtg                      <  4:1036807/149‑1 (MQ=255)
                                                                                                                                    tcatatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACt                       <  4:772660/136‑1 (MQ=255)
                                                                                                                                    ggAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGa                >  1:1449399/1‑149 (MQ=255)
                                                                                                                                     gAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGAtt                                                                      <  1:1222979/94‑1 (MQ=255)
                                                                                                                                       tatgTGGCAGCAGTCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATCGGGGaa                                                                >  4:1267022/4‑98 (MQ=255)
                                                                                                                                       tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATa                                              >  4:1229245/4‑116 (MQ=255)
                                                                                                                                       tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATa                                              <  3:1229245/113‑1 (MQ=255)
                                                                                                                                       tatgTGGCAGCAGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACATGATTGGGGaa                                                                <  3:1267022/95‑1 (MQ=255)
                                                                                                                                        aCGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACt                                                                           >  1:59653/1‑86 (MQ=255)
                                                                                                                                         cGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACtt                                                                          <  1:491292/86‑1 (MQ=255)
                                                                                                                                          gTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                             <  4:137906/114‑1 (MQ=255)
                                                                                                                                          gTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAg                                             >  3:137906/1‑114 (MQ=255)
                                                                                                                                          gTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGaa          >  1:880744/1‑149 (MQ=255)
                                                                                                                                           tGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGt                                            <  4:399016/114‑1 (MQ=255)
                                                                                                                                              ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTa                           <  1:203524/128‑1 (MQ=255)
                                                                                                                                              ctgctgCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCg      <  4:16529/149‑1 (MQ=255)
                                                                                                                                                  tgcCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGa                                                                                <  1:1038302/71‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGATAAAGGTTTAAGCGTGCTCAACTGGACGGCGGAAGCGGAACAGCTACGCTTGCGTTTGTTATGCGCCGCAAAGTGGTTGCCGGAATATGACTGGCCAGCGGTTGATGATGAAAGTTTATTGGCAGCGCTGGAAACGTGGCTGCTGCCACATATGACTGGCGTACATTCACTACGCGGCCTGAAATCACTCGACATTTATCAGGCACTACGCGGATTACTTGATTGGGGAATGCAGCAACGTCTGGATAGTGAATTGCCTGCGCATTACACTGTGCCGACGGGAAGCCGGATC  >  NC_000913/163955‑164249

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: