Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,569,437 T→C 7.2% H628R (CAC→CGC) ‡ 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,569,4370TC 7.2% 479.7 / 20.2 166H628R (CAC→CGC) ‡glgXlimit dextrin alpha‑1,6‑glucohydrolase
Reads supporting (aligned to +/- strand):  ref base T (68/86);  new base C (4/8);  total (72/94)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 169 total aligned reads displayed.

CGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAGCCAA  >  NC_000913/3569288‑3569586
                                                                                                                                                     |                                                                                                                                                     
cgcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGt                                                                                                                                                       >  2:1595002/1‑150 (MQ=255)
 gcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTg                                                                                                                                                      <  1:414684/150‑1 (MQ=255)
 gcgcCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTg                                                                                                                                                      <  2:626258/150‑1 (MQ=255)
    ccAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCa                                                                                                                                                   >  1:516939/1‑150 (MQ=255)
    ccAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCa                                                                                                                                                   >  2:252193/1‑150 (MQ=255)
     cAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCAc                                                                                                                                                  >  2:1012079/1‑150 (MQ=255)
      aaCATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACt                                                                                                                                                 >  2:446007/1‑150 (MQ=255)
        cATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTcc                                                                                                                                               <  1:278099/150‑1 (MQ=255)
         aTTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTccc                                                                                                                                              <  2:2081002/150‑1 (MQ=255)
         aTTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTccc                                                                                                                                              <  2:2216037/150‑1 (MQ=255)
          ttAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTcccc                                                                                                                                             <  1:1922165/150‑1 (MQ=255)
          ttAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTcccc                                                                                                                                             <  2:2380277/150‑1 (MQ=255)
           tAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCca                                                                                                                                            <  1:2195344/150‑1 (MQ=255)
             aGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTcccccat                                                                                                                                          <  1:1350430/150‑4 (MQ=255)
              gTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCcagca                                                                                                                                         <  1:1212835/150‑1 (MQ=255)
                taTCGTTGTTCTCTAAACTTACCATGACTATGTCGTTTTTTATCATCTCTGGATCTGACACAATCCGTGTGCAGGTCCCTGCCAGTCAGCCGTAATCACTGGGTTATACTCTTCTTCGAATGGGGGAATGGCGTGCCACTCACCAGCTgg                                                                                                                                       <  2:80376/149‑1 (MQ=255)
                 aTCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGt                                                                                                                                      >  2:631030/1‑150 (MQ=255)
                 aTCGTCCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCTTGTGCAGGTCGCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATTGCGTGCCACTCCCCAGCAGGt                                                                                                                                      <  2:337617/150‑1 (MQ=255)
                  tCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTccccca                                                                                                                                           <  2:222251/144‑3 (MQ=255)
                  tCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAGTGGCACGCCATTccccca                                                                                                                                           >  1:222251/1‑142 (MQ=255)
                  tCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTa                                                                                                                                     <  1:700341/150‑1 (MQ=255)
                  tCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTa                                                                                                                                     <  2:1358433/150‑1 (MQ=255)
                    gTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaaa                                                                                                                                   <  2:1979827/150‑1 (MQ=255)
                     ttcttcTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTaaaa                                                                                                                                  <  1:917602/150‑1 (MQ=255)
                     ttcttcTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCAGTGTGCAGGTCCCTGCCAGAAAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACACCCCAGCAGGTaaaa                                                                                                                                  >  1:1235722/1‑150 (MQ=255)
                      tcttctCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAAc                                                                                                                                 <  1:1131211/150‑1 (MQ=255)
                      tattctCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTTCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAAc                                                                                                                                 <  2:425397/148‑1 (MQ=255)
                        ttctCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACaa                                                                                                                               >  2:2095331/1‑150 (MQ=255)
                          ctctAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtc                                                                                                                             >  2:1125476/1‑150 (MQ=255)
                          ctctAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAAtc                                                                                                                             <  2:301191/150‑1 (MQ=255)
                          ctctAAACTAACCATGACTAACTACTTATAAATCATCACTGGAACACACACAATCCCTGTGCAGGTCCCTCCCAGACAGCCGTAATCACAGGCTAATACTCTCCAGCGAATGGGGGAAAGCCGTGCCACACCCCAGCAGGTAAAACaaac                                                                                                                             >  1:2423852/1‑148 (MQ=255)
                              aaaCTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTg                                                                                                                         <  1:1799007/150‑1 (MQ=255)
                              aaaCTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTg                                                                                                                         <  1:1595002/150‑1 (MQ=255)
                              aaaCTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTg                                                                                                                         >  2:2126963/1‑150 (MQ=255)
                               aaCTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGt                                                                                                                        <  2:660646/150‑1 (MQ=255)
                                 cTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTa                                                                                                                      <  2:1871233/150‑1 (MQ=255)
                                   aaCCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTGCCCAGCAGGTAAAACAATCTCTGTTAcc                                                                                                                    <  2:1180329/150‑1 (MQ=255)
                                   aaCCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTAcc                                                                                                                    <  2:1836759/150‑1 (MQ=255)
                                              aCTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGc                                                                                                         >  2:1221965/1‑150 (MQ=255)
                                              aCTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGc                                                                                                         <  1:2133065/150‑1 (MQ=255)
                                               cTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCg                                                                                                        <  1:384874/150‑1 (MQ=255)
                                                  cTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTa                                                                                                     >  1:1147237/1‑150 (MQ=255)
                                                       ttATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGc                                                                                                <  1:582830/150‑1 (MQ=255)
                                                          tcatcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGAt                                                                                             <  2:2374530/150‑1 (MQ=255)
                                                           catcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATc                                                                                            <  2:1474981/150‑1 (MQ=255)
                                                           catcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATc                                                                                            >  1:996883/1‑150 (MQ=255)
                                                            atcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCa                                                                                           >  1:2348276/1‑150 (MQ=255)
                                                            atcTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCa                                                                                           <  2:1851492/150‑1 (MQ=255)
                                                             tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaa                                                                                          >  1:1374776/1‑150 (MQ=255)
                                                             tctcTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaa                                                                                          >  2:421698/1‑150 (MQ=255)
                                                                ctGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCaaaaa                                                                                       >  2:1375253/1‑150 (MQ=255)
                                                                 tGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAg                                                                                      <  1:2102322/150‑1 (MQ=255)
                                                                  ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGc                                                                                     <  1:2058004/150‑1 (MQ=255)
                                                                  ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGc                                                                                     >  2:1211994/1‑150 (MQ=255)
                                                                  ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGc                                                                                     >  1:356649/1‑150 (MQ=255)
                                                                  ggAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGACACTCCCCAGGAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGAGATCaaaattg                                                                                     <  2:2452042/150‑4 (MQ=255)
                                                                    aaCACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGa                                                                                   <  2:1898554/150‑1 (MQ=255)
                                                                      cacacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATc                                                                                 >  2:1546476/1‑150 (MQ=255)
                                                                        cacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATccc                                                                               >  1:1973312/1‑149 (MQ=255)
                                                                        cacacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCg                                                                               <  2:1472404/150‑1 (MQ=255)
                                                                         acacaATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCACCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCga                                                                              >  2:1049027/1‑150 (MQ=255)
                                                                            caATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGagag                                                                           <  2:2348276/150‑1 (MQ=255)
                                                                                     tgCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTg                                                                  >  1:170821/1‑150 (MQ=255)
                                                                                       cAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCa                                                                <  2:533315/150‑1 (MQ=255)
                                                                                            cccTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGt                                                           >  1:500861/1‑150 (MQ=255)
                                                                                               tGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTc                                                        >  2:1223961/1‑150 (MQ=255)
                                                                                                gCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCg                                                       <  1:1953637/150‑1 (MQ=255)
                                                                                                 ccAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCgg                                                      >  1:101513/1‑150 (MQ=255)
                                                                                                   aGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGcc                                                    >  2:1629537/1‑150 (MQ=255)
                                                                                                   aGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGcc                                                    >  2:1017926/1‑150 (MQ=255)
                                                                                                   aGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGcc                                                    <  1:1777620/150‑1 (MQ=255)
                                                                                                     aCAGCCGTAATCACTTGGTTATCCTCTCCAGCGAATGGGGGAATGGCTTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCg                                                  <  2:28038/150‑1 (MQ=255)
                                                                                                      cAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGt                                                 >  1:1712081/1‑150 (MQ=255)
                                                                                                         ccGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGtttt                                              <  2:996883/150‑1 (MQ=255)
                                                                                                            tAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGcc                                           <  1:1101764/150‑1 (MQ=255)
                                                                                                                cACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTc                                       >  2:2237920/1‑150 (MQ=255)
                                                                                                                cACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTc                                       >  1:1751389/1‑150 (MQ=255)
                                                                                                                          tCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCtt                             >  2:1598487/1‑150 (MQ=255)
                                                                                                                               tcCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAgg                        <  1:2095331/150‑1 (MQ=255)
                                                                                                                                   gCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGa                    <  1:421698/150‑1 (MQ=255)
                                                                                                                                   gCGAATGGGGGAATGCCGTCCAACTCCCCAGCATGTAAAACAATCTCTGTTACCTCACGCGTGGCGTTAATTGCGATCAAAAAGCTATCCGAGAGCAGAATTTGCTTCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGa                    <  2:1796147/150‑1 (MQ=255)
                                                                                                                                    cGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAg                   <  1:143156/150‑1 (MQ=255)
                                                                                                                                    cGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTAAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCCTTTTGCCACTCATCCGTGCTTAAAGGTTGAg                   >  2:182310/1‑150 (MQ=255)
                                                                                                                                     gAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGc                  <  2:1174945/150‑1 (MQ=255)
                                                                                                                                     gAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGc                  >  2:577191/1‑150 (MQ=255)
                                                                                                                                      tgctggggAGTGGCACGCCATTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCa                 >  1:2397020/5‑150 (MQ=255)
                                                                                                                                      aaTGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCa                 <  2:516939/150‑1 (MQ=255)
                                                                                                                                      aaTGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCa                 <  1:1125476/150‑1 (MQ=255)
                                                                                                                                      aaTGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGAGCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCa                 <  1:1753184/150‑1 (MQ=255)
                                                                                                                                          ggggAGTGGCACGCCATTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGACACTCATCCGTGCTTAAAGGTTGAGCatata             >  2:1068243/1‑149 (MQ=255)
                                                                                                                                          ggggAATGGCGTGCCACTCCCCAGCACGTAAAACAATCTCTGTTACCTCAAGAGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATATGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATc             >  2:2133368/1‑150 (MQ=255)
                                                                                                                                            ggAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGa           >  1:2196418/1‑150 (MQ=255)
                                                                                                                                            ggAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGa           >  1:361311/1‑150 (MQ=255)
                                                                                                                                              aaTGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGAtt         <  2:957852/150‑1 (MQ=255)
                                                                                                                                               aTGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGAttt        >  1:681196/1‑150 (MQ=255)
                                                                                                                                                tGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTa       >  2:1114014/1‑150 (MQ=255)
                                                                                                                                                  gCGTGCCACTCACCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTAGCATATCGATTTAGcc    <  1:1049027/150‑1 (MQ=255)
                                                                                                                                                   cGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTTCTTAAAGGTTGAGCATATCGATTTAGcc    <  2:850981/150‑1 (MQ=255)
                                                                                                                                                   cGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAGcc    <  2:1752699/150‑1 (MQ=255)
                                                                                                                                                     tGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAGCCaa  <  1:496777/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGACTAACTCCTTTTTTATCATCTCTGGAACACACACAATCCGTGTGCAGGTCCCTGCCAGACAGCCGTAATCACTGGGTTATCCTCTCCAGCGAATGGGGGAATGGCGTGCCACTCCCCAGCAGGTAAAACAATCTCTGTTACCTCAAGCGTGGCGTTAATTGCGATCAAAAAGCGATCCGAGAGCAGAATTTGCAGCTGTTTCGGCCCGTTTTGCCACTCATCCGTGCTTAAAGGTTGAGCATATCGATTTAGCCAA  >  NC_000913/3569288‑3569586

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