Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,562,821 C→G 5.5% T270S (ACT→AGT)  glpD → aerobic glycerol 3‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,562,8210CG 5.5% 348.9 / 6.9 102T270S (ACT→AGT) glpDaerobic glycerol 3‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (48/48);  new base G (3/3);  total (51/51)
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.95e-01

Only 100 of 104 total aligned reads displayed.

GCCTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCGTCCGC  >  NC_000913/3562672‑3562967
                                                                                                                                                     |                                                                                                                                                  
gCCTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCAc                                                                                                                                                    >  2:350367/1‑150 (MQ=255)
gCCTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGGGCATACCCAGAAGCAACCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCAc                                                                                                                                                    <  2:1074838/150‑1 (MQ=255)
 ccTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACt                                                                                                                                                   <  1:1086707/150‑1 (MQ=255)
  cTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTa                                                                                                                                                  <  1:646803/150‑1 (MQ=255)
  cTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTa                                                                                                                                                  >  1:1174319/1‑150 (MQ=255)
   ttCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTAc                                                                                                                                                 <  1:1146190/150‑1 (MQ=255)
       ccTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGAt                                                                                                                                             >  1:87656/1‑150 (MQ=255)
        cTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATg                                                                                                                                            >  1:769295/1‑150 (MQ=255)
             ggCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAg                                                                                                                                       >  2:1206868/1‑150 (MQ=255)
             ggCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGGGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGTAGTGCCGATgatgga                                                                                                                                       <  1:177596/150‑6 (MQ=255)
              gCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGt                                                                                                                                      <  2:1065562/150‑1 (MQ=255)
               cATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGGGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTGCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGTAGTGCCGATgatggaaa                                                                                                                                     <  1:132393/150‑8 (MQ=255)
                    gCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAg                                                                                                                                >  1:841172/1‑150 (MQ=255)
                    gCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCAGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAg                                                                                                                                >  2:1190754/1‑150 (MQ=255)
                     ccTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAgg                                                                                                                               <  2:1053856/150‑1 (MQ=255)
                       tGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCg                                                                                                                             <  2:484647/150‑1 (MQ=255)
                       tGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCg                                                                                                                             <  1:529563/150‑1 (MQ=255)
                         aTCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGAt                                                                                                                           >  2:385302/1‑150 (MQ=255)
                           cAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATcc                                                                                                                         <  2:1320095/150‑1 (MQ=255)
                            aaaGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCg                                                                                                                        >  1:1527354/1‑150 (MQ=255)
                             aaGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGa                                                                                                                       >  1:523421/1‑150 (MQ=255)
                             aaGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGa                                                                                                                       >  2:898780/1‑150 (MQ=255)
                                gCAGCCATATTGTGGGGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAACATAAACGTATTGTGTTCGTGATCCCGTGCAGGGACGAGATTTCCCTCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAg                                                                                                                    <  2:270486/150‑1 (MQ=255)
                                  aGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCg                                                                                                                  <  1:728754/150‑1 (MQ=255)
                                   gCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCgg                                                                                                                 >  1:960922/1‑150 (MQ=255)
                                   gCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCgg                                                                                                                 <  1:244195/150‑1 (MQ=255)
                                   gCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCgg                                                                                                                 >  1:1137634/1‑150 (MQ=255)
                                      atatTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGa                                                                                                              >  1:1283615/1‑150 (MQ=255)
                                       tatTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGaa                                                                                                             <  1:717895/150‑1 (MQ=255)
                                        atTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAg                                                                                                            >  2:1009551/1‑150 (MQ=255)
                                        atTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAg                                                                                                            >  1:798167/1‑150 (MQ=255)
                                                  cgcgcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAagaga                                                                                                  >  1:1394667/1‑150 (MQ=255)
                                                  cgcgcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAagaga                                                                                                  >  1:1192883/1‑150 (MQ=255)
                                                  cgcgcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAagaga                                                                                                  >  1:839150/1‑150 (MQ=255)
                                                     gcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTg                                                                                               <  2:1170101/150‑1 (MQ=255)
                                                     gcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTg                                                                                               <  2:688296/150‑1 (MQ=255)
                                                     gcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTg                                                                                               <  2:1006768/150‑1 (MQ=255)
                                                     gcgTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTg                                                                                               <  1:1462401/150‑1 (MQ=255)
                                                             aCCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAAt                                                                                       >  1:1196989/1‑150 (MQ=255)
                                                              cccAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAAtt                                                                                      <  2:188706/150‑1 (MQ=255)
                                                              cccAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAAtt                                                                                      >  1:1074510/1‑150 (MQ=255)
                                                                      caagcCTACATTCTGCAAAACGAATATAATAGTATTGTGTTCGTTATACCTTGGATTGAAGAGTTTTCCATCATCTTCAGTTGGGATGTCGATTACAAAGGATATCCGAAATGGGTGAAGATTGAATAGAGTGAAATCAATTACctgctt                                                                              <  2:839150/150‑2 (MQ=255)
                                                                      caagcCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACctgctg                                                                              <  1:385302/150‑1 (MQ=255)
                                                                      caagcCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACctgctg                                                                              <  1:716908/150‑1 (MQ=255)
                                                                      caagcCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCAACATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACctgcgg                                                                              >  1:58715/1‑148 (MQ=255)
                                                                       aagcCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGa                                                                             >  1:623091/1‑150 (MQ=255)
                                                                        agcCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGGGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGaa                                                                            <  1:1341241/150‑1 (MQ=255)
                                                                          ccTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAAtg                                                                          <  2:1174319/150‑1 (MQ=255)
                                                                          ccTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAAtg                                                                          <  2:87656/150‑1 (MQ=255)
                                                                                ttCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGtata                                                                    <  2:235564/150‑1 (MQ=255)
                                                                                   tGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAaca                                                                 >  2:1512465/1‑150 (MQ=255)
                                                                                   tGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAaca                                                                 >  2:839318/1‑150 (MQ=255)
                                                                                   tGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAaca                                                                 <  1:427457/150‑1 (MQ=255)
                                                                                    gCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAacac                                                                >  1:217386/1‑150 (MQ=255)
                                                                                            aaGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGaa                                                                            <  2:528903/130‑1 (MQ=255)
                                                                                            aaGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGaa                                                                            >  1:528903/1‑130 (MQ=255)
                                                                                            aaGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTa                                                        <  2:1283615/150‑1 (MQ=255)
                                                                                              gATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTaaa                                                      <  1:426946/150‑1 (MQ=255)
                                                                                                tAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTaaaaa                                                    >  2:350971/1‑150 (MQ=255)
                                                                                                tAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTaaaaa                                                    >  2:1071035/1‑150 (MQ=255)
                                                                                                  aaCGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGATAGTGAAATCAATTACCTGCTGAAtgtg                                                                        <  2:454571/128‑1 (MQ=255)
                                                                                                  aaCGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAAtgtg                                                                        >  1:454571/1‑128 (MQ=255)
                                                                                                  aaCGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGc                                                  <  2:1192883/150‑1 (MQ=255)
                                                                                                    ggTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCTCTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAg                                                <  2:844649/149‑1 (MQ=255)
                                                                                                      tATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGtt                                              >  1:973362/1‑150 (MQ=255)
                                                                                                       aTTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTa                                             <  1:350367/150‑1 (MQ=255)
                                                                                                       aTTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTa                                             <  1:1206868/150‑1 (MQ=255)
                                                                                                       aTTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTa                                             <  2:523421/150‑1 (MQ=255)
                                                                                                           tgtTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGcc                                         <  2:584249/150‑1 (MQ=255)
                                                                                                            gtTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCg                                        <  2:229988/150‑1 (MQ=255)
                                                                                                             ttCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGt                                       >  1:1536990/1‑150 (MQ=255)
                                                                                                                gTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTaa                                                       <  1:1440779/131‑1 (MQ=255)
                                                                                                                gTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTaa                                                       >  2:1440779/1‑131 (MQ=255)
                                                                                                                            atggaCGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTgg                        <  2:769295/150‑1 (MQ=255)
                                                                                                                             tggaCGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGa                       <  1:1037674/150‑1 (MQ=255)
                                                                                                                              ggaCGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGAc                      <  1:504766/150‑1 (MQ=255)
                                                                                                                               gaCTAGTTTTCCATCATCTGCACTACCGATGTCGAGTACAAAGTCGATCCGAAAGCGGTGAATATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTTACGATATCGTCTGGAcc                     <  2:1527354/150‑1 (MQ=255)
                                                                                                                               gaCGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAACAAGCAGTTAAGCCGTGCCGATATCGTCTGGAcc                     >  2:497818/1‑150 (MQ=255)
                                                                                                                               gaCGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGAcc                     >  1:1511310/1‑150 (MQ=255)
                                                                                                                               gaCGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGAcc                     <  1:884890/150‑1 (MQ=255)
                                                                                                                                 cGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTa                   >  1:527491/1‑150 (MQ=255)
                                                                                                                                   aGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACt                 >  1:873953/1‑150 (MQ=255)
                                                                                                                                    gTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTc                >  1:1261151/1‑150 (MQ=255)
                                                                                                                                    gTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTc                >  1:864423/1‑150 (MQ=255)
                                                                                                                                    gTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTc                <  1:39477/150‑1 (MQ=255)
                                                                                                                                    gTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTc                <  1:283100/150‑1 (MQ=255)
                                                                                                                                      tttCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCg              >  1:1034077/1‑150 (MQ=255)
                                                                                                                                      tttCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCg              >  1:815269/1‑150 (MQ=255)
                                                                                                                                      tttCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCg              >  1:216170/1‑150 (MQ=255)
                                                                                                                                       ttCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCgg             >  2:61550/1‑150 (MQ=255)
                                                                                                                                        ctcgacatcGGTAGTGCCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGgt            >  1:1384436/6‑150 (MQ=255)
                                                                                                                                          catcatCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGgtgt          >  1:1151161/1‑150 (MQ=255)
                                                                                                                                           gacatcGGTAGTGCCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGtgtg         >  1:1211596/3‑150 (MQ=255)
                                                                                                                                           ataatcGCCACTACCGATGTCGAGTACAAACGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGtgtg         >  2:441741/4‑150 (MQ=255)
                                                                                                                                            acatcGGTAGTGCCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGc        >  1:598470/2‑150 (MQ=255)
                                                                                                                                             catcGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCg       <  1:350971/150‑1 (MQ=255)
                                                                                                                                                cGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCGTcc    <  1:1512465/150‑1 (MQ=255)
                                                                                                                                                 ggCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCGTCCg   >  1:346182/1‑150 (MQ=255)
                                                                                                                                                  gCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCGTCCGc  >  1:1403650/1‑150 (MQ=255)
                                                                                                                                                  gCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCGTCCGc  <  2:1137634/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCTTCGCCTTATGGCATTCGCCTGATCAAAGGCAGCCATATTGTGGTGCCGCGCGTGCATACCCAGAAGCAAGCCTACATTCTGCAAAACGAAGATAAACGTATTGTGTTCGTGATCCCGTGGATGGACGAGTTTTCCATCATCGGCACTACCGATGTCGAGTACAAAGGCGATCCGAAAGCGGTGAAGATTGAAGAGAGTGAAATCAATTACCTGCTGAATGTGTATAACACGCACTTTAAAAAGCAGTTAAGCCGTGACGATATCGTCTGGACCTACTCCGGTGTGCGTCCGC  >  NC_000913/3562672‑3562967

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