Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,582,565 C→A 5.6% Q943H (CAG→CAT) ‡ treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,582,5650CA 5.6% 504.5 / 13.1 153Q943H (CAG→CAT) ‡treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base C (65/78);  new base A (5/4);  total (70/83)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 156 total aligned reads displayed.

CCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCAC  >  NC_002947/4582417‑4582713
                                                                                                                                                    |                                                                                                                                                    
ccTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCt                                                                                                                                                     >  2:1662360/1‑150 (MQ=255)
 cTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtg                                                                                                                                                    >  2:3456047/1‑150 (MQ=255)
 cTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCtg                                                                                                                                                    >  1:2561641/1‑150 (MQ=255)
         ggCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCa                                                                                                                                            <  2:1232678/150‑1 (MQ=255)
           cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                          <  2:1338219/150‑1 (MQ=255)
           cTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAgg                                                                                                                                          <  1:500941/150‑1 (MQ=255)
            tgtgTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAggg                                                                                                                                         >  1:480257/1‑150 (MQ=255)
            tgtgTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAggg                                                                                                                                         >  1:1422860/1‑150 (MQ=255)
            tgtgTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCAAGCAGGCTCAggg                                                                                                                                         <  2:2933917/150‑1 (MQ=255)
               gtTTTGCCCGGCTCTCTTCAAGGGGACGTGCGGGTTCTCATTCAAAATCTATTAGGTATGCATCACCCTGCACCAAGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGTCCACCCTGAGCCTGATTGCTCAGGgtgt                                                                                                                                      <  2:2434106/150‑1 (MQ=255)
                ttttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTc                                                                                                                                     <  1:3987172/150‑1 (MQ=255)
                  ttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGa                                                                                                                                   >  1:3156811/1‑150 (MQ=255)
                  ttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGa                                                                                                                                   <  1:1163256/150‑1 (MQ=255)
                  ttGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGACACCCTGAGCAAGCAGGCTCAGGgtggcct                                                                                                                                   <  2:5313207/150‑5 (MQ=255)
                    gCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGAtg                                                                                                                                 >  1:3194715/1‑150 (MQ=255)
                     cccGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGgtggcctgct                                                                                                                                >  1:2942132/1‑143 (MQ=255)
                         gcgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTg                                                                                                                            <  1:960603/150‑1 (MQ=255)
                          cgcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                           <  1:2163423/150‑1 (MQ=255)
                           gcTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTggg                                                                                                                          >  2:2703492/1‑150 (MQ=255)
                            ctctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                         <  2:4711816/150‑1 (MQ=255)
                            ctctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                         <  2:4711822/150‑1 (MQ=255)
                             tctTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTgg                                                                                                                           <  1:3696455/147‑1 (MQ=255)
                                tCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAgg                                                                                                                     >  2:5109892/1‑150 (MQ=255)
                                 cAAGGGGCCGGGCGGGTTCACCTTCAATATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGc                                                                                                                    <  2:794039/150‑1 (MQ=255)
                                  aaGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGcc                                                                                                                   <  2:1860690/150‑1 (MQ=255)
                                      ggCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCa                                                                                                               <  2:3847767/150‑1 (MQ=255)
                                       gCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCAc                                                                                                              >  2:3028731/1‑150 (MQ=255)
                                           ggcgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCCGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTg                                                                                                          >  2:3554945/1‑150 (MQ=255)
                                             cgggTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGtt                                                                                                        <  2:4644203/150‑1 (MQ=255)
                                                 ttCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtgctg                                                                                                    <  1:5109892/150‑1 (MQ=255)
                                                    aCCTTCAAAATCGATGAGGTAGGCATCAACCTGCAACACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAg                                                                                                 <  1:33842/150‑1 (MQ=255)
                                                      cTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgt                                                                                               >  1:1867051/1‑150 (MQ=255)
                                                      cTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGAATAACCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgt                                                                                               >  1:3439036/1‑150 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGt                                                                                                         >  2:2162310/1‑139 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGt                                                                                                         <  1:2162310/139‑1 (MQ=255)
                                                       ttCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGgtc                                                                                              >  1:796427/1‑150 (MQ=255)
                                                           aaaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGtcgtcg                                                                                          >  2:1573977/1‑150 (MQ=255)
                                                            caTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATCAGCAGGCCACCCTGAGCGTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGa                                                                                         <  2:5450504/149‑1 (MQ=255)
                                                            aaTCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGa                                                                                         <  2:1867051/150‑1 (MQ=255)
                                                                gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                     <  1:4228174/150‑1 (MQ=255)
                                                                gATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                     <  2:2254715/150‑1 (MQ=255)
                                                                 aTGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                    <  1:5098931/150‑1 (MQ=255)
                                                                    aGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGAc                                                                                 <  1:1122348/150‑1 (MQ=255)
                                                                    aGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGAc                                                                                 <  1:1722683/150‑1 (MQ=255)
                                                                       tAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTgg                                                                              <  2:1457181/150‑1 (MQ=255)
                                                                            aTCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATc                                                                         <  1:996192/150‑1 (MQ=255)
                                                                             tCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTg                                                                                                          >  2:2856965/1‑116 (MQ=255)
                                                                             tCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTg                                                                                                          <  1:2856965/116‑1 (MQ=255)
                                                                               aCCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCt                                                                      <  1:3584778/150‑1 (MQ=255)
                                                                                  cTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTc                                                                   <  1:2490569/150‑1 (MQ=255)
                                                                                  cTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCACAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtg                                                                                                       <  2:293669/114‑1 (MQ=255)
                                                                                  cTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCACAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTtg                                                                                                       >  1:293669/1‑114 (MQ=255)
                                                                                    gCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCa                                                                 <  2:2576112/150‑1 (MQ=255)
                                                                                     caccacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAg                                                                >  2:3032402/1‑150 (MQ=255)
                                                                                      accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTGAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                               <  2:3862291/150‑1 (MQ=255)
                                                                                      accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                         >  2:1044916/1‑92 (MQ=255)
                                                                                      accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGc                                                                                                                         <  1:1044916/92‑1 (MQ=255)
                                                                                      accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                               <  2:4604476/150‑1 (MQ=255)
                                                                                      accacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                               <  1:4027122/150‑1 (MQ=255)
                                                                                        cacGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtgt                                                             >  2:4137926/1‑150 (MQ=255)
                                                                                         acGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgtgtg                                                            >  2:3126819/1‑150 (MQ=255)
                                                                                           gAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGTAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTc                                                          >  1:4939944/1‑150 (MQ=255)
                                                                                                  tGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCGGTTGCTGCAGGTCGTCGACCAAGGACTGGCTCTCGCTGTCCAGGGTGTCGCGATGc                                                   >  2:1991708/1‑150 (MQ=255)
                                                                                                   ggCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                               <  2:515289/137‑1 (MQ=255)
                                                                                                   ggCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAgt                                                               >  1:515289/1‑137 (MQ=255)
                                                                                                    gCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTg                                                 >  1:3670808/1‑150 (MQ=255)
                                                                                                     cccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTgg                                                <  1:4137926/150‑1 (MQ=255)
                                                                                                      ccAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTggg                                               >  2:1391302/1‑150 (MQ=255)
                                                                                                          gTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCcag                                           >  1:4537637/1‑150 (MQ=255)
                                                                                                            cttggTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAg                                                                                    >  2:3715837/4‑107 (MQ=255)
                                                                                                              aGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGGGTGTCGCGATGCTGGGCCAGCAgg                                       >  1:943166/1‑150 (MQ=255)
                                                                                                                gTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTc                                     >  1:2950529/1‑150 (MQ=255)
                                                                                                                   gCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAg                                  <  2:4929920/150‑1 (MQ=255)
                                                                                                                    ccGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAgg                                 <  1:2021630/150‑1 (MQ=255)
                                                                                                                    ccGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAgg                                 <  2:3617111/150‑1 (MQ=255)
                                                                                                                     cGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAggg                                <  2:4340034/150‑1 (MQ=255)
                                                                                                                     cGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAggg                                >  1:1295241/1‑150 (MQ=255)
                                                                                                                     cGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAggg                                >  1:1295059/1‑150 (MQ=255)
                                                                                                                        gtACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGt                             >  2:1382313/1‑150 (MQ=255)
                                                                                                                          acacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                           <  2:3156819/150‑1 (MQ=255)
                                                                                                                          acacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTgg                           <  2:1850317/150‑1 (MQ=255)
                                                                                                                           cacGCCTAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGGGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTggg                          >  2:4412657/1‑150 (MQ=255)
                                                                                                                           cacGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTggg                          <  2:2535499/150‑1 (MQ=255)
                                                                                                                             cGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtc                        >  1:5427729/1‑150 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCaa                                                                                     <  2:2382051/88‑1 (MQ=255)
                                                                                                                              gCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtca                       <  2:2970749/150‑1 (MQ=255)
                                                                                                                                   aGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCa                                                                                                  <  1:1306826/70‑1 (MQ=255)
                                                                                                                                   aGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCa                                                                                                  >  2:1306826/1‑70 (MQ=255)
                                                                                                                                    gCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcag                 >  2:3215808/1‑150 (MQ=255)
                                                                                                                                     cAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcagt                >  1:3076016/1‑150 (MQ=255)
                                                                                                                                     atcgacACCCTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGtcagttcagt                >  2:529803/6‑150 (MQ=255)
                                                                                                                                        gacACCCTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCt             >  2:302214/3‑150 (MQ=255)
                                                                                                                                          cACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTTGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTTTAGCTGTACAGTGTGTCGCGATGCTTGGCCAGCAGGTCGTTGGAGTGTGTCAGTTCATTTCTGa           <  2:3439036/150‑1 (MQ=255)
                                                                                                                                           aCCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGAt          >  1:4822604/1‑150 (MQ=255)
                                                                                                                                            cccTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATc         >  2:5386231/1‑150 (MQ=255)
                                                                                                                                            cccTGAGCAAGCAGGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCGGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGGGGGTCAGTGCAGTGCTGATc         >  2:3970648/1‑150 (MQ=255)
                                                                                                                                             ccTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGa           <  1:2696829/147‑1 (MQ=255)
                                                                                                                                             ccTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGa           >  2:2696829/1‑147 (MQ=255)
                                                                                                                                                 aGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGc    <  1:3554945/150‑1 (MQ=255)
                                                                                                                                                   ccTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGcac  >  1:1194353/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTTGTAAGGGCTGTGTTTTGCCCGGCGCTCTTCAAGGGGCCGGGCGGGTTCACCTTCAAAATCGATGAGGTAGGCATCACCCTGCACCACGAGCACCTGGCCCAGGTGCAGGTCGCCGTGTACACGCATAAGCAGGCCACCCTGAGCCTGCTTGCTCAGGGTGTCGATGTGCTGGGCCAGGCCGTCACGCTGTTGCTGCAGGTCGTCGACCAAGGACTGGCTATCGCTGTCCAGTGTGTCGCGATGCTGGGCCAGCAGGTCGAGGGCGTGGGTCAGTTCAGTGCTGATCTGTGCAC  >  NC_002947/4582417‑4582713

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