Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,645,927 T→A 9.4% Q491L (CAG→CTG)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,645,9270TA 9.4% 394.3 / 24.1 128Q491L (CAG→CTG) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (56/60);  new base A (5/7);  total (61/67)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 138 total aligned reads displayed.

CCACTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGG  >  NC_002947/1645778‑1646076
                                                                                                                                                     |                                                                                                                                                     
ccACTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggct                                                                                                                                                       <  1:1209345/150‑1 (MQ=255)
ccACTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggct                                                                                                                                                       >  1:533025/1‑150 (MQ=255)
ccACTGCGACTTGTGGGGTGGGGATCTCGAGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGCTACGATGATCGCAGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGTTTCggct                                                                                                                                                       <  2:3676042/150‑1 (MQ=255)
  aCTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggctgg                                                                                                                                                     >  1:4883137/1‑150 (MQ=255)
   cTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCggctggc                                                                                                                                                    <  1:1329487/150‑1 (MQ=255)
        ccTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCg                                                                                                                                               >  2:1209875/1‑150 (MQ=255)
         cTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGa                                                                                                                                              >  2:1011200/1‑150 (MQ=255)
            gtggggtgggGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccc                                                                                                                                           <  1:3661175/150‑3 (MQ=255)
             tggggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccct                                                                                                                                          <  1:493590/150‑4 (MQ=255)
              ggggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTcc                                                                                                                                         <  1:2111406/150‑1 (MQ=255)
              ggggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccctg                                                                                                                                         <  2:2029938/150‑5 (MQ=255)
              ggggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccctg                                                                                                                                         <  1:5037813/150‑5 (MQ=255)
               gggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCa                                                                                                                                        >  2:2622694/1‑150 (MQ=255)
               gggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCAGCCAGCCGAccctgg                                                                                                                                        <  1:971388/150‑6 (MQ=255)
                ggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAt                                                                                                                                       <  2:1980889/150‑1 (MQ=255)
                ggtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAt                                                                                                                                       <  1:3367928/150‑1 (MQ=255)
                 gtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATg                                                                                                                                      <  2:1462357/150‑1 (MQ=255)
                 gtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATg                                                                                                                                      >  1:4715645/1‑150 (MQ=255)
                 gtggggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATg                                                                                                                                      >  1:4715641/1‑150 (MQ=255)
                    gggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTcc                                                                                                                                   <  1:4015709/150‑1 (MQ=255)
                    gggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTcc                                                                                                                                   >  1:3838011/1‑150 (MQ=255)
                    gggATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTcc                                                                                                                                   >  1:831561/1‑150 (MQ=255)
                      gATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCca                                                                                                                                 <  2:3598540/150‑1 (MQ=255)
                       aTCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCcac                                                                                                                                <  2:1806625/150‑1 (MQ=255)
                           cgcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGt                                                                                                                            <  1:2840612/150‑1 (MQ=255)
                            gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                           >  1:1501962/1‑150 (MQ=255)
                            gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                           >  2:3993369/1‑150 (MQ=255)
                            gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTc                                                                                                                           >  1:4927125/1‑150 (MQ=255)
                            gcgGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCACCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACGGGTc                                                                                                                           >  2:2358727/1‑150 (MQ=255)
                                    gCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTggtggg                                                                                                                   >  2:2827634/1‑150 (MQ=255)
                                         aTTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                              <  2:3215771/150‑1 (MQ=255)
                                         aTTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAAc                                                                                                              <  1:4361246/150‑1 (MQ=255)
                                              cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                         <  1:1228459/150‑1 (MQ=255)
                                              cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                         <  2:1982982/150‑1 (MQ=255)
                                              cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                         <  2:1095309/150‑1 (MQ=255)
                                              cGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATa                                                                                                         <  1:961598/150‑1 (MQ=255)
                                               gTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAg                                                                                                        >  2:785665/1‑150 (MQ=255)
                                                   cGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTcc                                                                                                    <  1:4707040/150‑1 (MQ=255)
                                                    gTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCa                                                                                                   <  2:1871328/150‑1 (MQ=255)
                                                       gcggcgGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCg                                                                                                >  1:2348908/1‑150 (MQ=255)
                                                             gTCGGGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCtt                                                                                          >  2:1736263/1‑150 (MQ=255)
                                                                    aCGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTg                                                                                   <  1:361729/150‑1 (MQ=255)
                                                                    aCGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTg                                                                                   >  2:3376154/1‑150 (MQ=255)
                                                                    aCGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTg                                                                                   <  2:4126458/150‑1 (MQ=255)
                                                                        aCAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCa                                                                               <  2:4749932/150‑1 (MQ=255)
                                                                        aCAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCa                                                                               >  1:2499816/1‑150 (MQ=255)
                                                                        aCAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCa                                                                               <  2:5012259/150‑1 (MQ=255)
                                                                           ggCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGt                                                                            >  1:2321415/1‑150 (MQ=255)
                                                                           ggCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGt                                                                            <  1:16476/150‑1 (MQ=255)
                                                                            gCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTc                                                                           >  1:5034555/1‑150 (MQ=255)
                                                                              tGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCaa                                                                         <  1:3602870/150‑1 (MQ=255)
                                                                               gCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCATGCCGGTGGCCAGGTCAAg                                                                        <  2:1282285/150‑1 (MQ=255)
                                                                                ccTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAgg                                                                       >  1:1128105/1‑150 (MQ=255)
                                                                                 cttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAggg                                                                      >  2:166602/1‑150 (MQ=255)
                                                                                  ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGc                                                                     >  2:4989604/1‑150 (MQ=255)
                                                                                  ttgcttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGc                                                                     >  2:2375862/1‑150 (MQ=255)
                                                                                      ttgGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGAcc                                                                 >  1:1972720/1‑150 (MQ=255)
                                                                                         gTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACg                                                              >  1:321689/1‑150 (MQ=255)
                                                                                          tCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGc                                                             <  1:1900197/150‑1 (MQ=255)
                                                                                            ggTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCg                                                           >  2:1958301/1‑150 (MQ=255)
                                                                                                cGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGgcgc                                                       <  1:3506211/150‑1 (MQ=255)
                                                                                                cGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGgcgc                                                       <  2:588232/150‑1 (MQ=255)
                                                                                                 gatgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCt                                                      >  2:1337031/1‑150 (MQ=255)
                                                                                                   tgatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGt                                                    >  1:4560035/1‑150 (MQ=255)
                                                                                                    gatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTa                                                   <  1:3942041/150‑1 (MQ=255)
                                                                                                    gatCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTa                                                   <  1:2433395/150‑1 (MQ=255)
                                                                                                     atCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGc                                                             >  1:3630444/1‑139 (MQ=255)
                                                                                                     atCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTAc                                                  <  2:4848215/150‑1 (MQ=255)
                                                                                                      tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtggt                                                                                                                      >  2:692216/1‑81 (MQ=255)
                                                                                                      tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtggt                                                                                                                      <  1:692216/81‑1 (MQ=255)
                                                                                                      tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACt                                                 >  1:3566298/1‑150 (MQ=255)
                                                                                                      tCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCAGGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACt                                                 >  2:1410180/1‑150 (MQ=255)
                                                                                                          cGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTtctc                                             <  2:281186/150‑1 (MQ=255)
                                                                                                          cGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTtctc                                             <  1:4404260/150‑1 (MQ=255)
                                                                                                           ggTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCg                                            <  2:4887128/150‑1 (MQ=255)
                                                                                                               tCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCGGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGgcgc                                        >  2:3257412/1‑150 (MQ=255)
                                                                                                                 acacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCgtggtg                                                                                                                     >  1:3349998/1‑71 (MQ=255)
                                                                                                                 acacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGccc                                      >  1:1491181/1‑150 (MQ=255)
                                                                                                                  cacCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCg                                     >  1:3791279/1‑150 (MQ=255)
                                                                                                                    ccATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCggg                                   <  1:56207/150‑1 (MQ=255)
                                                                                                                    ccATCGTCGGTTGTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCggg                                   <  2:4221266/150‑1 (MQ=255)
                                                                                                                         gtcgGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCg                              >  1:3431988/1‑150 (MQ=255)
                                                                                                                              ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCgg                         >  2:2179338/1‑150 (MQ=255)
                                                                                                                              ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCgg                         >  2:2179392/1‑150 (MQ=255)
                                                                                                                              ttttCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCgg                         >  2:1094078/1‑150 (MQ=255)
                                                                                                                                   aGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGccc                    <  2:357846/150‑1 (MQ=255)
                                                                                                                                    ggTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGcccc                   <  2:3838011/150‑1 (MQ=255)
                                                                                                                                      tgtACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCAc                 <  2:4715641/150‑1 (MQ=255)
                                                                                                                                      tgtACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCAc                 <  2:4715645/150‑1 (MQ=255)
                                                                                                                                       gtACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACt                <  1:3850882/150‑1 (MQ=255)
                                                                                                                                        tACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTg               <  1:520221/150‑1 (MQ=255)
                                                                                                                                        tACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTg               <  2:1297999/150‑1 (MQ=255)
                                                                                                                                        tACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGGTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTg               <  1:520410/150‑1 (MQ=255)
                                                                                                                                         atggacgTCGGCAGCCAGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTgg              >  2:1453193/7‑150 (MQ=255)
                                                                                                                                          ccAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGt             >  2:1240066/1‑150 (MQ=255)
                                                                                                                                          ccAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGt             >  1:2192937/1‑150 (MQ=255)
                                                                                                                                            gacgTCGGCAGCCAGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTcc           >  2:3201364/4‑150 (MQ=255)
                                                                                                                                              ggTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCgg         <  1:3096717/150‑1 (MQ=255)
                                                                                                                                                     tggctgCCGACGTCCATGTCCCACAGGTCGTGGTGGTTGAACTGATAGTTCCAGCGGGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGGTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTgg  >  1:30993/1‑150 (MQ=255)
                                                                                                                                                     tggctgCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTgg  <  1:4795800/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACTGCGCCTTGTGGGGTGGGGATCTCGCGGATCGGCACGATTGGCGTACCGTCGCGGCGGTCGAGCACGTACAGGCTGCCTTGCTTGGTCGGTACGATGATCGCCGGTTTCACACCATCGTCGGTTTTCAGGTGTACCAGGGTCGGCTGGCTGCCGACGTCCATGTCCCACAGGTCGTGGTGGGTGAACTGATAGTTCCAGCGTGCCTTGCCGGTGGCCAGGTCAAGGGCGACCACGCCGGCGCTGTACTTCTCGGCGCCCGGGGTGCGGTCGGCGCCCCACTGGTCCGGGGTCTGG  >  NC_002947/1645778‑1646076

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