Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,724 Δ1 bp 5.7% coding (596/1974 nt) glgX ← limit dextrin alpha‑1,6‑glucohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,7240A. 5.7% 347.3 / 12.7 106coding (596/1974 nt)glgXlimit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base A (51/48);  new base . (2/4);  total (53/53)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 112 total aligned reads displayed.

TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCT  >  NC_000913/3570575‑3570872
                                                                                                                                                      |                                                                                                                                                    
tGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgca                                                                                                                                                      >  1:1128947/1‑150 (MQ=255)
 gATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcag                                                                                                                                                     >  2:1391117/1‑150 (MQ=255)
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcagc                                                                                                                                                    <  1:925086/150‑1 (MQ=255)
  aTTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGgcagc                                                                                                                                                    >  2:527344/1‑150 (MQ=255)
      cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtt                                                                                                                                                <  1:1625422/150‑1 (MQ=255)
      cATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGtt                                                                                                                                                <  2:694961/150‑1 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTcc                                                                                                                                              <  1:1650787/150‑1 (MQ=255)
        tCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTcc                                                                                                                                              <  1:894815/150‑1 (MQ=255)
         cgcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgc                                                                                                                                             <  1:995785/150‑4 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgcc                                                                                                                                            <  1:1231222/150‑5 (MQ=255)
          gcgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAg                                                                                                                                            >  1:1538250/1‑150 (MQ=255)
           cgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAgc                                                                                                                                           >  2:1608834/1‑150 (MQ=255)
           cgAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAgc                                                                                                                                           >  1:217346/1‑150 (MQ=255)
            gAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCCTACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGCGCAGtttgccag                                                                                                                                          <  2:911039/150‑7 (MQ=255)
              aaCTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGcgcg                                                                                                                                        >  1:1279799/1‑150 (MQ=255)
               aCTCATCCAGCGCCGTTCCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCgg                                                                                                                                       >  1:83714/1‑150 (MQ=255)
                    tCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATg                                                                                                                                  <  1:766240/150‑1 (MQ=255)
                        gcgcCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCATACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCa                                                                                                                              >  2:1318410/1‑150 (MQ=255)
                             gTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                          >  2:1651080/1‑149 (MQ=255)
                             gTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                                                                                          <  1:1651080/149‑1 (MQ=255)
                             gTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGt                                                                                                                         >  2:865006/1‑150 (MQ=255)
                                 cTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                                                                                     <  2:1613156/150‑1 (MQ=255)
                                   ggCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaaa                                                                                                                   >  2:253833/1‑150 (MQ=255)
                                     cGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATa                                                                                                                 >  1:994987/1‑150 (MQ=255)
                                        gCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                              <  2:796173/150‑1 (MQ=255)
                                        gCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                                                                              >  2:890719/1‑150 (MQ=255)
                                         cAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTg                                                                                                             <  1:1474052/150‑1 (MQ=255)
                                            gCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatc                                                                                                          <  2:1279799/150‑1 (MQ=255)
                                             cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                                                                                         <  1:1378872/150‑1 (MQ=255)
                                             cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                                                                                         >  2:577268/1‑150 (MQ=255)
                                             cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                                                                                         <  1:493605/150‑1 (MQ=255)
                                              aTACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcat                                                                                                        <  2:994987/150‑1 (MQ=255)
                                               tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                                                                                       >  2:1626018/1‑150 (MQ=255)
                                               tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGAAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCTGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                                                                                       <  2:181659/150‑1 (MQ=255)
                                                   ccGGATTCAGCTCAAACATAGCGAACGGGTTGTAACCCCAGTAGTTACTTTTCCCCATGCGTTGCAGACGTGGTTCTCTGGCATACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTGAAATAGTTGATCATCACCg                                                                                                   <  2:1157590/150‑1 (MQ=255)
                                                   ccGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                                                                                   >  1:1148113/1‑150 (MQ=255)
                                                   ccGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                                                                                   <  2:1524772/150‑1 (MQ=255)
                                                          cAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGccc                                                                                            >  2:708580/1‑150 (MQ=255)
                                                              gcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGCCGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCATCGGATGCCCGAgg                                                                                        >  2:1022690/1‑150 (MQ=255)
                                                               cAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAggg                                                                                       <  2:702646/150‑1 (MQ=255)
                                                                   cATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCttt                                                                                   >  1:1189024/1‑150 (MQ=255)
                                                                   cATCGCCACCCGGTAGTCACCCCAGAAGTTACTTAGCCCCCTGCGATGCAGCCGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCgtt                                                                                   >  2:1516982/1‑150 (MQ=255)
                                                                    aTCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGTATGCCCGAGGGCTTta                                                                                  <  1:62110/150‑1 (MQ=255)
                                                                    aTCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTta                                                                                  <  1:1119940/150‑1 (MQ=255)
                                                                    aTCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTta                                                                                  <  1:458657/150‑1 (MQ=255)
                                                                     tCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtat                                                                                 <  1:865006/150‑1 (MQ=255)
                                                                     tCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtat                                                                                 >  2:1142504/1‑150 (MQ=255)
                                                                     tCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtat                                                                                 >  1:694342/1‑150 (MQ=255)
                                                                      cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAg                                                                                         >  2:87184/1‑141 (MQ=255)
                                                                      cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAg                                                                                         <  1:87184/141‑1 (MQ=255)
                                                                      cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtata                                                                                <  1:198365/150‑1 (MQ=255)
                                                                        ccACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAg                                                                              <  2:1111550/150‑1 (MQ=255)
                                                                        ccACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAg                                                                              >  1:706226/1‑150 (MQ=255)
                                                                             gggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGCCGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCa                                                                         >  2:10801/1‑150 (MQ=255)
                                                                             gggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCa                                                                         <  2:509190/150‑1 (MQ=255)
                                                                              ggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc                                                                        <  2:217346/150‑1 (MQ=255)
                                                                              ggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc                                                                        <  1:1391117/150‑1 (MQ=255)
                                                                              ggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc                                                                        >  1:926471/1‑150 (MQ=255)
                                                                               gTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACg                                                                       <  1:391927/150‑1 (MQ=255)
                                                                                   tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                      <  2:403129/99‑1 (MQ=255)
                                                                                   tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTc                                                                                                                      >  1:403129/1‑99 (MQ=255)
                                                                                     aCCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtctc                                                                 <  1:597928/150‑1 (MQ=255)
                                                                                      ccccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCg                                                                <  1:996238/150‑1 (MQ=255)
                                                                                          agtagtTACTTAGCCCCATGCCTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCg                                                            >  2:1645039/1‑150 (MQ=255)
                                                                                           gtagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCgg                                                           >  2:850802/1‑150 (MQ=255)
                                                                                                tACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGAtct                                                      >  2:1032213/1‑150 (MQ=255)
                                                                                                  cTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTcc                                                    >  1:1598688/1‑150 (MQ=255)
                                                                                                   ttAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCg                                                   <  1:1667391/150‑1 (MQ=255)
                                                                                                    tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCgg                                                  >  1:1064722/1‑150 (MQ=255)
                                                                                                     aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCggg                                                 <  2:905429/150‑1 (MQ=255)
                                                                                                       ccccATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTg                                               <  2:647457/150‑1 (MQ=255)
                                                                                                        cccATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGc                                              >  2:986609/1‑150 (MQ=255)
                                                                                                         ccATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCa                                             >  1:363751/1‑150 (MQ=255)
                                                                                                          cATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCaa                                            >  1:716267/1‑150 (MQ=255)
                                                                                                           aTGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAg                                           <  1:1641931/150‑1 (MQ=255)
                                                                                                             gCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTa                                         <  2:1538250/150‑1 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa                                                                     >  1:613133/1‑120 (MQ=255)
                                                                                                               gTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa                                                                     <  2:613133/120‑1 (MQ=255)
                                                                                                                 tGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCATCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAGCGGGATCTCCGTGTGCAAGTACGtt                                     <  2:1148113/150‑1 (MQ=255)
                                                                                                                 tGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGtt                                     <  1:687290/150‑1 (MQ=255)
                                                                                                                 tGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGtt                                     <  2:61356/150‑1 (MQ=255)
                                                                                                                   cAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGCTGCAAGTACGTTaa                                   >  2:383709/1‑150 (MQ=255)
                                                                                                                    aGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAAt                                  <  1:1645039/150‑1 (MQ=255)
                                                                                                                     gACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATc                                 >  1:714289/1‑150 (MQ=255)
                                                                                                                       cGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCt                               >  1:462008/1‑150 (MQ=255)
                                                                                                                         tGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCttt                             >  1:713500/1‑150 (MQ=255)
                                                                                                                           gTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGa                           >  1:780260/1‑150 (MQ=255)
                                                                                                                                 tGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                     <  1:530346/150‑1 (MQ=255)
                                                                                                                                 gctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATggg                     >  1:599725/6‑150 (MQ=255)
                                                                                                                                  ggCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATCGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                    >  1:253451/1‑150 (MQ=255)
                                                                                                                                  ctggaaCTGCTGCCAGTGGCGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGc                    >  1:1478259/5‑150 (MQ=255)
                                                                                                                                        cTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATa              <  1:1318410/150‑1 (MQ=255)
                                                                                                                                        cTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATa              >  2:1157386/1‑150 (MQ=255)
                                                                                                                                         tGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaa             >  1:1309459/1‑150 (MQ=255)
                                                                                                                                          gcgcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATaaa            >  1:136805/1‑150 (MQ=255)
                                                                                                                                           cgcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAat           >  1:503937/1‑150 (MQ=255)
                                                                                                                                             gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAatg          >  2:1074960/1‑150 (MQ=255)
                                                                                                                                                   ggCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCCTTGACATGGGCTTCATAAATGATGGtg    >  1:1591443/1‑150 (MQ=255)
                                                                                                                                                     cagcagTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtgct  >  1:311163/1‑150 (MQ=255)
                                                                                                                                                     cagcagTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGtgct  >  2:1392673/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGATTGCATCGCGAAACTCATCCAGCGCCGTTTCTGGCGAGCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCGCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACCGGGATCTCCGGGTGCAAGTACGTTAATCCTTTGACATGGGCTTCATAAATGATGGTGCT  >  NC_000913/3570575‑3570872

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: