Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 435,402 T→A 5.7% L89Q (CTG→CAG)  nusB → transcription antitermination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913435,4020TA 5.7% 437.5 / 13.4 123L89Q (CTG→CAG) nusBtranscription antitermination protein
Reads supporting (aligned to +/- strand):  ref base T (59/57);  new base A (4/3);  total (63/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 124 total aligned reads displayed.

AAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACAAA  >  NC_000913/435256‑435550
                                                                                                                                                  |                                                                                                                                                    
aaaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGt                                                                                                                                                    >  2:266941/1‑149 (MQ=255)
 aaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTa                                                                                                                                                   >  4:457876/1‑149 (MQ=255)
 aaaGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTa                                                                                                                                                   >  4:642494/1‑149 (MQ=255)
     aCGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTActgc                                                                                                                                               <  3:204682/149‑4 (MQ=255)
       gTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACt                                                                                                                                             >  2:476162/1‑149 (MQ=255)
               ccTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                                                                                                                     <  3:635724/149‑1 (MQ=255)
               ccTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                                                                                                                     <  3:420522/149‑1 (MQ=255)
                 tGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGt                                                                                                                                   >  1:260735/1‑149 (MQ=255)
                      ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGACAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGa                                                                                                                              >  1:293939/1‑149 (MQ=255)
                      ttCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGa                                                                                                                              >  3:380332/1‑149 (MQ=255)
                         cGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtgt                                                                                                                           >  1:677361/1‑149 (MQ=255)
                           tGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGc                                                                                                                         <  3:454489/149‑1 (MQ=255)
                                         gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                           <  3:159616/149‑1 (MQ=255)
                                              gCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                                                                                      <  1:266941/149‑1 (MQ=255)
                                                gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                           >  2:232292/1‑110 (MQ=255)
                                                gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg                                                                                                    >  2:260231/1‑149 (MQ=255)
                                                gACTAATACCGCATACATCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGt                                                                                                                                           <  1:232292/110‑1 (MQ=255)
                                                   tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAGGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                                                                                 >  4:522983/1‑149 (MQ=255)
                                                   tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                                                                                 >  2:446315/1‑149 (MQ=255)
                                                   tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCAATGCGCAGTa                                                                                                                                                   <  4:715425/99‑1 (MQ=255)
                                                   tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCAATGCGCAGTa                                                                                                                                                   >  3:715425/1‑99 (MQ=255)
                                                   tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAct                                                                                                                                                 <  3:4153/101‑2 (MQ=255)
                                                   tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACAGCGCAATGCGCAGTAct                                                                                                                                                 >  4:4153/1‑100 (MQ=255)
                                                        ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                                                                                 <  3:514418/144‑1 (MQ=255)
                                                        ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg                                                                                                 >  4:514418/1‑144 (MQ=255)
                                                        ccGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc                                                                                            <  4:380332/149‑1 (MQ=255)
                                                            aTACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAc                                                                                        <  3:394424/149‑1 (MQ=255)
                                                               ccTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTgg                                                                                     <  3:735648/149‑1 (MQ=255)
                                                                cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGc                                                                                    <  4:456366/149‑1 (MQ=255)
                                                                cTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGc                                                                                    <  4:713330/149‑1 (MQ=255)
                                                                 tCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCg                                                                                   >  3:778680/1‑149 (MQ=255)
                                                                 tCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCg                                                                                   >  2:798956/1‑149 (MQ=255)
                                                                 tCGACAGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAGAGTGGCCATTAACGAAGCGATCGAACTGGCg                                                                                   >  1:230624/1‑149 (MQ=255)
                                                                      ggACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATc                                                                              >  2:500169/1‑149 (MQ=255)
                                                                       gACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTg                                                                                                              <  4:75681/116‑1 (MQ=255)
                                                                         cTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGtt                                                                           >  4:590940/1‑149 (MQ=255)
                                                                          tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                                                                                                             >  1:344461/1‑114 (MQ=255)
                                                                          tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                                                                                                             <  2:344461/114‑1 (MQ=255)
                                                                          tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTc                                                                          <  3:551104/149‑1 (MQ=255)
                                                                            atgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa                                                                                                  >  3:384035/1‑123 (MQ=255)
                                                                            atgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg                                                                        >  1:551643/1‑149 (MQ=255)
                                                                            atgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCgg                                                                        >  1:744247/1‑149 (MQ=255)
                                                                              gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                                                                                >  3:771518/1‑139 (MQ=255)
                                                                              gaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGaaa                                                                                <  4:771518/139‑1 (MQ=255)
                                                                                   cATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaa                                                                 <  3:372846/149‑1 (MQ=255)
                                                                                    aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTTCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaag                                                                >  1:697423/1‑149 (MQ=255)
                                                                                     tACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaaga                                                               >  2:387970/1‑149 (MQ=255)
                                                                                     tACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCagaaga                                                               <  2:445774/149‑1 (MQ=255)
                                                                                        cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                                                                               <  4:585618/66‑1 (MQ=255)
                                                                                        cTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                                                                                                               >  3:585618/1‑66 (MQ=255)
                                                                                         tGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGc                                                           >  1:83946/1‑149 (MQ=255)
                                                                                         tGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGc                                                           >  3:466062/1‑149 (MQ=255)
                                                                                         tGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAACATTACTGCGCTTTGCGCTGTACGAACTGTCTAAACGTAACGATTTGCCATACAAAGTTGCCATTAACGAAGCGATCGAACTAGCGAAAACGTTCGGCGCAGAAGACAGc                                                           >  4:257411/1‑149 (MQ=255)
                                                                                             ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                           <  4:424547/97‑1 (MQ=255)
                                                                                             ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                                                                                                           >  3:424547/1‑97 (MQ=255)
                                                                                             ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACt                                                                                       <  4:341004/117‑1 (MQ=255)
                                                                                             ccGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATa                                                       >  3:625259/1‑149 (MQ=255)
                                                                                                     tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                               >  3:110091/1‑149 (MQ=255)
                                                                                                     tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                               >  1:572478/1‑149 (MQ=255)
                                                                                                     tgGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTc                                               >  1:379247/1‑149 (MQ=255)
                                                                                                       gaagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                                                                                      >  2:100479/1‑92 (MQ=255)
                                                                                                       gaagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTaa                                                                                                      <  1:100479/92‑1 (MQ=255)
                                                                                                        aagaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAAc                                            <  4:778680/149‑1 (MQ=255)
                                                                                                         agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                           <  3:703198/149‑1 (MQ=255)
                                                                                                         agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                           <  3:331348/149‑1 (MQ=255)
                                                                                                         agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                           >  4:331348/1‑149 (MQ=255)
                                                                                                         agaaCTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACg                                           >  4:703198/1‑149 (MQ=255)
                                                                                                              tGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                                      <  3:409571/133‑1 (MQ=255)
                                                                                                              tGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATaa                                                      >  4:409571/1‑133 (MQ=255)
                                                                                                              tGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTa                                      >  1:552311/1‑149 (MQ=255)
                                                                                                                ggACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACt                                    >  1:698479/1‑149 (MQ=255)
                                                                                                                 gACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTc                                   <  4:166782/149‑1 (MQ=255)
                                                                                                                   cAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGa                                 <  1:625949/149‑1 (MQ=255)
                                                                                                                     ggTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATa                               <  2:28912/149‑1 (MQ=255)
                                                                                                                      gTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATaa                              >  2:622106/1‑149 (MQ=255)
                                                                                                                        aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                             <  4:710170/116‑1 (MQ=255)
                                                                                                                        aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACa                                                             >  3:710170/1‑116 (MQ=255)
                                                                                                                        aGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAag                            <  2:697423/149‑1 (MQ=255)
                                                                                                                         gAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagc                           >  2:777802/1‑149 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          <  2:744247/149‑1 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          <  2:552311/149‑1 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          >  2:738396/1‑149 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          <  2:379247/149‑1 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          <  3:120835/149‑1 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          <  2:551643/149‑1 (MQ=255)
                                                                                                                          aaaaaGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAagca                          <  2:572478/149‑1 (MQ=255)
                                                                                                                               ttcgTACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCAcc                     >  4:784555/4‑149 (MQ=255)
                                                                                                                                cAGTACTGCGCATTGCGCTGTACGATCTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCt                    <  3:522983/149‑1 (MQ=255)
                                                                                                                                   tACAGCGCAATGCGCAGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCtgtg                 >  4:413715/1‑149 (MQ=255)
                                                                                                                                      tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt              >  2:108268/1‑149 (MQ=255)
                                                                                                                                      tGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGAtt              >  3:303087/1‑149 (MQ=255)
                                                                                                                                       gcgcATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTc             >  1:362803/1‑149 (MQ=255)
                                                                                                                                          cATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGcc          <  2:579534/149‑1 (MQ=255)
                                                                                                                                           aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc         <  4:110091/149‑1 (MQ=255)
                                                                                                                                           aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc         >  3:97583/1‑149 (MQ=255)
                                                                                                                                           aTTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGccc         <  3:590940/149‑1 (MQ=255)
                                                                                                                                            ttGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCt        >  2:239176/1‑149 (MQ=255)
                                                                                                                                              gcgcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTaa      <  1:777802/149‑1 (MQ=255)
                                                                                                                                                gcTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACa    >  2:24638/1‑149 (MQ=255)
                                                                                                                                                  tGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGCTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACaaa  >  4:620807/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAAAGACGTTGACGTCCTGTACTTCCGTGAGCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAACTGGCGAAATCGTTCGGCGCAGAAGACAGCCATAAGTTCGTCAACGGCGTACTCGATAAAGCAGCACCTGTGATTCGCCCTAACAAA  >  NC_000913/435256‑435550

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: