Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,570,714:1 +T 15.2% coding (606/1974 nt) glgX ← glycogen debranching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,570,7141.T15.2% 644.0 / 130.2 204coding (606/1974 nt)glgXglycogen debranching enzyme
Reads supporting (aligned to +/- strand):  ref base . (91/82);  new base T (18/13);  total (109/95)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 215 total aligned reads displayed.

GCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACC  >  NC_000913/3570615‑3570813
                                                                                                    |                                                                                                   
 naGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCACATGCGTTGCAGACGTGGTTCACTGGCAAACTgc‑g                                                                                                    >  2:1385607/2‑100 (MQ=255)
 cAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTgc‑g                                                                                                    <  1:947415/100‑1 (MQ=255)
  aGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTAACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGGAAGCTgc‑gc                                                                                                   <  1:2806638/100‑1 (MQ=255)
   ggCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑Gcc                                                                                                  <  2:203481/100‑1 (MQ=255)
     cATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑Gcc                                                                                                  >  1:4241219/1‑98 (MQ=255)
       tACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTg                                                                                              <  2:932113/100‑1 (MQ=255)
        aCGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTgg                                                                                             >  2:1595342/1‑100 (MQ=255)
        aCGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTgg                                                                                             <  1:1221398/100‑1 (MQ=255)
         cGCCGGATGCAGCGCAAACATCGCCACCTGGTTGTAACCCCATTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgc                                                                                            <  1:760565/100‑1 (MQ=255)
         cGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTGACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgc                                                                                            <  2:3446465/100‑1 (MQ=255)
           ccGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgcag                                                                                          <  1:4398133/100‑1 (MQ=255)
            cGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGgcagc                                                                                         <  2:3947653/100‑1 (MQ=255)
              gATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGcagcag                                                                                       <  1:258381/100‑1 (MQ=255)
               aTGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGt                                                                                      <  1:1035815/100‑1 (MQ=255)
               aTGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGt                                                                                      <  2:2903391/100‑1 (MQ=255)
                tGCAGCGCAAACATCCCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtt                                                                                     <  1:510157/100‑1 (MQ=255)
                  cAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttg                                                                                   <  1:3939384/100‑3 (MQ=255)
                  cAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttg                                                                                   <  2:1871197/100‑3 (MQ=255)
                  cAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTcc                                                                                   >  1:2826448/1‑100 (MQ=255)
                   aGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCa                                                                                  >  2:102483/1‑100 (MQ=255)
                    gcgcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTTACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgcc                                                                                 <  2:1559075/100‑5 (MQ=38)
                    gcgcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgcc                                                                                 <  2:3119347/100‑5 (MQ=38)
                    gcgcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgcc                                                                                 <  1:787382/100‑5 (MQ=38)
                    acgcAAACATCGCCACCGGGTTGTAACCCCAGAAGTTACTTAGACCCATGCGTTGCAGACGAAGTTCAATGGCAAACTAC‑ACCACTGGCAGCAGTTCCa                                                                                  >  1:1383185/2‑99 (MQ=255)
                      gcAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAgcg                                                                               >  2:3361574/1‑100 (MQ=255)
                         aacaacGCCACCAGGTTGTAAACCCAGTAGTTACTTAGCCCCATGCGTTGCAGACTTGGTTCACTGGCAAACTGCTGCCAGTGGC‑GCAGtttgcca                                                                                <  1:370350/96‑6 (MQ=25)
                         aaCATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCgg                                                                            <  1:4145089/100‑1 (MQ=255)
                              cGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAACCTGC‑GCCACCGGCCGCAGTCCCAGCGCGGt                                                                           >  1:3373388/1‑96 (MQ=255)
                               gCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGc                                                                      >  1:2422704/1‑100 (MQ=255)
                                 cACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGccc                                                                    <  2:1736133/100‑1 (MQ=255)
                                 cACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGccc                                                                    <  2:3376000/100‑1 (MQ=255)
                                 cACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGAAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGccc                                                                    >  1:3349226/1‑100 (MQ=255)
                                  aCCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCa                                                                   >  2:602961/1‑100 (MQ=255)
                                   ccGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCaa                                                                  >  1:1361683/1‑100 (MQ=255)
                                   ccGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACGGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCaa                                                                  >  1:2965352/1‑100 (MQ=255)
                                     gggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAAtt                                                                >  2:4087954/1‑100 (MQ=255)
                                      ggTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTg                                                               <  1:1455584/100‑1 (MQ=255)
                                         tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                            <  1:3629738/100‑1 (MQ=255)
                                         tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                            <  2:2331905/100‑1 (MQ=255)
                                         tGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGttt                                                            <  1:1739435/100‑1 (MQ=255)
                                           tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                          <  1:2742487/100‑1 (MQ=255)
                                           tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                          >  1:1325816/1‑100 (MQ=255)
                                           tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                          <  1:1548216/100‑1 (MQ=255)
                                           tAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCa                                                          >  2:404579/1‑100 (MQ=255)
                                              ccccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCaact                                                       >  2:180854/1‑98 (MQ=255)
                                                ccAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGTTAGCAGTTCAAACTCGGTAATGCCCAATTGTTTCAAATa                                                      >  1:2151051/1‑99 (MQ=255)
                                                 cAGTAGTTACTTAGCCCCATGCGTGGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGt                                                    <  1:1481326/100‑1 (MQ=255)
                                                 cAGTAGTTACTTAGCCCCATGCGATCCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGTCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGt                                                    <  1:1211299/100‑1 (MQ=255)
                                                  agtagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGtt                                                   <  1:2225319/100‑1 (MQ=255)
                                                   gtagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTg                                                  <  2:1219801/100‑1 (MQ=255)
                                                   gtagtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTg                                                  >  2:886613/1‑100 (MQ=255)
                                                     agtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGat                                                <  1:2090157/100‑1 (MQ=255)
                                                      gtTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatc                                               >  1:4165708/1‑100 (MQ=255)
                                                       ttACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatca                                              <  2:2518517/100‑1 (MQ=255)
                                                         acatAGCCCCATGCGTGGCAGACGTTGTTGAGTGGCAAACTGC‑GCCACTGGCAGCAGGTCCAGCGCGGTAATGCCCAATTATGTCAAATAGTTGat                                                <  1:3534850/93‑1 (MQ=255)
                                                         aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGTGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                            >  2:4200532/1‑100 (MQ=255)
                                                         aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGatcatc                                            >  2:38655/1‑100 (MQ=255)
                                                         aCTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTAC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGat                                                <  1:4377934/96‑1 (MQ=255)
                                                          cTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGAtcatca                                           >  1:293302/1‑100 (MQ=255)
                                                            tAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCAcc                                         >  2:3589980/1‑100 (MQ=255)
                                                             aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                        <  2:579561/100‑1 (MQ=255)
                                                             aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                        <  1:1149256/100‑1 (MQ=255)
                                                             aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                        >  2:832057/1‑100 (MQ=255)
                                                             aGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCg                                        <  1:404578/100‑1 (MQ=255)
                                                              gCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCgg                                       >  1:3649268/1‑100 (MQ=255)
                                                              gCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCgg                                       >  1:3411221/1‑100 (MQ=255)
                                                                nccATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGAt                                     >  1:4212501/2‑100 (MQ=255)
                                                                  cATGCGTTGCAGACGTGGTTCACTGGCAAGCTGC‑GGCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGc                                   >  2:3028307/1‑100 (MQ=255)
                                                                  cATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGc                                   <  2:1484707/100‑1 (MQ=255)
                                                                  cATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGc                                   >  2:3258495/1‑100 (MQ=255)
                                                                  cATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGc                                   >  2:2646709/1‑100 (MQ=255)
                                                                      cGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGa                               >  2:3322524/1‑100 (MQ=255)
                                                                          gCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGc                           >  2:3333186/1‑100 (MQ=255)
                                                                            aGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGACtt                         <  1:4066870/100‑1 (MQ=255)
                                                                              aCGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTta                       <  2:2793517/100‑1 (MQ=255)
                                                                                gTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGTGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtata                     >  2:2410778/1‑100 (MQ=255)
                                                                                gTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTtata                     <  2:376878/100‑1 (MQ=255)
                                                                                  ggTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAg                   >  1:2395027/1‑100 (MQ=255)
                                                                                  ggTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAg                   >  2:1308916/1‑100 (MQ=255)
                                                                                    ttCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTg                 >  2:3057660/1‑100 (MQ=255)
                                                                                    ttCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTg                 <  1:1199973/100‑1 (MQ=255)
                                                                                        cTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc             >  2:314941/1‑100 (MQ=255)
                                                                                        cTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCAc             >  1:2785572/1‑100 (MQ=255)
                                                                                         tGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACg            >  2:1335809/1‑100 (MQ=255)
                                                                                         tGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCAGTAATGCCCAAATGGTTCAAATAGTTGAACATCCCCCGATGCCCGACGGCTGTACAAGTGCCAc             >  1:1111696/1‑99 (MQ=255)
                                                                                          ggCAAACTGT‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACgg           >  1:1842357/1‑100 (MQ=255)
                                                                                          ctggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACgg           >  1:681909/5‑100 (MQ=38)
                                                                                           tggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa          >  1:1287597/4‑100 (MQ=255)
                                                                                           tggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCGTTATAAGTGCCACGGa          >  1:4035550/4‑100 (MQ=255)
                                                                                           gCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGa          <  2:4343999/100‑1 (MQ=255)
                                                                                            ggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAt         >  2:2839958/3‑100 (MQ=255)
                                                                                            ggaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAt         >  2:467022/3‑100 (MQ=255)
                                                                                            cAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAt         >  1:550100/1‑100 (MQ=255)
                                                                                            cAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAt         >  1:1708646/1‑100 (MQ=255)
                                                                                             gaaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtc        >  2:2431474/2‑100 (MQ=255)
                                                                                              aaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGAtct       >  1:3314384/1‑100 (MQ=255)
                                                                                              aaCTGCTGCCAGTGGC‑GCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGCACGCCCGACGGCCTTATAAGTccccc             >  1:3749552/1‑92 (MQ=255)
                                                                                                  gc‑gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGCTCTCGAc   >  1:3178133/1‑100 (MQ=255)
                                                                                                  gc‑gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc   >  2:40760/1‑100 (MQ=255)
                                                                                                  gc‑gcCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGAc   >  1:1930243/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
GCAGGCATACGCCGGATGCAGCGCAAACATCGCCACCGGGTTGTAACCCCAGTAGTTACTTAGCCCCATGCGTTGCAGACGTGGTTCACTGGCAAACTGC‑GCCACTGGCAGCAGTTCCAGCGCGGTAATGCCCAATTGTTTCAAATAGTTGATCATCACCGGATGCCCGAGGGCTTTATAAGTGCCACGGATCTCGACC  >  NC_000913/3570615‑3570813

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