Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,014,867 A→G 44.5% I547V (ATA→GTA)  fliF → flagellar basal‑body MS‑ring and collar protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,014,8670AG44.5% 31.0 / 190.3 146I547V (ATA→GTA) fliFflagellar basal‑body MS‑ring and collar protein
Reads supporting (aligned to +/- strand):  ref base A (32/49);  new base G (19/46);  total (51/95)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-02

Only 100 of 160 total aligned reads displayed.

CGCGCTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACC  >  NC_000913/2014768‑2014964
                                                                                                   |                                                                                                 
cgcgcTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGa                                                                                                   >  1:2827724/1‑100 (MQ=255)
 gcgcTAACCAACGTCTGGGGGCAGAAGTCATGAGCTAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGat                                                                                                  <  2:1504813/100‑1 (MQ=255)
 gcgcTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGat                                                                                                  <  2:3715110/100‑1 (MQ=255)
  cgcTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCGGTGGata                                                                                                 <  2:2241260/100‑1 (MQ=255)
   gcTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGgtaa                                                                                                >  2:3077687/1‑100 (MQ=255)
    cTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGataaa                                                                                               >  2:1228953/1‑100 (MQ=255)
      aaCCAACGTCTGGGGGCAGAAGTCATGAGTCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCATGGTGGCGCTGGTCATTCGCCAGTGGataaata                                                                                             <  1:2570872/100‑1 (MQ=255)
        ccAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAAc                                                                                           >  2:2321182/1‑100 (MQ=255)
            cGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATc                                                                                       <  2:4154273/100‑1 (MQ=255)
            cGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATc                                                                                       <  2:392701/100‑1 (MQ=255)
             gTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCa                                                                                      <  2:3134308/100‑1 (MQ=255)
              tCTGGGGGCAGAAGTCATGAGCCGGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCAt                                                                                     <  1:1813288/100‑1 (MQ=255)
              tCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCAt                                                                                     >  2:902903/1‑100 (MQ=255)
               cTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATg                                                                                    <  2:1044813/100‑1 (MQ=255)
                 gggggCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGan                                                                                  <  1:2685875/100‑2 (MQ=255)
                 gggggCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAg                                                                                  >  2:4266065/1‑100 (MQ=255)
                      cAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAAcc                                                                             <  1:3713944/100‑1 (MQ=255)
                      cAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAAcc                                                                             <  1:4000493/100‑1 (MQ=255)
                        gAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTg                                                                           >  1:1346572/1‑100 (MQ=255)
                         aaGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGa                                                                          >  1:3882564/1‑100 (MQ=255)
                           gTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACa                                                                        >  2:2791506/1‑100 (MQ=255)
                             cATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAgg                                                                      >  1:1604534/1‑100 (MQ=255)
                              aTGAGCCAGCGTATCCTTTAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGc                                                                     <  2:14426/100‑1 (MQ=255)
                              aTGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGc                                                                     <  2:3994588/100‑1 (MQ=255)
                              aTGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGc                                                                     <  2:419806/100‑1 (MQ=255)
                                gAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGTCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGAcn                                                                        <  1:2121312/95‑2 (MQ=255)
                                gAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCAc                                                                   <  2:1774245/100‑1 (MQ=255)
                                 aGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCAcc                                                                  >  1:2166604/1‑100 (MQ=255)
                                 aGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTCGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCAcc                                                                  <  2:1431389/100‑1 (MQ=255)
                                 aGCCAGCATATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCAcc                                                                  >  1:3057655/1‑100 (MQ=255)
                                  gCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGTTCATGAGTAACCTGACAGGCACCg                                                                 <  1:704148/100‑1 (MQ=255)
                                  gCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCAc                                                                   >  1:478726/1‑98 (MQ=255)
                                   ccAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGa                                                                <  1:2827847/100‑1 (MQ=255)
                                   ccAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGa                                                                >  2:2137277/1‑100 (MQ=255)
                                     aGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATa                                                              >  2:4191494/1‑100 (MQ=255)
                                      gCGTAGCCGTGAAACGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATaa                                                             <  1:2906707/100‑1 (MQ=255)
                                        gTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAGATAACGATCATGAGCAACCTGACAGGCACCGATaaaa                                                           >  1:2889456/1‑100 (MQ=255)
                                        gTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATaaaa                                                           >  2:2526509/1‑100 (MQ=255)
                                         tATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAg                                                          <  2:2822944/100‑1 (MQ=255)
                                         tATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAg                                                          <  1:483865/100‑1 (MQ=255)
                                           tCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCg                                                        <  2:1750506/100‑1 (MQ=255)
                                            ccGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGt                                                       >  1:4118190/1‑100 (MQ=255)
                                            ccGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGt                                                       <  2:2956393/100‑1 (MQ=255)
                                             cGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTc                                                      >  1:669404/1‑100 (MQ=255)
                                             cGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTc                                                      >  1:2713907/1‑100 (MQ=255)
                                              gTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCa                                                     >  1:4095818/1‑100 (MQ=255)
                                              gTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGAGCAGGAGTAACCTGACAGGCACCGATAAAAGCGTCa                                                     >  2:3213889/1‑100 (MQ=255)
                                                 aaaTGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATcc                                                  <  2:231803/100‑1 (MQ=255)
                                                 aaaTGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATcc                                                  <  1:79483/100‑1 (MQ=255)
                                                 aaaTGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATcc                                                  <  2:3329389/100‑1 (MQ=255)
                                                 aaaTGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATcc                                                  <  2:785378/100‑1 (MQ=255)
                                                  aaTGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCct                                                 <  2:1168261/100‑1 (MQ=255)
                                                   aTGTCTGATAACGATCCGCGCGTGGTGGCTCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctg                                                <  1:2561725/100‑1 (MQ=255)
                                                    tGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgc                                               <  2:350215/100‑1 (MQ=255)
                                                     gTCTGATAGCGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgct                                              <  2:1489571/100‑1 (MQ=255)
                                                     gTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgct                                              <  1:2388325/100‑1 (MQ=255)
                                                     gTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgct                                              <  1:1150544/100‑1 (MQ=255)
                                                     gTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgct                                              <  2:1899165/100‑1 (MQ=255)
                                                     gTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACTATAATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgc                                               <  1:1609983/99‑1 (MQ=255)
                                                     gTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgct                                              >  1:2521534/1‑100 (MQ=255)
                                                      tCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgctg                                             <  2:3231818/100‑1 (MQ=255)
                                                      tCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgctg                                             <  1:1024545/100‑1 (MQ=255)
                                                      tCTGATAACGAGCCGCGCGTGGTGGCGCTGGTCAGTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCctgctg                                             <  2:4036408/100‑1 (MQ=255)
                                                       cTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCtga                                            >  2:1963045/1‑100 (MQ=255)
                                                        tGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCtgat                                           >  1:3897244/1‑100 (MQ=255)
                                                         gATAACTATCCGCGCCTGGTGGCCCTCCTCATCCGCCACTGGATAAATAACGATCAAAAATAACCTGACAGGCACCGATAAAAGCGTCATCCTGCtgatt                                          <  1:3555325/100‑2 (MQ=255)
                                                         gATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCtgatg                                          <  2:2734813/100‑1 (MQ=255)
                                                           tAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGAc                                        <  2:4128707/100‑1 (MQ=255)
                                                             aCGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTTGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCa                                      <  2:880964/100‑1 (MQ=255)
                                                             aCGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCa                                      >  1:1123571/1‑100 (MQ=255)
                                                                   cgcgcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGGTCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACAATTGGCg                                <  1:390949/100‑1 (MQ=255)
                                                                   cgcgcgTGGTGCCGCTGGTAATTCGACAGTGGATAAATCACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTAATCCTGCTGATGACCATTGGc                                 >  1:850934/1‑99 (MQ=255)
                                                                    gcgcgTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGa                               <  1:3289299/100‑1 (MQ=255)
                                                                      gcgTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAg                             <  2:2973030/100‑1 (MQ=255)
                                                                      gcgTGGTGGCGCTGGTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAg                             <  1:2935244/100‑1 (MQ=255)
                                                                      gcgTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAg                             <  2:751220/100‑1 (MQ=255)
                                                                         tggtggCGCTGGTCATTTGCCAGTGGGAAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTAATCCTGCTGATGACCATTGGCGAAGAc                           >  1:758712/1‑99 (MQ=255)
                                                                                  tGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCag                 <  2:2472018/100‑1 (MQ=255)
                                                                                   ggTCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCaga                <  2:3833525/100‑1 (MQ=255)
                                                                                     tCATTCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACGGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAgg              <  2:1560627/100‑1 (MQ=255)
                                                                                     tCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAgg              <  1:589463/100‑1 (MQ=255)
                                                                                     tCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAgg              >  1:938000/1‑100 (MQ=255)
                                                                                       aTTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGgtg            <  2:241254/100‑1 (MQ=255)
                                                                                        ttCGCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGgtgt           <  1:2791495/100‑1 (MQ=255)
                                                                                           gCCAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCa        <  1:967920/100‑1 (MQ=255)
                                                                                           gCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCa        >  2:2174897/1‑100 (MQ=255)
                                                                                           gCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCa        <  2:1806546/100‑1 (MQ=255)
                                                                                           gCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCa        <  2:1254922/100‑1 (MQ=255)
                                                                                           gCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCa        >  1:1820978/1‑100 (MQ=255)
                                                                                           gCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCAGCCTGCAGATGACCATTGACGAAGACCGGGCCTCAGAGgtag           >  1:888211/1‑95 (MQ=255)
                                                                                           gCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGGCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTc         >  1:3651323/1‑99 (MQ=255)
                                                                                            ccAGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCaa       <  1:3077675/100‑1 (MQ=255)
                                                                                            ccAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCaa       <  2:3640727/100‑1 (MQ=255)
                                                                                              aGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGc     <  2:1703461/100‑1 (MQ=255)
                                                                                              aGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGc     <  1:2321176/100‑1 (MQ=255)
                                                                                              aGTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGc     <  2:841132/100‑1 (MQ=255)
                                                                                               gTGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGgtgcttgaaca    <  2:1676650/100‑9 (MQ=255)
                                                                                               gTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCa    <  1:1573251/100‑1 (MQ=255)
                                                                                                tGGGTAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAc   >  2:3423084/1‑100 (MQ=255)
                                                                                                 ggATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCAcc  <  1:2526502/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
CGCGCTAACCAACGTCTGGGGGCAGAAGTCATGAGCCAGCGTATCCGTGAAATGTCTGATAACGATCCGCGCGTGGTGGCGCTGGTCATTCGCCAGTGGATAAATAACGATCATGAGTAACCTGACAGGCACCGATAAAAGCGTCATCCTGCTGATGACCATTGGCGAAGACCGGGCGGCAGAGGTGTTCAAGCACC  >  NC_000913/2014768‑2014964

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