Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,573,238 G→A 5.9% H89Y (CAT→TAT)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,573,2380GA 5.9% 210.6 / 2.3 67H89Y (CAT→TAT) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base A (2/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

TACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCGGGT  >  NC_000913/3573091‑3573383
                                                                                                                                                   |                                                                                                                                                 
tACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGcc                                                                                                                                                 <  1:930105/150‑1 (MQ=255)
  ccTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAg                                                                                                                                               <  2:812439/150‑1 (MQ=255)
   cTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGa                                                                                                                                              >  2:399961/1‑150 (MQ=255)
     gtgACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGaca                                                                                                                                            >  1:327184/1‑150 (MQ=255)
     gtgACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGaca                                                                                                                                            <  2:1261225/150‑1 (MQ=255)
      tgACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGacaa                                                                                                                                           >  1:645472/1‑150 (MQ=255)
         cGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAg                                                                                                                                        >  2:1339178/1‑150 (MQ=255)
               ccatAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACt                                                                                                                                  <  1:138454/150‑1 (MQ=255)
               ccatAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACt                                                                                                                                  <  2:119778/150‑1 (MQ=255)
               ccatAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACt                                                                                                                                  <  1:583219/150‑1 (MQ=255)
                     tATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGc                                                                                                                            <  2:162015/150‑1 (MQ=255)
                     tATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACAAAAACTGTAAGGATCATCATTCAGTTTTTGCTGACCATGCCAGACAACAGCCAAGTGATAGc                                                                                                                            <  2:634207/150‑1 (MQ=255)
                      aTCTGCATGCGCGCCTAAGGTTTGATACGGGCGCAGGTGAGTACGTTGAGATAATAGCCAGGCATCGATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATGAATGAGGTTTTGCTGACCATGCCAGACAAGAGGGAAGTGATAGCg                                                                                                                           <  2:1053598/150‑1 (MQ=255)
                         tGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGaa                                                                                                                        >  2:696424/1‑150 (MQ=255)
                         tGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATACCAGACAACAGCCAACTGATAGCGGaa                                                                                                                        >  2:1307676/1‑150 (MQ=255)
                         tGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATACCAGACAACAGCCAACTGATAGCGGaa                                                                                                                        >  2:1312979/1‑150 (MQ=255)
                          gCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGaaa                                                                                                                       >  2:1121821/1‑150 (MQ=255)
                               cgcgcCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAAtt                                                                                                                  >  2:597388/1‑150 (MQ=255)
                                gcgcCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTc                                                                                                                 <  2:817533/150‑1 (MQ=255)
                                   ccTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTa                                                                                                              <  2:645472/150‑1 (MQ=255)
                                           ttCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCgg                                                                                                      >  2:22527/1‑150 (MQ=255)
                                             cATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGaa                                                                                                    <  2:409362/150‑1 (MQ=255)
                                                aCGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGa                                                                                                 >  2:1314560/1‑150 (MQ=255)
                                                   ggCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGAcgc                                                                                              >  2:971383/1‑150 (MQ=255)
                                                     cgcAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGAcgccg                                                                                            >  2:442511/1‑150 (MQ=255)
                                                       cAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGac                                                                                                                                             <  2:1346446/99‑1 (MQ=255)
                                                       cAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGac                                                                                                                                             >  1:1346446/1‑99 (MQ=255)
                                                           tGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAg                                                                                      <  2:416653/150‑1 (MQ=255)
                                                                 ccTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGTCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATccc                                                                                >  2:89010/1‑150 (MQ=255)
                                                                 ccTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATccc                                                                                >  1:116546/1‑150 (MQ=255)
                                                                 ccTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATccc                                                                                <  1:696424/150‑1 (MQ=255)
                                                                   ttCAGATTATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGCCAACAGCCAACTGATAGAGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCg                                                                              >  2:1258454/1‑150 (MQ=255)
                                                                   ttCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATc                                                                                  <  2:643474/146‑1 (MQ=255)
                                                                   ttCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATc                                                                                  >  1:643474/1‑146 (MQ=255)
                                                                     cAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTg                                                                            <  1:639904/150‑1 (MQ=255)
                                                                       gATAATAGCCACGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCACGTTTTGCTGACCATGTCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAg                                                                          >  2:1366437/1‑150 (MQ=255)
                                                                        ataataGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAG‑ATCCCCGTGAGTc                                                                        >  1:958422/1‑150 (MQ=255)
                                                                              gCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGac                                                                   >  1:565277/1‑150 (MQ=255)
                                                                                  ggCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTc                                                               >  2:196379/1‑150 (MQ=255)
                                                                                       ccATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGtt                                                          >  2:675766/1‑150 (MQ=255)
                                                                                       ccATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCAAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACT‑CAGttt                                                         >  2:616975/1‑150 (MQ=255)
                                                                                        cATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGttt                                                         >  2:901907/1‑150 (MQ=255)
                                                                                                      gCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGcc                                           <  1:399961/150‑1 (MQ=255)
                                                                                                           ccAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGtt                                      >  1:625591/1‑150 (MQ=255)
                                                                                                             aaaaCGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGtttt                                    >  1:653517/1‑150 (MQ=255)
                                                                                                                 cGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGt                                >  1:1222666/1‑150 (MQ=255)
                                                                                                                   gTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTc                              <  1:901907/150‑1 (MQ=255)
                                                                                                                     aaGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCaa                            >  1:798065/1‑150 (MQ=255)
                                                                                                                      aGGATCATCAATCAGGTTTTTCTGACCATGCCAGACATCAGGCAAGTGATAGCGGAATAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGATACACTCCAGTTTTGCGTGTTTGCGCCCGGTTTTCGGTTCAAt                           <  2:1146613/150‑1 (MQ=255)
                                                                                                                      aGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAAt                           <  2:1266599/150‑1 (MQ=255)
                                                                                                                       ggATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATc                          >  2:178452/1‑150 (MQ=255)
                                                                                                                         atcatcAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCAc                        >  1:307932/1‑150 (MQ=255)
                                                                                                                           catcaATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCAccc                      <  2:473202/150‑1 (MQ=255)
                                                                                                                             tcaatcaGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCcac                    <  2:307932/150‑1 (MQ=255)
                                                                                                                             tcaatcaGGTTTTGCTGACCATACCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCcac                    <  1:1307676/150‑1 (MQ=255)
                                                                                                                             tcaatcaGGTTTTGCTGACCATACCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCcac                    <  1:1312979/150‑1 (MQ=255)
                                                                                                                              caatcaGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGCTTTTCGGTTCAATCACCcaca                   >  1:416872/1‑150 (MQ=255)
                                                                                                                               aatcaGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCcacac                  >  1:1386954/1‑150 (MQ=255)
                                                                                                                                  caGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATc               <  1:1311413/150‑1 (MQ=255)
                                                                                                                                  caGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATc               <  1:746448/150‑1 (MQ=255)
                                                                                                                                   aGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCg              <  1:1121821/150‑1 (MQ=255)
                                                                                                                                      ttttGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTa           <  1:178452/150‑1 (MQ=255)
                                                                                                                                             gACCATGTCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCGCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCgg    <  1:1366437/150‑1 (MQ=255)
                                                                                                                                              aCCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCggg   <  1:1314560/150‑1 (MQ=255)
                                                                                                                                              aCCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCggg   >  2:1390527/1‑150 (MQ=255)
                                                                                                                                              aCCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCggg   <  1:22527/150‑1 (MQ=255)
                                                                                                                                               ccATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCGGGt  <  1:1390527/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGCGGAATGACGCCGCTAAAGAATCCCCGTGAGTCGAGACACTCCAGTTTTGCGAGTTTGCGCCCGGTTTTCGGTTCAATCACCCACACATCGGTAGCGTCGGGT  >  NC_000913/3573091‑3573383

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: ‑