Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,151,112 A→C D44E (GAT→GAG)  exbD ← membrane spanning protein in TonB‑ExbB‑ExbD complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,151,1120AC93.6% 481.5 / 19.2 156D44E (GAT→GAG) exbDmembrane spanning protein in TonB‑ExbB‑ExbD complex
Reads supporting (aligned to +/- strand):  ref base A (5/5);  new base C (68/78);  total (73/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 161 total aligned reads displayed.

CGATAAACATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCG  >  NC_000913/3151015‑3151208
                                                                                                 |                                                                                                
cGATAAACATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTc                                                                                                <  1:3781042/100‑1 (MQ=255)
 gATAAACATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACctct                                                                                               >  2:2780388/1‑100 (MQ=255)
    aaaCATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTAcc                                                                                            <  1:242067/100‑1 (MQ=255)
     aaCATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGCGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCg                                                                                           >  2:2972115/1‑100 (MQ=255)
      aCATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGt                                                                                          >  2:533735/1‑100 (MQ=255)
       cATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGTTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTc                                                                                         >  2:3624830/1‑100 (MQ=255)
       cATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTc                                                                                         >  2:872585/1‑100 (MQ=255)
       cATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAAGCAAGTTCACCTTCACATCTACCGTc                                                                                         <  1:767711/100‑1 (MQ=255)
        aTCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCg                                                                                        <  2:2528745/100‑1 (MQ=255)
        aTCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCg                                                                                        <  1:3515113/100‑1 (MQ=255)
        aTCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTGCCGGGCGCGGCTGCGGTGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCg                                                                                        <  1:57438/100‑1 (MQ=255)
         tCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGc                                                                                       <  1:318265/100‑1 (MQ=255)
          cGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCt                                                                                      >  1:2913443/1‑100 (MQ=255)
          cGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCt                                                                                      <  1:1788671/100‑1 (MQ=255)
           gAGTTGTCTGCCTTCAACGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTa                                                                                     <  1:473729/100‑1 (MQ=255)
             gTTGTCTGCCTTCACCGACAGATAATCCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAAc                                                                                   <  2:2562759/100‑1 (MQ=255)
             gTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAAc                                                                                   >  2:1111713/1‑100 (MQ=255)
               tGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACgg                                                                                 <  2:344271/100‑1 (MQ=255)
                gTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGt                                                                                >  2:1765961/1‑100 (MQ=255)
                 tCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTg                                                                               >  1:1862493/1‑100 (MQ=255)
                 tCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTg                                                                               <  2:2780183/100‑1 (MQ=255)
                 tCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTg                                                                               <  1:1351150/100‑1 (MQ=255)
                 tCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTg                                                                               <  1:3720361/100‑1 (MQ=255)
                    gCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGATAACGGTgccg                                                                            >  1:2665265/1‑100 (MQ=255)
                     ccgtCACCGACAGATAAACCGGTGTTGCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTgccgc                                                                           <  1:396881/97‑1 (MQ=255)
                     ccTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTgccgc                                                                           >  2:2467655/1‑100 (MQ=255)
                      cTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTgccgcc                                                                          <  1:1371940/100‑1 (MQ=255)
                         cACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGccacc                                                                       <  1:544573/100‑1 (MQ=255)
                           ccGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCAt                                                                     >  1:1151838/1‑100 (MQ=255)
                           ccGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACATCTACCGTCGCTAACGGTGCCGCCACCAt                                                                     >  2:1196904/1‑100 (MQ=255)
                           ccGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTAACCTCTACCGTCGCTAACGGTGCCGCCACCAt                                                                     >  2:2444597/1‑100 (MQ=255)
                            cGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATa                                                                    >  1:3415245/1‑100 (MQ=255)
                              aCAGATAAACCGGTTTTTCCGGGCGCGGTTTCTGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATaaa                                                                  <  1:3691408/100‑1 (MQ=255)
                              aCAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATaaa                                                                  <  2:1940068/100‑1 (MQ=255)
                              aCAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATaaa                                                                  <  1:1771787/100‑1 (MQ=255)
                               cAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAg                                                                 <  1:952883/100‑1 (MQ=255)
                               cAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAg                                                                 <  1:84124/100‑1 (MQ=255)
                                   tAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATaa                                                             <  2:2178908/100‑1 (MQ=255)
                                    aaaCCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAAt                                                            >  2:3778983/1‑100 (MQ=255)
                                      aCCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATca                                                          >  1:1096486/1‑100 (MQ=255)
                                         ggTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATcagca                                                       <  1:2395767/100‑1 (MQ=255)
                                         ggTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATcagca                                                       >  1:1605493/1‑100 (MQ=255)
                                         ggTTTTTCCGGGCGCGCCTACGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACAATAAAGATAATcagca                                                       <  2:1836974/100‑1 (MQ=255)
                                         ggTTTTTCCGGCCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATcagca                                                       <  2:3194807/100‑1 (MQ=255)
                                          gTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATcagcag                                                      <  1:2298621/100‑1 (MQ=255)
                                          gTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATcagcag                                                      <  2:2041185/100‑1 (MQ=255)
                                           tttttCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGa                                                     >  1:3499341/1‑100 (MQ=255)
                                           tttttCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGa                                                     <  2:3473580/100‑1 (MQ=255)
                                              ttCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAAcc                                                  <  1:3081856/100‑1 (MQ=255)
                                               tCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCCAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGGTAATCAGCAGAACCa                                                 <  2:3043432/100‑1 (MQ=255)
                                               tCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCa                                                 <  2:2956941/100‑1 (MQ=255)
                                                ccGTGCCCGCCTGCGCCGTCCTGGCAGAAGCAGGCAAGTCCACGTTCACCTCTACCGTCGCTAACGCTGCCGCCACCATAAAGATAACCAGCAGAACCaa                                                >  2:583388/1‑100 (MQ=255)
                                                 cGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACTGTGCCGCCACCATAAAGATAATCAGCAGAACCAAc                                               >  2:3477660/1‑100 (MQ=255)
                                                   ggCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACAt                                             >  1:2678389/1‑100 (MQ=255)
                                                     cgcgGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCa                                           <  1:1649046/100‑1 (MQ=255)
                                                     cgcgGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCa                                           >  2:506521/1‑100 (MQ=255)
                                                     cgcgGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCa                                           >  2:1608143/1‑100 (MQ=255)
                                                     cgcgGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCa                                           >  2:1129384/1‑100 (MQ=255)
                                                        ggCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGt                                        <  1:632580/100‑1 (MQ=255)
                                                         gCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTc                                       <  2:885140/100‑1 (MQ=255)
                                                         gCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTc                                       >  2:3277449/1‑100 (MQ=255)
                                                            gccgcgTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGAt                                    <  2:852924/97‑1 (MQ=255)
                                                             cggcgTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATa                                   <  1:441699/100‑1 (MQ=255)
                                                              ggcgTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATaa                                  >  2:3671353/1‑100 (MQ=255)
                                                              ggcgTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATaa                                  <  2:147900/100‑1 (MQ=255)
                                                                cgTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAAc                                <  1:3000772/100‑1 (MQ=255)
                                                                 gTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACg                               <  1:342007/100‑1 (MQ=255)
                                                                 gTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACg                               >  1:1241506/1‑100 (MQ=255)
                                                                 gTGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACg                               <  1:516586/100‑1 (MQ=255)
                                                                  tGCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACgg                              <  1:169087/100‑1 (MQ=255)
                                                                   gCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGc                             >  2:1773674/1‑100 (MQ=255)
                                                                   gCTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGc                             >  2:1275585/1‑100 (MQ=255)
                                                                    cTGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCg                            <  2:1717342/100‑1 (MQ=255)
                                                                     tGGTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGt                           >  2:416578/1‑100 (MQ=255)
                                                                      ggTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTc                          >  1:2725787/1‑100 (MQ=255)
                                                                      ggTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTc                          >  1:1546224/1‑100 (MQ=255)
                                                                      ggTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTc                          >  2:1296853/1‑100 (MQ=255)
                                                                       gTAGAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCa                         >  2:345555/1‑100 (MQ=255)
                                                                          gAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGt                      <  2:152130/100‑1 (MQ=255)
                                                                          gAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGt                      >  1:1528182/1‑100 (MQ=255)
                                                                          gAAGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGt                      <  1:1969360/100‑1 (MQ=255)
                                                                           atctaGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGtt                     <  1:745045/96‑1 (MQ=255)
                                                                            aGCAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTg                    >  1:909908/1‑100 (MQ=255)
                                                                              cAGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGat                  <  1:63744/100‑1 (MQ=255)
                                                                               aGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGata                 <  2:758008/100‑1 (MQ=255)
                                                                               aGGCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGata                 <  2:2729169/100‑1 (MQ=255)
                                                                                ggCAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGatat                <  2:1074413/100‑1 (MQ=255)
                                                                                  cAAGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCa              <  1:1299421/100‑1 (MQ=255)
                                                                                    aGTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATACTCACCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATg            >  2:1710249/1‑100 (MQ=255)
                                                                                     gTTCACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGc           >  2:3599019/1‑100 (MQ=255)
                                                                                      ttcaccttcaccTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCa          >  1:1243667/1‑100 (MQ=255)
                                                                                        cACCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCAtt        >  2:1294322/1‑100 (MQ=255)
                                                                                         aCCTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCAttt       >  2:1011139/1‑100 (MQ=255)
                                                                                          ccTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTc      <  1:2168559/100‑1 (MQ=255)
                                                                                          ccTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTc      >  1:360325/1‑100 (MQ=255)
                                                                                          ccTTCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTc      <  1:3355922/100‑1 (MQ=255)
                                                                                           cTTCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCg     >  2:3491370/1‑100 (MQ=255)
                                                                                            ttCACCTCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGc    <  1:3805698/100‑1 (MQ=255)
                                                                                              caccTCTACCGTCGCTAACGGTGCCGCCAACATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCg  <  2:3695813/100‑1 (MQ=255)
                                                                                              cacaTCTACAGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCg  >  2:3113704/1‑100 (MQ=255)
                                                                                                 |                                                                                                
CGATAAACATCGAGTTGTCTGCCTTCACCGACAGATAAACCGGTTTTTCCGGGCGCGGCTGCGGCGTGCTGGTAGAAGCAGGCAAGTTCACCTTCACATCTACCGTCGCTAACGGTGCCGCCACCATAAAGATAATCAGCAGAACCAACATCACGTCGATAAACGGCGTCACGTTGATATCATGCATTTCGCCG  >  NC_000913/3151015‑3151208

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: ‑