Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,131,310 A→T 5.4% N98I (AAT→ATT)  flgB → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,131,3100AT 5.4% 318.7 / 7.7 92N98I (AAT→ATT) flgBflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (46/41);  new base T (2/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.37e-01

Only 100 of 103 total aligned reads displayed.

CTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTT  >  NC_000913/1131162‑1131458
                                                                                                                                                    |                                                                                                                                                    
cTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTaa                                                                                                                                                      >  1:562275/1‑149 (MQ=255)
cTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTaa                                                                                                                                                      >  2:308629/1‑149 (MQ=255)
 ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAAt                                                                                                                                                     >  1:441757/1‑149 (MQ=255)
 ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAAt                                                                                                                                                     >  4:393258/1‑149 (MQ=255)
 ttAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCACAACTGCAATACAGTATTCCGGACCAGCATTCGCTTGACGGTAAt                                                                                                                                                     >  2:24070/1‑149 (MQ=255)
   aaaaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                                                                                                                                   <  3:440569/149‑1 (MQ=255)
   aaaaaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                                                                                                                                   >  4:10553/1‑149 (MQ=255)
      aaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGt                                                                                                                                                >  2:92318/1‑149 (MQ=255)
         gTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaa                                                                                                                                             <  2:315883/149‑3 (MQ=255)
         gTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGa                                                                                                                                             <  3:516775/149‑1 (MQ=255)
         gTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGa                                                                                                                                             <  1:327233/149‑1 (MQ=255)
          tCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaag                                                                                                                                            <  3:574111/149‑4 (MQ=255)
            aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGatat                                                                                                                                          <  2:11240/149‑1 (MQ=255)
            aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGatat                                                                                                                                          >  4:242635/1‑149 (MQ=255)
              gCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATgg                                                                                                                                        <  1:124486/149‑1 (MQ=255)
               cAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGa                                                                                                                                       >  1:218317/1‑149 (MQ=255)
                aaCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGAt                                                                                                                                      <  2:481449/149‑1 (MQ=255)
                 aCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATc                                                                                                                                     >  1:561629/1‑149 (MQ=255)
                    tGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATcgcg                                                                                                                                  <  3:287884/149‑1 (MQ=255)
                      gACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                                                                                                                                <  1:597755/149‑1 (MQ=255)
                           agtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc                                                                                                                           >  4:332071/3‑149 (MQ=255)
                            aTGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAcc                                                                                                                          >  3:559233/1‑149 (MQ=255)
                                aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGt                                                                                                                      >  2:358046/1‑149 (MQ=255)
                                aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGt                                                                                                                      >  4:279250/1‑149 (MQ=255)
                                 aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt                                                                                                                     >  4:236693/1‑149 (MQ=255)
                                 aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt                                                                                                                     >  2:308737/1‑149 (MQ=255)
                                  ccAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt                                                                                                                    >  4:153512/1‑149 (MQ=255)
                                    aGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATTGATCGCGAACGCACCCAGtttac                                                                                                                  >  1:14073/1‑147 (MQ=255)
                                    aGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGc                                                                                                                  <  3:393258/149‑1 (MQ=255)
                                    aGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGc                                                                                                                  <  3:10553/149‑1 (MQ=255)
                                          gTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATaa                                                                                                            >  2:143050/1‑149 (MQ=255)
                                          gTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATaa                                                                                                            >  3:153836/1‑149 (MQ=255)
                                               aCTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGcc                                                                                                       >  1:191210/1‑149 (MQ=255)
                                               aCTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGcc                                                                                                       <  3:314048/149‑1 (MQ=255)
                                                 tGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTg                                                                                                     >  2:246420/1‑149 (MQ=255)
                                                    cGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCaa                                                                                                  <  3:279250/149‑1 (MQ=255)
                                                       tGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATAc                                                                                               <  1:358046/149‑1 (MQ=255)
                                                         aCCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCa                                                                                             >  3:182964/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           <  2:581570/149‑1 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           >  1:19759/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           <  3:137059/149‑1 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           >  1:406763/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           >  1:407465/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           >  4:306785/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           >  1:437292/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGa                                                                                           >  4:24865/1‑149 (MQ=255)
                                                           cTCAACGCAACACATTCCGGCGCAGCGGCTGACGCCTCCTACCGCAGAACTGCAATACCCTATCCCGGAGCAGCCTTCGCCTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCCGTTTACCCATAACAGCCTGCAATACCAGa                                                                                           <  4:491191/149‑1 (MQ=255)
                                                            tCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGcc                                                                                                                 >  3:328596/1‑126 (MQ=255)
                                                            tCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGcc                                                                                                                 <  4:328596/126‑1 (MQ=255)
                                                              aaCGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGa                                                                                        <  1:308629/149‑1 (MQ=255)
                                                                   aaCACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCtt                                                                                   >  2:454762/1‑149 (MQ=255)
                                                                    acacATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTa                                                                                  <  1:493058/149‑1 (MQ=255)
                                                                     cacaTTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAg                                                                                 >  3:136499/1‑149 (MQ=255)
                                                                     cacaTTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAg                                                                                 <  3:236693/149‑1 (MQ=255)
                                                                      acaTTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAgc                                                                                >  4:546240/1‑149 (MQ=255)
                                                                      acaTTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAgc                                                                                <  4:136499/149‑1 (MQ=255)
                                                                           ccGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGtt                                                                           <  1:454762/149‑1 (MQ=255)
                                                                               cgcAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGc                                                                       <  4:138528/149‑1 (MQ=255)
                                                                                 cAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCgg                                                                     <  3:204958/149‑1 (MQ=255)
                                                                                          aCGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACGGCCTGCAATAc                                                                                               <  3:333377/114‑1 (MQ=255)
                                                                                          aCGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATAc                                                                                               >  4:333377/1‑114 (MQ=255)
                                                                                           cGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           <  2:19759/149‑1 (MQ=255)
                                                                                           cGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           <  1:152882/149‑1 (MQ=255)
                                                                                           cGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           <  3:15349/149‑1 (MQ=255)
                                                                                            gCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAg                                                          >  4:270156/1‑149 (MQ=255)
                                                                                              ctcctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGc                                                        >  1:550673/1‑149 (MQ=255)
                                                                                               tcctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCa                                                       <  2:428702/149‑1 (MQ=255)
                                                                                               tcctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCa                                                       <  2:184938/149‑1 (MQ=255)
                                                                                                 ctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatg                                                     <  1:28225/149‑1 (MQ=255)
                                                                                                  tACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatga                                                    <  3:532079/149‑1 (MQ=255)
                                                                                                          aaCTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           >  2:382066/1‑134 (MQ=255)
                                                                                                          aaCTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           >  2:90176/1‑134 (MQ=255)
                                                                                                          aaCTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           <  1:90176/134‑1 (MQ=255)
                                                                                                          aaCTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                           <  1:382066/134‑1 (MQ=255)
                                                                                                           aCTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGttt                                           >  2:548189/1‑149 (MQ=255)
                                                                                                              gCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTAc                                        >  1:260375/1‑149 (MQ=255)
                                                                                                              gCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTAc                                        <  2:411279/149‑1 (MQ=255)
                                                                                                                aaTACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACag                                      >  4:96906/1‑149 (MQ=255)
                                                                                                                aaTACCGTATTCCGGACCAGCCTTCGCTGGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGATTTACag                                      >  3:503192/1‑149 (MQ=255)
                                                                                                                       tATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc                                                                                                                           <  1:247247/57‑1 (MQ=255)
                                                                                                                       tATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc                                                                                                                           >  2:247247/1‑57 (MQ=255)
                                                                                                                       tATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGaa                               >  3:166512/1‑149 (MQ=255)
                                                                                                                         ttCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAAt                             <  2:407465/149‑1 (MQ=255)
                                                                                                                               aCCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACgg                       >  2:220658/1‑149 (MQ=255)
                                                                                                                                 cAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGAt                     <  3:497377/149‑1 (MQ=255)
                                                                                                                                   gCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATgg                   >  1:403609/1‑149 (MQ=255)
                                                                                                                                   gCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATgg                   <  4:33360/149‑1 (MQ=255)
                                                                                                                                   gCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAACTAACGGATgg                   <  2:478561/149‑1 (MQ=255)
                                                                                                                                          atcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                         >  3:293937/4‑120 (MQ=255)
                                                                                                                                          atcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                         <  4:293937/117‑1 (MQ=255)
                                                                                                                                           tcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCActgctg           >  3:385486/3‑149 (MQ=255)
                                                                                                                                              aCGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAat        >  1:429577/1‑149 (MQ=255)
                                                                                                                                                   aaTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAtttt   <  2:348050/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGATCGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  3:99303/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  1:429593/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  1:548189/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  2:218317/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  2:260375/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  3:24865/149‑1 (MQ=255)
                                                                                                                                                    aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt  <  3:512056/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTT  >  NC_000913/1131162‑1131458

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