Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,249,628 G→A 100% G85D (GGC→GAC)  yqjD → membrane‑anchored ribosome‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,249,6280GA85.1% 225.0 / 29.7 94G85D (GGC→GAC) yqjDmembrane‑anchored ribosome‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (6/8);  new base A (39/41);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGT  >  NC_000913/3249532‑3249723
                                                                                                |                                                                                               
tccTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACAt                                                                                              <  1:1873087/99‑1 (MQ=255)
gCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACAt                                                                                              >  2:2314004/1‑100 (MQ=255)
gCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCATGCCGCTGAGTATGTGCGCGCACATCCGTGCCCGGGCGTGGACa                                                                                               >  1:1004018/1‑99 (MQ=255)
 ccTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACAtt                                                                                             <  2:1786228/100‑1 (MQ=255)
 ccTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACAta                                                                                             >  2:1873096/1‑99 (MQ=255)
   tGGGTGAAACCGGTGATGCCATTCCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTgg                                                                                           <  2:2542137/100‑1 (MQ=255)
        gAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTg                                                                                      <  1:2741816/100‑1 (MQ=255)
         aaaCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGc                                                                                     >  2:2004156/1‑100 (MQ=255)
         aaaCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGc                                                                                     <  2:1188585/100‑1 (MQ=255)
           aCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCaa                                                                                   >  2:2090962/1‑100 (MQ=255)
           aCCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCaa                                                                                   <  2:293436/100‑1 (MQ=255)
            ccGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAAt                                                                                  >  2:1234463/1‑100 (MQ=255)
               gTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCgg                                                                               <  1:2004145/100‑1 (MQ=255)
               gTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCgg                                                                               <  2:1211190/100‑1 (MQ=255)
               gTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCgg                                                                               <  2:1805324/100‑1 (MQ=255)
               gTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCgg                                                                               <  1:2342938/100‑1 (MQ=255)
                tGATGCCATTGCCAAACAAACCCGTGTCGCGGAGGCGCGCGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGgt                                                                              <  1:2047663/100‑1 (MQ=255)
                 gATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGgtg                                                                             <  2:1832883/100‑1 (MQ=255)
                 gATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGgtg                                                                             <  2:1254433/100‑1 (MQ=255)
                   taccATTGCCAAAAAAACCCGTGTCGCGGCGCCGCGTGCCGACGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTa                                                                           >  2:1967644/3‑100 (MQ=255)
                   tGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTa                                                                           >  1:656441/1‑100 (MQ=255)
                     ccATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTa                                                                           >  1:566264/1‑98 (MQ=255)
                     ccATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGt                                                                         <  2:212367/100‑1 (MQ=255)
                       aTTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGc                                                                       <  1:2753036/100‑1 (MQ=255)
                       aTTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGc                                                                       <  2:13711/100‑1 (MQ=255)
                           ccAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCgg                                                                   <  2:143003/100‑1 (MQ=255)
                           ccAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCgg                                                                   >  2:1594458/1‑100 (MQ=255)
                              aacaaacCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGt                                                                >  1:2765835/1‑100 (MQ=255)
                               acaaacCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGtt                                                               >  2:1410571/1‑100 (MQ=255)
                                caaacCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTc                                                              >  1:932676/1‑100 (MQ=255)
                                caaacCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTc                                                              >  2:2622192/1‑100 (MQ=255)
                                  aacCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctg                                                            >  2:991874/1‑100 (MQ=255)
                                   acCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                           <  1:2197231/100‑1 (MQ=255)
                                   acCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                           <  2:810787/100‑1 (MQ=255)
                                   acCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                           <  1:846084/100‑1 (MQ=255)
                                   acCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                           <  2:2713091/100‑1 (MQ=255)
                                   acCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGGGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                           <  1:2733907/100‑1 (MQ=255)
                                    cccGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgc                                                           >  1:2787703/1‑99 (MQ=255)
                                     ccGTGTCGCGGCGGCGCGTGCCGATGGGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTctgctg                                                         <  1:743407/100‑1 (MQ=255)
                                        tgtCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTcg                                                      <  1:940339/100‑1 (MQ=255)
                                        tgtCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTcg                                                      <  2:1019901/100‑1 (MQ=255)
                                         gtCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTcgc                                                     >  2:2693037/1‑100 (MQ=255)
                                           cgcgGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgt                                                   <  2:2436238/100‑1 (MQ=255)
                                           cgcgGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgt                                                   >  2:978258/1‑100 (MQ=255)
                                           cgcgGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGcgt                                                   >  1:8969/1‑100 (MQ=255)
                                                cggcgCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGTGTTCGGCGGTTGcttcggaa                                              >  1:2347233/1‑97 (MQ=255)
                                                 gtttctttccGATGATTATGTGCTCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTaa                                             <  2:2728082/92‑1 (MQ=255)
                                                   cgcgTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAAtt                                           >  2:951234/1‑100 (MQ=255)
                                                    gcgTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTa                                          <  2:1228334/100‑1 (MQ=255)
                                                    gcgTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTa                                          <  1:2622176/100‑1 (MQ=255)
                                                      gTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATg                                        >  1:2617624/1‑100 (MQ=255)
                                                       tGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATgg                                       >  1:53201/1‑100 (MQ=255)
                                                        gCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGc                                      <  1:2274517/100‑1 (MQ=255)
                                                         ccGATGAGTATGTGCGTGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCg                                     <  1:1700862/100‑1 (MQ=255)
                                                         ccGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCg                                     <  2:1417391/100‑1 (MQ=255)
                                                         ccGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCg                                     >  1:2421628/1‑100 (MQ=255)
                                                         ccGATGAGTAGGTGCGCGACCATCCGTTCGCGGGCCTGTACCTTGGCGCTACCATCCGCGTACTGCTCGGCGTTCTGCTCTCGCGTCGCTAATTATGGCg                                     <  1:1585742/100‑1 (MQ=255)
                                                          cGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCgg                                    >  2:168127/1‑100 (MQ=255)
                                                          cGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATGATGGCgg                                    >  1:2854719/1‑100 (MQ=255)
                                                           gATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCgg                                    >  1:1359181/1‑99 (MQ=255)
                                                           gATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGa                                   >  2:2478087/1‑100 (MQ=255)
                                                           gATGAGTATGTACGGGAAAATCCGTGGACGGGCGTGGACATTGGCTCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGGCGTTAATTATGGCg                                     <  1:1391166/98‑1 (MQ=255)
                                                            aTGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGac                                  >  2:903618/1‑100 (MQ=255)
                                                              gAGTATGTGCGCGAAAATCCGTGGACCGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGacac                                <  2:2378985/100‑1 (MQ=255)
                                                                gTATGTGCGCGAAAATCCGTGGACGGGTGTGGGCATTCGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGacac                                <  1:2127608/98‑1 (MQ=255)
                                                                 tATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtca                             <  2:134599/100‑1 (MQ=255)
                                                                 tATGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtca                             <  1:1579847/100‑1 (MQ=255)
                                                                  aTGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtcat                            >  2:2729648/1‑100 (MQ=255)
                                                                  aTGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtcat                            >  1:1444157/1‑100 (MQ=255)
                                                                  aTGTGCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtcat                            >  1:1409030/1‑100 (MQ=255)
                                                                    gtgCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACtcatca                          >  2:692333/1‑100 (MQ=255)
                                                                     tgCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCAc                         >  2:922630/1‑100 (MQ=255)
                                                                     tgCGCGAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCAc                         >  1:940612/1‑100 (MQ=255)
                                                                       cgcgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGc                       >  2:1254698/1‑100 (MQ=255)
                                                                        gcgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGca                      >  2:95283/1‑100 (MQ=255)
                                                                         cgAAAATCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGcac                     >  1:1054049/1‑100 (MQ=255)
                                                                            aaaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAg                  <  2:915855/100‑1 (MQ=255)
                                                                            aaaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAg                  >  2:1793689/1‑100 (MQ=255)
                                                                             aaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACaa                   >  1:2902194/1‑98 (MQ=255)
                                                                             aaTCCGTGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAgg                 >  2:287161/1‑100 (MQ=255)
                                                                                   tGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCgg           <  1:903615/100‑1 (MQ=255)
                                                                                   tGGACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCgg           <  1:1156412/100‑1 (MQ=255)
                                                                                    ggACGGGCGTGGACATTGGCGCTGCAATCGTTGTAGTGCTCGCCGTTCTGCTGTCGCTCCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCgg           <  1:970824/99‑1 (MQ=255)
                                                                                    ggACGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGt          <  1:888371/100‑1 (MQ=255)
                                                                                    ggACGGGCGTGGACAATGGCGCTGCAAGCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGt          <  1:421117/100‑1 (MQ=255)
                                                                                       cGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTaaa       <  1:2692930/100‑1 (MQ=255)
                                                                                       cGGGCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTaaa       <  2:2040849/100‑1 (MQ=255)
                                                                                         ggCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg     <  2:318334/100‑1 (MQ=255)
                                                                                         ggCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg     <  2:318252/100‑1 (MQ=255)
                                                                                         ggCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg     >  2:2186325/1‑100 (MQ=255)
                                                                                         ggCGTGGACATTGGCGCTGCAAACGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAg     >  1:1688906/1‑100 (MQ=255)
                                                                                          gCGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGc    <  2:1155012/100‑1 (MQ=255)
                                                                                           cGTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTTGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCg   >  1:1972686/1‑100 (MQ=255)
                                                                                            gTGGACATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGt  >  1:1865418/1‑100 (MQ=255)
                                                                                                |                                                                                               
GCCTGGGTGAAACCGGTGATGCCATTGCCAAACAAACCCGTGTCGCGGCGGCGCGTGCCGATGAGTATGTGCGCGAAAATCCGTGGACGGGCGTGGGCATTGGCGCTGCAATCGGTGTAGTGCTCGGCGTTCTGCTGTCGCGTCGTTAATTATGGCGGACACTCATCACGCACAAGGGCCCGGTAAAAGCGT  >  NC_000913/3249532‑3249723

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