Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,396 T→A 5.6% 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.6% 316.0 / 7.5 89I87N (ATT→AAT) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (53/30);  new base A (3/2);  total (57/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCT  >  NC_000913/435253‑435544
                                                                                                                                               |                                                                                                                                                    
tGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                                                                 >  4:349082/1‑149 (MQ=255)
tGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                                                                 >  1:616653/1‑149 (MQ=255)
tGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTgcgc                                                                                                                                                 >  4:323257/1‑149 (MQ=255)
 gTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCt                                                                                                                                                >  3:77502/1‑149 (MQ=255)
    aaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATTCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTa                                                                                                                                             >  2:674195/1‑149 (MQ=255)
    aaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTa                                                                                                                                             <  2:139398/149‑1 (MQ=255)
         cGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAAc                                                                                                                                        <  3:445048/149‑1 (MQ=255)
          gTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                                                                                       >  3:85981/1‑149 (MQ=255)
            tGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                     <  1:682202/149‑1 (MQ=255)
                 tCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                                                                                >  1:485935/1‑149 (MQ=255)
                 tCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                                                                                                >  3:68361/1‑149 (MQ=255)
                   cTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACg                                                                                                                              >  4:573621/1‑149 (MQ=255)
                         ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGa                                                                                                                        >  2:507556/1‑149 (MQ=255)
                         ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGa                                                                                                                        >  2:263089/1‑149 (MQ=255)
                          tCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAt                                                                                                                       <  2:576592/149‑1 (MQ=255)
                                aGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCa                                                                                                                 <  2:426072/149‑1 (MQ=255)
                                   tgctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATAc                                                                                                              >  4:39429/1‑149 (MQ=255)
                                      tgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaaa                                                                                                           <  2:485935/149‑1 (MQ=255)
                                      tgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaaa                                                                                                           <  4:68361/149‑1 (MQ=255)
                                          cGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                                                                                                       <  2:352384/149‑1 (MQ=255)
                                            gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                     <  1:670056/149‑1 (MQ=255)
                                                   gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                     >  3:398668/1‑110 (MQ=255)
                                                   gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                     <  4:398668/110‑1 (MQ=255)
                                                     cTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                                                                            <  1:408002/149‑1 (MQ=255)
                                                      tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                                                                           >  1:497716/1‑149 (MQ=255)
                                                      tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                                                                           >  4:543774/1‑149 (MQ=255)
                                                         tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGa                                                                                        >  4:389027/1‑149 (MQ=255)
                                                          aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATGAACGAAGCGAt                                                                                       <  3:349082/149‑1 (MQ=255)
                                                             gCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGa                                                                                    >  4:153379/1‑149 (MQ=255)
                                                               aTACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAATCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                                                                  >  2:465962/1‑149 (MQ=255)
                                                               aTACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                                                                  >  3:189577/1‑149 (MQ=255)
                                                               aTACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                                                                  >  4:467462/1‑149 (MQ=255)
                                                                tACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACt                                                                                 <  3:323257/149‑1 (MQ=255)
                                                                   cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATgcgcag                                                                                                                                               >  1:166527/1‑82 (MQ=255)
                                                                   cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATgcgcag                                                                                                                                               <  2:166527/84‑3 (MQ=255)
                                                                             tgatgaAGCCATAGCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                                                                    <  4:437812/149‑1 (MQ=255)
                                                                             tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                                                                    <  1:519674/149‑1 (MQ=255)
                                                                             tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                                                                    <  3:467462/149‑1 (MQ=255)
                                                                             tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                                                                    >  3:453467/1‑149 (MQ=255)
                                                                             tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                                                                    <  3:39429/149‑1 (MQ=255)
                                                                               atgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTgg                                                                               >  4:519108/1‑136 (MQ=255)
                                                                               atgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTgg                                                                               <  3:519108/136‑1 (MQ=255)
                                                                                tgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgc                                                                 >  1:687399/1‑149 (MQ=255)
                                                                                 gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgcg                                                                >  2:332247/1‑149 (MQ=255)
                                                                                 gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGgcg                                                                >  1:658750/1‑149 (MQ=255)
                                                                                 gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTa                                                                                                                                             <  3:373152/72‑1 (MQ=39)
                                                                                 gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTa                                                                                                                                             >  4:373152/1‑72 (MQ=255)
                                                                                   aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACt                                                                                 >  4:661037/1‑130 (MQ=255)
                                                                                   aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACt                                                                                 <  3:661037/130‑1 (MQ=255)
                                                                                     ccATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCaga                                                            <  3:242914/149‑1 (MQ=255)
                                                                                         aCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGAc                                                        >  2:330195/1‑149 (MQ=255)
                                                                                           cTGTCCCGCCTGCTGGAAGAACTTGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAg                                                      <  1:389817/149‑1 (MQ=255)
                                                                                           cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAg                                                      >  3:259699/1‑149 (MQ=255)
                                                                                           cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAg                                                      >  2:22777/1‑149 (MQ=255)
                                                                                           cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAg                                                      >  3:343450/1‑149 (MQ=255)
                                                                                              tCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCa                                                   >  2:123280/1‑149 (MQ=255)
                                                                                              tCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCACAGAAGACAGCCa                                                   <  3:153379/149‑1 (MQ=255)
                                                                                                ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATa                                                 <  1:123280/149‑1 (MQ=255)
                                                                                                  gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAg                                               >  3:411812/1‑149 (MQ=255)
                                                                                                  gCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAg                                               >  1:629428/1‑149 (MQ=255)
                                                                                                        tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                         >  2:543938/1‑149 (MQ=255)
                                                                                                         ggAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCa                                        >  3:296489/1‑149 (MQ=255)
                                                                                                         ggAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCa                                        >  1:339921/1‑149 (MQ=255)
                                                                                                          gaagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCaa                                       >  2:70518/1‑149 (MQ=255)
                                                                                                              aaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGTGCAGAAGACAGCCATAAGTTCGTCAACGGc                                   <  2:339921/149‑1 (MQ=255)
                                                                                                               aCTGGGACAGGTAGAAAAAGCAGTACTGCGCAGTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTGGTCAACGGCg                                  >  4:441136/1‑149 (MQ=255)
                                                                                                                  gggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTAc                               >  2:602393/1‑149 (MQ=255)
                                                                                                                     aCAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGTGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCg                            >  4:140819/1‑149 (MQ=255)
                                                                                                                     aCAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCg                            >  2:667136/1‑149 (MQ=255)
                                                                                                                      cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa                           <  2:629428/149‑1 (MQ=255)
                                                                                                                         gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATaa                        >  4:55330/1‑149 (MQ=255)
                                                                                                                         gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATaa                        <  3:445054/149‑1 (MQ=255)
                                                                                                                            gAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagc                     >  3:33802/1‑149 (MQ=255)
                                                                                                                             aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                    >  1:546656/1‑149 (MQ=255)
                                                                                                                             aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                    >  1:75514/1‑149 (MQ=255)
                                                                                                                              aaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcag                   >  2:295004/1‑149 (MQ=255)
                                                                                                                               aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAGCGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcagc                  >  4:93554/1‑149 (MQ=255)
                                                                                                                               aaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagcagc                  >  4:491179/1‑149 (MQ=255)
                                                                                                                                aaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAgcagca                 >  1:628787/1‑149 (MQ=255)
                                                                                                                                 aGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAc                >  3:332181/1‑149 (MQ=255)
                                                                                                                                  ttcgTACGGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAcc               >  2:373096/4‑149 (MQ=255)
                                                                                                                                     gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgt            <  1:465962/149‑1 (MQ=255)
                                                                                                                                     gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgt            >  4:671376/1‑149 (MQ=255)
                                                                                                                                     gTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgt            <  4:77502/149‑1 (MQ=255)
                                                                                                                                         tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt        >  3:268828/1‑149 (MQ=255)
                                                                                                                                           cgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCg      >  2:442783/1‑149 (MQ=255)
                                                                                                                                             cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc    <  1:507556/149‑1 (MQ=255)
                                                                                                                                             cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc    <  1:70518/149‑1 (MQ=255)
                                                                                                                                              aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc   >  4:281679/1‑149 (MQ=255)
                                                                                                                                               ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCt  >  1:181330/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
TGTAAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCT  >  NC_000913/435253‑435544

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: