Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,573,145 G→A R120C (CGC→TGC)  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,1450GA98.5% 234.9 / ‑4.2 68R120C (CGC→TGC) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (31/36);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

GACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGC  >  NC_000913/3573013‑3573281
                                                                                                                                    |                                                                                                                                        
gACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGa                                                                                                                                   <  1:1129018/140‑1 (MQ=255)
 aCGGCTACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAg                                                                                                                                  <  1:302354/140‑1 (MQ=255)
 aCGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAg                                                                                                                                  <  1:971953/140‑1 (MQ=255)
  cGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGTGCGCCTAAGGTTTCATACGGGCACAGGTGAGt                                                                                                                                 <  2:1092429/140‑1 (MQ=255)
  cGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGt                                                                                                                                 <  1:1146251/140‑1 (MQ=255)
  cGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGt                                                                                                                                 <  2:486219/140‑1 (MQ=255)
   ggCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTa                                                                                                                                >  2:313399/1‑140 (MQ=255)
    gCGACCGTCCCAGTAGTTGAGTTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAc                                                                                                                               <  1:321684/140‑1 (MQ=255)
    gCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAc                                                                                                                               <  2:228078/140‑1 (MQ=255)
     cGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAcc                                                                                                                              <  1:514607/140‑1 (MQ=255)
       aCCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCtt                                                                                                                            <  2:1174391/140‑1 (MQ=255)
       aCCGTCCCAGTAGTTGAATTCCCCAACCACCGAGACCCGACGGGCGGTTGGAGCCCAGCCAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCtt                                                                                                                            >  2:258129/1‑140 (MQ=255)
           tCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAg                                                                                                                                  <  2:290372/130‑1 (MQ=255)
           tCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAg                                                                                                                                  >  1:290372/1‑130 (MQ=255)
           tCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGa                                                                                                                        >  1:1029702/1‑140 (MQ=255)
           tCCCAGTAGTTGAACTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGa                                                                                                                        >  1:337352/1‑140 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAc                                                                                                                               >  2:379462/1‑122 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAc                                                                                                                               <  1:379462/122‑1 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAc                                                                                                                               >  2:550927/1‑122 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTAc                                                                                                                               <  1:550927/122‑1 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCtt                                                                                                                            <  2:58089/125‑1 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCtt                                                                                                                            >  1:58089/1‑125 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAgg                                                                                                             >  2:266165/1‑140 (MQ=255)
                      gAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACAGGCACAGGTGAGTACCTTCAGATAATAGCCAgg                                                                                                             >  2:718215/1‑140 (MQ=255)
                           gCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATcc                                                                                                        <  2:413227/140‑1 (MQ=255)
                            cccAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAgg                                                                                                             >  1:1155205/1‑134 (MQ=255)
                            cccAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAgg                                                                                                             <  2:1155205/134‑1 (MQ=255)
                            cccAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCa                                                                                                       >  2:357334/1‑140 (MQ=255)
                                     cGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGa                                                                                              >  1:286374/1‑140 (MQ=255)
                                      gagaCCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGAt                                                                                             >  1:81978/1‑140 (MQ=255)
                                         aCCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAg                                                                                          <  2:706025/140‑1 (MQ=255)
                                         aCCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAg                                                                                          <  2:337352/140‑1 (MQ=255)
                                           ccGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCg                                                                                        <  1:313399/140‑1 (MQ=255)
                                              aCGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGAc                                                                                     <  2:216280/140‑1 (MQ=255)
                                              aCGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGAc                                                                                     >  1:528545/1‑140 (MQ=255)
                                                 ggCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCaa                                                                                  <  1:357334/140‑1 (MQ=255)
                                                 ggCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCaa                                                                                  <  1:1163722/140‑1 (MQ=255)
                                                 ggCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCaa                                                                                  <  1:266165/140‑1 (MQ=255)
                                                  gCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCaaa                                                                                 <  2:1029702/140‑1 (MQ=255)
                                                      ttGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACgg                                                                             <  1:1201548/140‑1 (MQ=255)
                                                              cagacagaGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGatc                                                                     <  1:718215/140‑1 (MQ=255)
                                                                    gagaAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCatcaat                                                               >  2:892555/1‑140 (MQ=255)
                                                                         aCGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGt                                                          >  1:619686/1‑140 (MQ=255)
                                                                          cgcgTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGtt                                                         >  2:44058/1‑140 (MQ=255)
                                                                               aCCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCt                                                    >  1:1007074/1‑140 (MQ=255)
                                                                               aCCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCt                                                    >  1:491760/1‑140 (MQ=255)
                                                                               aCCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCt                                                    >  1:116266/1‑140 (MQ=255)
                                                                                  tgtgACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGAc                                                 <  1:546089/140‑1 (MQ=255)
                                                                                  tgtgACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGAc                                                 <  2:491760/140‑1 (MQ=255)
                                                                                      aCGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATg                                             <  2:528545/140‑1 (MQ=255)
                                                                                       cGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGc                                            >  2:921021/1‑140 (MQ=255)
                                                                                           atccatAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGa                                        <  1:609322/140‑1 (MQ=255)
                                                                                                  gTATCTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGc                                 <  2:116266/140‑1 (MQ=255)
                                                                                                      cTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCTaat                             >  2:187899/1‑139 (MQ=255)
                                                                                                      cTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAAc                             >  2:44117/1‑140 (MQ=255)
                                                                                                      cTGCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAAc                             >  2:1032305/1‑140 (MQ=255)
                                                                                                        gCATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTg                           >  2:81768/1‑140 (MQ=255)
                                                                                                         cATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGaca                                      >  1:1102628/1‑128 (MQ=255)
                                                                                                         cATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGaca                                      <  2:1102628/128‑1 (MQ=255)
                                                                                                         cATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGaca                                      >  1:749345/1‑128 (MQ=255)
                                                                                                         cATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGaca                                      <  2:749345/128‑1 (MQ=255)
                                                                                                         cATGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGa                          <  2:1007074/140‑1 (MQ=255)
                                                                                                           tGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATa                        >  1:893473/1‑140 (MQ=255)
                                                                                                           tGCGCGCCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATa                        >  1:1119719/1‑140 (MQ=255)
                                                                                                               cgcCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGa                          >  1:641002/1‑134 (MQ=255)
                                                                                                               cgcCTAAGGTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGa                          <  2:641002/134‑1 (MQ=255)
                                                                                                                      ggTTTCATACGGGCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAAt             >  2:380727/1‑140 (MQ=255)
                                                                                                                                gggCACAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTcg   <  2:893473/140‑1 (MQ=255)
                                                                                                                                 ggCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTcgc  >  1:617886/1‑140 (MQ=255)
                                                                                                                                    |                                                                                                                                        
GACGGCGACCGTCCCAGTAGTTGAATTGCCCAACCACCGAGACCCGACGGGCGTTTGGAGCCCAGACAGAGAAACGCGTACCTGTGACGCCATCCATAGTATCTGCATGCGCGCCTAAGGTTTCATACGGGCGCAGGTGAGTACCTTCAGATAATAGCCAGGCATCCATTTCCTGGATTAGCGGACCAAAACGGTAAGGATCATCAATCAGGTTTTGCTGACCATGCCAGACAACAGCCAACTGATAGCGGAAAAAATTCTTACGTCGC  >  NC_000913/3573013‑3573281

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