Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,131,401 A→C 5.2% K128N (AAA→AAC 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,4010AC 5.2% 232.1 / 4.5 77K128N (AAA→AACflgBflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base C (2/2);  total (35/42)
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.67e-01

CGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGTGACCGGTCCC  >  NC_000913/1131253‑1131548
                                                                                                                                                    |                                                                                                                                                   
cGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCaaa                                                                                                                                                     <  3:36594/149‑1 (MQ=255)
 gCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCGGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAg                                                                                                                                                    >  8:264184/1‑149 (MQ=255)
 gCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAg                                                                                                                                                    >  2:89006/1‑149 (MQ=255)
 gCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAg                                                                                                                                                    <  1:89006/149‑1 (MQ=255)
      ctACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCa                                                                                                                                                 <  8:216286/147‑1 (MQ=255)
         ccGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatgat                                                                                                                                             <  2:368/148‑1 (MQ=255)
         ccGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatgat                                                                                                                                             <  5:3459/148‑1 (MQ=255)
         ccGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatga                                                                                                                                              <  5:341769/147‑1 (MQ=255)
         ccGCAGAACTGCAAAACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATAGCGAACGCACCCAGTTTGCCGATAACACCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCatgat                                                                                                                                             >  7:284463/1‑148 (MQ=255)
            cAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAAc                                                                                                                                         >  8:328330/1‑149 (MQ=255)
               aaCTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtt                                                                                                                                      >  8:385541/1‑149 (MQ=255)
               aaCTGCAATACCGTATTCCGGACCAGCCGTCGCGTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtt                                                                                                                                      <  3:300703/149‑1 (MQ=255)
                 cTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGtttt                                                                                                                                    <  5:156153/149‑1 (MQ=255)
                  tGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                                                                                                                   >  7:380923/1‑149 (MQ=255)
                    cactACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTGGATATGTATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGTTGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                                                                                                                   <  8:284463/144‑1 (MQ=255)
                    cAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                                                                                                                   >  1:209412/1‑147 (MQ=255)
                    cAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTAc                                                                                                                                  <  2:197683/148‑1 (MQ=255)
                     aaTACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTa                                                                                                                                   >  1:200420/1‑146 (MQ=255)
                     aaTACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACag                                                                                                                                <  5:296805/149‑1 (MQ=255)
                       tACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACagag                                                                                                                              >  3:279314/1‑149 (MQ=255)
                         ccGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATACCAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCg                                                                                                                            <  2:256289/149‑1 (MQ=255)
                                  ggACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAAc                                                                                                                   <  5:151979/149‑1 (MQ=255)
                                   gACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACg                                                                                                                  <  1:91830/149‑1 (MQ=255)
                                         ccTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATgg                                                                                                             >  1:1385/1‑148 (MQ=255)
                                         ccTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATgg                                                                                                             >  4:276463/1‑148 (MQ=255)
                                           ccatatcgACGGTATTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCAc                                                                                                          >  7:324582/8‑149 (MQ=255)
                                             cGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACGGAGCGGAAATTAACGGCTGGCAc                                                                                                          <  6:308465/147‑1 (MQ=255)
                                              gCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAACGGCATGATGAACGTTTTACAGAGCGGAAATTAACTGATGGCActgc                                                                                                       >  2:271125/1‑149 (MQ=255)
                                              gCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAACGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCActgc                                                                                                       >  2:265681/1‑149 (MQ=255)
                                              gCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCActgc                                                                                                       >  3:216214/1‑149 (MQ=255)
                                               cTTGATGGTAATACCGTCGATATGGATCGCGAACTCACCCATTTTGCCGATAACTGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCTTGATGAACGTTTTa                                                                                                                                   <  6:170034/120‑1 (MQ=255)
                                               cTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATCAACGTTTTa                                                                                                                                   >  5:170034/1‑120 (MQ=255)
                                                       tAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGTAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttt                                                                                              >  6:87297/1‑149 (MQ=255)
                                                         ttACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt                                                                                            <  8:324582/148‑1 (MQ=255)
                                                         aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt                                                                                            <  8:380923/149‑1 (MQ=255)
                                                         aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt                                                                                            <  4:216214/149‑1 (MQ=255)
                                                         aTACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAtttt                                                                                             <  3:223027/148‑1 (MQ=255)
                                                            ccGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt                                                                                            <  2:369151/146‑1 (MQ=255)
                                                            ccGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATAttttt                                                                                            >  1:369151/1‑146 (MQ=255)
                                                            ccGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGa                                                                                          <  4:224163/148‑1 (MQ=255)
                                                            ccGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGa                                                                                          <  7:385541/148‑1 (MQ=255)
                                                            ccGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGa                                                                                          >  4:338959/1‑148 (MQ=255)
                                                               tCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAAAACCAGAAGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTATTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATc                                                                                      >  1:203142/1‑149 (MQ=255)
                                                                 gATATGGATCTCGAACGCACCCAGTTTCCCGATAACAGCCTGCAATACCAGATGAGCATGAGCGCGTTGAGCGGGCAAATCAAAGGCATCATGAACGTTCAACAGAGCGGAAATTAACGGATGGCACTGCTGAATATATTTGataactc                                                                                    >  5:267704/1‑145 (MQ=255)
                                                                       gATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGt                                                                               <  3:314577/148‑1 (MQ=255)
                                                                       gATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTc                                                                              >  6:159116/1‑149 (MQ=255)
                                                                         tCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCgg                                                                            <  2:200420/149‑1 (MQ=255)
                                                                             gAACGCACCCTGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGtt                                                                        <  7:264184/149‑1 (MQ=255)
                                                                             gAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGTAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGtt                                                                        <  5:87297/149‑1 (MQ=255)
                                                                                ggCACCCAGTTTGGCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAAc                                                                     <  2:203142/148‑1 (MQ=255)
                                                                                   aCCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGc                                                                  >  2:290985/1‑149 (MQ=255)
                                                                                            gCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTccc                                                          <  8:18754/148‑1 (MQ=255)
                                                                                              cGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTACTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAgc                                                       >  3:11375/1‑149 (MQ=255)
                                                                                               gATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCCGTCCCAgcg                                                      >  6:182058/1‑149 (MQ=255)
                                                                                                   aCAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCt                                                   >  5:317027/1‑148 (MQ=255)
                                                                                                    cAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTg                                                  >  2:388205/1‑148 (MQ=255)
                                                                                                      gCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGc                                                                                    <  6:103541/112‑1 (MQ=255)
                                                                                                      gCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGc                                                                                    >  5:103541/1‑112 (MQ=255)
                                                                                                        cTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAACGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGt                                             <  1:265681/149‑1 (MQ=255)
                                                                                                        cTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAACGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGt                                             <  1:271125/149‑1 (MQ=255)
                                                                                                         tGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTaa                                                                      <  6:348106/123‑1 (MQ=255)
                                                                                                           cAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTg                                            <  5:159116/147‑1 (MQ=255)
                                                                                                              tACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGTCATGATGCACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTggcggc                                       <  1:377489/149‑1 (MQ=255)
                                                                                                              tACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTc                                                                              >  1:19591/1‑110 (MQ=255)
                                                                                                              tACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTc                                                                              <  2:19591/110‑1 (MQ=255)
                                                                                                              tACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTggcggc                                       <  3:276463/149‑1 (MQ=255)
                                                                                                                ccAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCa                                     <  2:1385/149‑1 (MQ=255)
                                                                                                                  aGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTc                                                            >  6:2685/1‑124 (MQ=255)
                                                                                                                  aGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTc                                                            <  5:2685/124‑1 (MQ=255)
                                                                                                                        gCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTg                             <  1:290985/149‑1 (MQ=255)
                                                                                                                                  gTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATggt                    >  3:117041/1‑146 (MQ=255)
                                                                                                                                   ttGAGCTGGCAAATAATAGTTATGATGAACGTTTTTCTGAGGGTATATTAAGGGTTGTCACTGCTGAATTTTTTTGATATCGCCGGGTAGGAGTTATCTGCCCAGTACCAGCGACTTATCGTGGCGGGCAGTAATCTGGCGAATGCTGa                  <  4:279314/149‑1 (MQ=255)
                                                                                                                                   ttGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGa                  >  3:230634/1‑149 (MQ=255)
                                                                                                                                     gAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGAt                 <  2:209412/148‑1 (MQ=255)
                                                                                                                                         gggCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGt            >  5:196629/1‑149 (MQ=255)
                                                                                                                                         gggCAAATCAAAGGCATGATGAACGTGTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGGGGCGGCCAGTAATCTGGCGAATGCTGATAGCGt            <  2:127848/149‑1 (MQ=255)
                                                                                                                                           gCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGGCAGTAATCTGGCGAATGCTGATAGCGTGa          >  8:374538/1‑149 (MQ=255)
                                                                                                                                                tCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGTGACCGGt     >  8:141695/1‑149 (MQ=255)
                                                                                                                                                  aaaGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGTGACCGGTcc   >  1:322212/1‑149 (MQ=255)
                                                                                                                                                   aaGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGTGACCGGTccc  >  4:288024/1‑149 (MQ=255)
                                                                                                                                                   aaGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGTGACCGGTccc  >  8:341827/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCGATAACAGCCTGCAATACCAGATGAGCCTTAGCGCGTTGAGCGGGCAAATCAAAGGCATGATGAACGTTTTACAGAGCGGAAATTAACGGATGGCACTGCTGAATATTTTTGATATCGCCGGGTCGGCGTTAACTGCCCAGTCCCAGCGCCTGAACGTGGCGGCCAGTAATCTGGCGAATGCTGATAGCGTGACCGGTCCC  >  NC_000913/1131253‑1131548

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