Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,396 T→A 5.3% I87N (ATT→AAT)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913435,3960TA 5.3% 764.3 / 22.2 206I87N (ATT→AAT) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (101/94);  new base A (6/5);  total (107/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

Only 100 of 206 total aligned reads displayed.

TAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCC  >  NC_000913/435255‑435543
                                                                                                                                             |                                                                                                                                                   
tAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                                                                <  2:990388/147‑1 (MQ=255)
 aaaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGt                                                                                                                                             >  2:401093/1‑149 (MQ=255)
 aaaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGt                                                                                                                                             >  4:323029/1‑149 (MQ=255)
  aaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTa                                                                                                                                            >  3:407729/1‑149 (MQ=255)
      aCGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                                                                        <  1:89815/149‑1 (MQ=255)
                 cTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                    >  1:828708/1‑142 (MQ=255)
                 cTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                    <  2:828708/142‑1 (MQ=255)
                 cTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTa                                                                                                                                            >  4:1193589/1‑134 (MQ=255)
                 cTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTa                                                                                                                                            <  3:1193589/134‑1 (MQ=255)
                    tACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAg                                                                                                                          >  3:979646/1‑149 (MQ=255)
                    tACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAg                                                                                                                          <  3:363265/149‑1 (MQ=255)
                       ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                                                                        >  4:1153629/1‑132 (MQ=255)
                             gAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAc                                                                                                                                           >  1:1042656/1‑123 (MQ=255)
                                ctgctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATa                                                                                                              >  2:1332090/1‑149 (MQ=255)
                                 tgctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                                                                        >  1:1015284/1‑122 (MQ=255)
                                 tgctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATAc                                                                                                             <  4:407729/149‑1 (MQ=255)
                                   ctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgt                                                                                                                    <  4:754473/140‑1 (MQ=255)
                                   ctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaa                                                                                                           >  3:1409299/1‑149 (MQ=255)
                                     ggCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                                                                               <  1:38562/127‑1 (MQ=255)
                                     ggCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                                                                               >  2:38562/1‑127 (MQ=255)
                                      gCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                                                                         <  1:285961/132‑1 (MQ=255)
                                      gCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                                                                         >  2:1155602/1‑132 (MQ=255)
                                      gCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                                                                         <  1:1155602/132‑1 (MQ=255)
                                      gCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                                                                         >  2:285961/1‑132 (MQ=255)
                                          gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                    <  3:744092/149‑1 (MQ=255)
                                                 gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGCAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg                                                                                             >  1:957250/1‑149 (MQ=255)
                                                    tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                                                                                                         >  1:254141/1‑118 (MQ=255)
                                                      aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcg                                                                                                                                                 <  2:1040557/92‑1 (MQ=255)
                                                      aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCt                                                                                                                                  <  3:1141395/107‑1 (MQ=255)
                                                       tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGa                                                                                       <  1:401093/149‑1 (MQ=255)
                                                        aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                                                                                      >  2:1067324/1‑149 (MQ=255)
                                                          cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTAc                                                                                                                                           <  4:960645/94‑1 (MQ=255)
                                                          cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGc                                                                                                                  <  1:499034/119‑1 (MQ=255)
                                                          cGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGc                                                                                                                  >  2:499034/1‑119 (MQ=255)
                                                            cATACCTCGACGGAATGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTAc                                                                                                                                           >  3:108249/1‑92 (MQ=255)
                                                             aTACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCt                                                                                                                                               <  2:427949/87‑1 (MQ=255)
                                                             aTACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                                                                 >  3:846983/1‑149 (MQ=255)
                                                                 cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                                                                                      <  3:434117/140‑1 (MQ=255)
                                                                 cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                                                                                      >  4:434117/1‑140 (MQ=255)
                                                                     aCGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                                                                         <  1:1067324/149‑1 (MQ=255)
                                                                               gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCg                                                                            >  1:1211592/1‑136 (MQ=255)
                                                                               gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgcg                                                               >  2:265664/1‑149 (MQ=255)
                                                                                   ccATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                    >  1:722119/1‑108 (MQ=255)
                                                                                   ccATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                    <  2:722119/108‑1 (MQ=255)
                                                                                     aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGc                                                                                                                  <  4:887284/92‑1 (MQ=255)
                                                                                     aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGc                                                                                                                  >  3:887284/1‑92 (MQ=255)
                                                                                     aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                    <  1:675921/106‑1 (MQ=255)
                                                                                     aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                    >  2:675921/1‑106 (MQ=255)
                                                                                        ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc                                                                                              >  2:975992/1‑109 (MQ=255)
                                                                                        ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                                                                           >  4:693112/1‑112 (MQ=255)
                                                                                        ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                      >  4:474259/1‑149 (MQ=255)
                                                                                        ccTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                      >  1:804960/1‑149 (MQ=255)
                                                                                         cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt                                                                                      <  2:1457851/116‑1 (MQ=255)
                                                                                         cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAg                                                     >  1:647480/1‑149 (MQ=255)
                                                                                         cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAct                                                                                                                                          >  1:1108022/1‑63 (MQ=37)
                                                                                         cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAct                                                                                                                                          <  2:1108022/64‑2 (MQ=37)
                                                                                           gTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGcc                                                   <  2:647480/149‑1 (MQ=255)
                                                                                              ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTg                                                                                                                                                   <  1:872846/50‑1 (MQ=255)
                                                                                              ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgtg                                                                                                                   <  2:930570/82‑1 (MQ=255)
                                                                                              ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaaa                                                                                                          >  1:936890/1‑91 (MQ=255)
                                                                                               cGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                               <  3:474259/149‑1 (MQ=255)
                                                                                                     ctgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGt                                         <  1:241255/149‑1 (MQ=255)
                                                                                                        gaagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCaa                                      <  2:804960/149‑1 (MQ=255)
                                                                                                          agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgt                                                                                                                    >  4:1337406/1‑69 (MQ=255)
                                                                                                          agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgt                                                                                                                    <  3:1337406/69‑1 (MQ=255)
                                                                                                          agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTc                                           >  2:173315/1‑142 (MQ=255)
                                                                                                          agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                    <  4:378003/149‑1 (MQ=255)
                                                                                                          agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                    <  2:1294473/149‑1 (MQ=255)
                                                                                                           gaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGtt                                                                    >  1:373033/1‑116 (MQ=255)
                                                                                                           gaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGtt                                                                    <  2:373033/116‑1 (MQ=255)
                                                                                                            aaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGc                                  <  3:663080/149‑1 (MQ=255)
                                                                                                              cTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg                                                                  >  4:181472/1‑115 (MQ=255)
                                                                                                              cTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCg                                                                  <  3:181472/115‑1 (MQ=255)
                                                                                                                    cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa                          <  4:1409299/149‑1 (MQ=255)
                                                                                                                       gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                            <  4:1126797/144‑1 (MQ=255)
                                                                                                                       gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                            <  3:733359/144‑1 (MQ=255)
                                                                                                                       gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                            >  4:733359/1‑144 (MQ=255)
                                                                                                                        tAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg                                                                                             <  1:159640/78‑1 (MQ=255)
                                                                                                                        tAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATaaa                      <  2:907670/149‑1 (MQ=255)
                                                                                                                          gAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc                                                                                              >  1:560643/1‑75 (MQ=255)
                                                                                                                          gAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc                                                                                              <  2:560643/75‑1 (MQ=255)
                                                                                                                           aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                               <  1:316480/121‑1 (MQ=255)
                                                                                                                           aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                               <  1:976022/121‑1 (MQ=255)
                                                                                                                           aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                               >  2:316480/1‑121 (MQ=255)
                                                                                                                            aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                      >  1:990050/1‑113 (MQ=255)
                                                                                                                             aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcagc                 >  4:854972/1‑149 (MQ=255)
                                                                                                                              aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa                                                          >  1:1354510/1‑107 (MQ=255)
                                                                                                                              aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAgcagca                >  2:1010980/1‑149 (MQ=255)
                                                                                                                              aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCAATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAgcagca                <  3:959793/149‑1 (MQ=255)
                                                                                                                               aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTc                                           >  2:949894/1‑121 (MQ=255)
                                                                                                                                gCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAcc              >  4:559983/1‑149 (MQ=255)
                                                                                                                                  aGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtg            >  1:1425943/1‑149 (MQ=255)
                                                                                                                                    tacagcgcAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                                                                         <  1:1296381/82‑1 (MQ=255)
                                                                                                                                      cTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcag                  >  1:914179/1‑139 (MQ=255)
                                                                                                                                      cTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAt        >  1:751187/1‑149 (MQ=255)
                                                                                                                                      cTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGGGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAt        <  2:247389/149‑1 (MQ=255)
                                                                                                                                       tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt       <  2:751187/149‑1 (MQ=255)
                                                                                                                                           cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc   <  2:980430/149‑1 (MQ=255)
                                                                                                                                           cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc   <  2:881017/149‑1 (MQ=255)
                                                                                                                                            aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc  >  4:829533/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
TAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCC  >  NC_000913/435255‑435543

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