Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,965,065 G→A 6.0% A189V (GCA→GTA)  rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,0650GA 6.0% 235.5 / 3.3 83A189V (GCA→GTA) rhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (40/37);  new base A (3/2);  total (44/39)
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.08e-01

GGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGC  >  NC_000913/3964918‑3965212
                                                                                                                                                   |                                                                                                                                                   
gggTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTTGCGGAGAACAGCATGTTGAGGCGCTGGTTTGc                                                                                                                                                    >  2:697494/1‑149 (MQ=255)
 ggTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCa                                                                                                                                                   <  2:484319/149‑1 (MQ=255)
  gTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCATTATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAg                                                                                                                                                  <  2:353442/149‑1 (MQ=255)
  gTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAg                                                                                                                                                  >  2:568789/1‑149 (MQ=255)
      aaaaGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg                                                                                                                                              <  1:516701/149‑1 (MQ=255)
      aaaaGCTCTTCTTTAATACGGTGGCCCGGTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg                                                                                                                                              <  1:446410/149‑1 (MQ=255)
      aaaaGCTCTTCTTTAATAAGGGGGCCCGTTTTCTGTTCCGGATCCACTTCAATATCTTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGAGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg                                                                                                                                              <  2:474156/149‑1 (MQ=255)
       aaaGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCTTGTTGAGGCGCTGGTTTGCAggcggc                                                                                                                                             >  2:968639/1‑149 (MQ=255)
        aaGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTACCACTTCAATATAATCGGCATTGAACATCTGCTCAAACGACAGTTCAAGTACCCGGTACGAAAGCGAGGCGGAGAACAGCATCTTGAGGCGCTGGTTTGCAggagcct                                                                                                                                            >  1:745234/1‑144 (MQ=255)
          gCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGAATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGAATGTTGAGGCGCTGGTTTGCAGGCGGCATa                                                                                                                                          >  1:827287/1‑149 (MQ=255)
               tcttTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcg                                                                                                                                     >  1:735718/1‑149 (MQ=255)
                cttTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATAcggcgg                                                                                                                                    <  1:1114864/149‑1 (MQ=255)
                   tAATACGGGGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTGGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAc                                                                                                                                 <  2:871960/149‑1 (MQ=255)
                    aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGaccc                                                                                                                                >  2:687823/1‑146 (MQ=255)
                    aaTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAca                                                                                                                                <  2:421623/149‑1 (MQ=255)
                      tACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGATcagc                                                                                                                              <  1:295948/149‑1 (MQ=255)
                      tACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGATTGCAGGCGGCATACGGCGGAAcagc                                                                                                                              <  1:163364/149‑1 (MQ=255)
                         ggTGGCCCGTTTTCTGTTCCGGATCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAAcagccag                                                                                                                           <  1:735620/149‑1 (MQ=255)
                             gCCCGTTTTCTGATCCTGTTCCACTTCAATAAATTCGGCATTGTTCATCTGCTCTAACCCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGACGAGATGGTTTACATGCGGCAAACCGAGTAACAGCCAGAGGa                                                                                                                       >  2:3630/1‑149 (MQ=255)
                               ccGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTACAGGCGGCATACGGCGGCACAGCCAGCGGata                                                                                                                     >  2:8990/1‑149 (MQ=255)
                               ccGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTACAGGCGGCATACGGCGGAACAGCCAGCGGata                                                                                                                     >  2:4023/1‑149 (MQ=255)
                               ccGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGATTGCAGGCGGCATACGGCGGAACAGCCAGCGGaaa                                                                                                                     >  1:350666/1‑147 (MQ=255)
                                cGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGATCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCAACTGGAGGCGCAGGAATCCAGGCGGACAACCGCAGAACAGACAGCgcaaat                                                                                                                    >  1:809882/1‑144 (MQ=255)
                                  ttttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACTTAACCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGCAACAGCCAGCGGATctct                                                                                                                  >  2:993783/1‑149 (MQ=255)
                                  ttttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCt                                                                                                                  >  1:330073/1‑149 (MQ=255)
                                   tttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                                                                                                                 <  2:984460/149‑1 (MQ=255)
                                   tttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                                                                                                                 <  1:57720/149‑1 (MQ=255)
                                   tttCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCtt                                                                                                                 <  2:1163847/149‑1 (MQ=255)
                                    ttCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCttt                                                                                                                >  2:626353/1‑149 (MQ=255)
                                    ttCTGTGCCGCTTCCACTTCAATATAATCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCttt                                                                                                                <  2:360259/149‑1 (MQ=255)
                                        gTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGACCAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTtaata                                                                                                            >  1:72757/1‑149 (MQ=255)
                                           ccGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGCACAGCATGTTGAGGCGCTGGATTGCCGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATaaat                                                                                                         >  1:332698/1‑148 (MQ=255)
                                                     ttCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGt                                                                                               <  1:136691/149‑1 (MQ=255)
                                                      tCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGGTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTa                                                                                              <  2:730734/149‑1 (MQ=255)
                                                        aaTATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACa                                                                                            <  1:678118/149‑1 (MQ=255)
                                                         atatatTCGCCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACAt                                                                                           <  1:87208/149‑1 (MQ=255)
                                                            tatTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCg                                                                                        >  1:509795/1‑149 (MQ=255)
                                                             atTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGa                                                                                       >  1:404358/1‑149 (MQ=255)
                                                              ttCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACGCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCATGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGAt                                                                                      >  1:1164847/1‑149 (MQ=255)
                                                                  gCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGc                                                                                  >  1:706348/1‑149 (MQ=255)
                                                                     ttGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGATTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTc                                                                               <  1:607811/149‑1 (MQ=255)
                                                                        ttCATATGCTCGTACGCCAGTTCACGTACCCGGTTCGAAATCGTGGCGGAGAACAGCATTTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTc                                                                            <  1:1054772/149‑1 (MQ=255)
                                                                           aTCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAg                                                                         <  1:961724/149‑1 (MQ=255)
                                                                           aTCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAg                                                                         <  2:330073/149‑1 (MQ=255)
                                                                           aTCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAg                                                                         >  2:425977/1‑149 (MQ=255)
                                                                            tcagCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCACAGATCGTACATGCGATCGGCTTCGTCCAGt                                                                        <  1:21563/146‑1 (MQ=255)
                                                                              tGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTCCATGCGATCGGCTTCGTCCAGTac                                                                      >  2:1023230/1‑149 (MQ=255)
                                                                               gCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTacc                                                                     <  2:53319/149‑1 (MQ=255)
                                                                                cTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTacca                                                                    >  2:602242/1‑149 (MQ=255)
                                                                                 tCGAACGCCAGGTCACGTACCGGGTACGAAAGCGTGGCGGAGAACAGCATGGTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGACCGGCTTCGTCCAGTaccac                                                                   <  2:670717/149‑1 (MQ=255)
                                                                                   gAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTaccacca                                                                 >  2:1086789/1‑149 (MQ=255)
                                                                                   gAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTaccacca                                                                 >  2:37898/1‑149 (MQ=255)
                                                                                    aaCGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTaccaccac                                                                <  2:132513/149‑1 (MQ=255)
                                                                                       gCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTg                                                             >  1:1076944/1‑149 (MQ=255)
                                                                                         cAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGaa                                                           <  1:9011/149‑1 (MQ=255)
                                                                                          aGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAAt                                                          <  2:1143544/149‑1 (MQ=255)
                                                                                               aCGTACCCGGTACGAAAGCGTGGCGGAGAACAGCAAGTTGAGGCGCTGGTTTGCAGGCGGCATACGCGGGAACAGCCAGCGGATTTCATTAATAAAGCCCAGAGCGTACCTGCGATCGGCTTCGACCACTACCACCACCTGAATgtccc                                                     >  1:909550/1‑145 (MQ=255)
                                                                                                 gTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCg                                                   >  1:881957/1‑149 (MQ=255)
                                                                                                    cccGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAgg                                                <  2:910896/149‑1 (MQ=255)
                                                                                                     ccGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATTCGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTAGATGCGATCGGATTCGTCCAGTACCACCACCTGAATGGCACCGAGGt                                               <  1:993783/149‑1 (MQ=255)
                                                                                                     ccGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTCCCTGCGATCGGCTTCGTCCAGTCCCACCACCTGAATGGCACCGCGGt                                               >  2:686856/1‑149 (MQ=255)
                                                                                                       ggTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTa                                             >  1:398546/1‑149 (MQ=255)
                                                                                                        gTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTaa                                            <  1:720391/149‑1 (MQ=255)
                                                                                                             aacgCGTGGCGGAGAACAGCATGTGGAGGCGCTGGGTTGCAGGCGGCATACGGCGGAACAGCCAGCGGAGATCTTTAATAAAGCCCAGATCGTACAGGCGCTCGGCTTCGTCCAGTACCCCCACCGGAATGGCACAGAGGTTAATGTgg                                       >  1:408654/4‑149 (MQ=255)
                                                                                                             aaaGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACCTGGGATCGGCTTCGTCCAGTACCACCGCCTGAATGGAACCGAGGTTAAGGTgg                                       >  2:990094/1‑149 (MQ=255)
                                                                                                             aaaGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTgg                                       >  1:408653/1‑149 (MQ=255)
                                                                                                             aaaGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGAAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGACCAGATCGTACATGCGATCGGCTTCGTCCAGTACAACCACCTGAATGGCACCGAGGTTAATGTgg                                       >  1:1052463/1‑149 (MQ=255)
                                                                                                                   tGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGc                                 >  1:71500/1‑149 (MQ=255)
                                                                                                                   tGGCGGAGAACAGCATGTTGAGGCGCTGGTTTACAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGc                                 <  1:8990/149‑1 (MQ=255)
                                                                                                                   tGGCGGAGAACAGCATGTTGAGGCGCTGGTTTACAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGc                                 <  1:4023/149‑1 (MQ=255)
                                                                                                                        gagaACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGGCCAGATCGTACATGCGATCGGCTTCGTCCAGTAGCACCACCTGAATGGCACCGAGTTTAATGTGGTTCTGCTTGGc                            <  1:3630/149‑1 (MQ=255)
                                                                                                                        gagaACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGc                            >  1:507272/1‑149 (MQ=255)
                                                                                                                            aCAGCATGTTGAGGCGCTGGTTTGCAGGCGGCTTACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAg                        <  1:459272/149‑1 (MQ=255)
                                                                                                                                caagTTGATGCGCGGGGATGCAGGCGGCATACGGCGGAACAGCCAGCGGATAGCTATAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCTCCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCaa                    <  1:697494/146‑1 (MQ=255)
                                                                                                                                 aTGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAAt                   <  2:284956/149‑1 (MQ=255)
                                                                                                                                   gTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCCCCGAGGGTAATGTGGTTCTGCTTGGCGTAGTCAATTa                 >  2:298839/1‑149 (MQ=255)
                                                                                                                                   gTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTa                 >  2:343115/1‑149 (MQ=255)
                                                                                                                                      gAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAAc              >  1:615824/1‑149 (MQ=255)
                                                                                                                                      gAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTACGTCCAGTACCACCACCTGCATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAAc              >  2:965851/1‑149 (MQ=255)
                                                                                                                                        ggCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTAAATAAAGCCCAGATCGTACATGAGATCGGCTTCGACCAGTACCACCACCTGAATGGCACCGATGTTAATGTTGTTCTCCTTGGCGGAGTCAATTAAAcac            >  2:1003122/1‑147 (MQ=255)
                                                                                                                                              ggTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTAGCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAACCGCCCCgtg      >  2:778205/1‑149 (MQ=255)
                                                                                                                                                 ttGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCgtggtg   >  1:490896/1‑149 (MQ=255)
                                                                                                                                                  tGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGc  <  1:602242/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCATACGGCGGAACAGCCAGCGGATATCTTTAATAAAGCCCAGATCGTACATGCGATCGGCTTCGTCCAGTACCACCACCTGAATGGCACCGAGGTTAATGTGGTTCTGCTTGGCGTAGTCAATTAAACGCCCCGTGGTGC  >  NC_000913/3964918‑3965212

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