Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 711,370 G→C 6.8% A390P (GCG→CCG)  ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947711,3700GC 6.8% 189.8 / 4.1 59A390P (GCG→CCG) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (27/28);  new base C (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

ACCTGCAAGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGCGTG  >  NC_002947/711223‑711511
                                                                                                                                                    |                                                                                                                                             
aCCTGCAAGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTgc                                                                                                                                              >  2:251845/1‑149 (MQ=255)
    gCAAGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGAATGTCTACTCGCTGCCATTCTTCCGCGGCCAGTTCAACTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGACGGTGCGCt                                                                                                                                           >  7:296212/1‑148 (MQ=255)
     cAAGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGa                                                                                                                                         >  2:320752/1‑149 (MQ=255)
      aaGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGTGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCTCTGAt                                                                                                                                        <  5:267198/149‑1 (MQ=255)
      aaGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGAt                                                                                                                                        <  8:86949/149‑1 (MQ=255)
        gcgcTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGc                                                                                                                                      <  6:345524/149‑1 (MQ=255)
             aCGGCATGGTCCACGACGACATCCTTACCCCGCTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACAcc                                                                                                                                 <  4:285536/149‑1 (MQ=255)
              cGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCacac                                                                                                                                  <  1:320752/147‑1 (MQ=255)
                gCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCTATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGa                                                                                                                               >  8:256657/1‑148 (MQ=255)
                gCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGaa                                                                                                                              <  6:167750/149‑1 (MQ=255)
                gCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGa                                                                                                                               >  4:348221/1‑148 (MQ=255)
                     gTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGGCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTg                                                                                                                                               <  4:160690/127‑1 (MQ=255)
                     gTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTg                                                                                                                                               >  3:160690/1‑127 (MQ=255)
                      tCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATc                                                                                                                        <  8:53339/149‑1 (MQ=255)
                       cAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCAt                                                                                                                         >  3:256172/1‑147 (MQ=255)
                          cgacgaCATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGcc                                                                                                                    <  4:256172/149‑1 (MQ=255)
                             cgacATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCca                                                                                                                   >  1:21190/1‑147 (MQ=255)
                                 aTCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCTTAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCt                                                                                                              <  1:314686/149‑1 (MQ=255)
                                  tCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTAc                                                                                                            <  4:192778/149‑1 (MQ=255)
                                   ccTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCCACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGCTGCGCTGATGCACCCCGAAACCATCAGCCACAGCTTCa                                                                                                           >  4:101048/1‑149 (MQ=255)
                                       aCCCCGGTGCAGAGCCACGGCGTGTATGTCGACTCGCTGCCAGTCTTCGGCGGCCAGTTCATCTGGACGGCGCACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGGACAGGGCAGCCATCAGCGAGAGCTCCAAGCa                                                                                                       >  6:248018/1‑149 (MQ=255)
                                        ccccGGTGCAGAGTAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCAc                                                                                                      <  1:92190/149‑1 (MQ=255)
                                                  gagCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGgcgc                                                                                             >  3:134547/1‑148 (MQ=255)
                                                  gagCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGcc                                                                                            <  4:327527/149‑1 (MQ=255)
                                                   agCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCca                                                                                           >  6:76288/1‑149 (MQ=255)
                                                    gCAACGGCGTGTATGTCGACTCGCTTCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCTGCCTTCGTCGAGAAGCTGAGCGAAGTCGGTGCTCTGATGCACACCGAAACCATCTGCCACAGCTACATGCACTGCTGGCGCcac                                                                                          <  8:296212/149‑1 (MQ=255)
                                                            gtgtATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGAcccc                                                                                  >  4:331354/1‑149 (MQ=255)
                                                               tATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTCCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCt                                                                               >  6:30373/1‑149 (MQ=255)
                                                               tATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTCCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCt                                                                               >  6:24946/1‑149 (MQ=255)
                                                                   tCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATc                                                                           <  1:251845/149‑1 (MQ=255)
                                                                          gCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACcgcg                                                                    >  1:183749/1‑149 (MQ=255)
                                                                          gCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACcgcg                                                                    >  7:64071/1‑149 (MQ=255)
                                                                                      cggcggCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCgcgc                                                           <  3:331354/146‑1 (MQ=255)
                                                                                       ggcggcCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTg                                                       <  6:312260/149‑1 (MQ=255)
                                                                                       ggcggcCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTg                                                       <  5:299672/149‑1 (MQ=255)
                                                                                        gcggcCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACaa                                                                                        <  7:354793/115‑1 (MQ=255)
                                                                                        gcggcCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACaa                                                                                        >  8:354793/1‑115 (MQ=255)
                                                                                           gcCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTc                                                   >  1:235741/1‑149 (MQ=255)
                                                                                              aGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTg                                                >  2:247409/1‑149 (MQ=255)
                                                                                                   aTCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCa                                            <  2:235741/148‑1 (MQ=255)
                                                                                                   aTCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCAt                                           <  3:348221/149‑1 (MQ=255)
                                                                                                    tCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATg                                          >  4:288277/1‑149 (MQ=255)
                                                                                                     cTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATg                                          >  8:212280/1‑148 (MQ=255)
                                                                                                          aGGCAAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAag                                    <  6:165329/149‑1 (MQ=255)
                                                                                                           ggCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGCCCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGCCAag                                    >  2:324270/1‑148 (MQ=255)
                                                                                                                  cccGGCCAGCGTCGAGAAGCTGAGCGAAGGCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGa                             <  4:134547/149‑1 (MQ=255)
                                                                                                                  cccGCCCATCGGCGAGACGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTTCTGGCGCCCCAATCCCCAGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTTGGCATTGACAAGCAGCCGa                             <  8:171643/149‑1 (MQ=255)
                                                                                                                    cGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGc                           <  4:381721/149‑1 (MQ=255)
                                                                                                                     ggCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACATGGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCa                          >  1:274833/1‑149 (MQ=255)
                                                                                                                             tcgAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAg                  >  2:58942/1‑149 (MQ=255)
                                                                                                                              cgAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGc                 >  2:163028/1‑149 (MQ=255)
                                                                                                                               gagaAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGcc                >  5:188658/1‑149 (MQ=255)
                                                                                                                               gagaAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGc                 <  5:248018/148‑1 (MQ=255)
                                                                                                                                 gaAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGc              >  2:326186/1‑149 (MQ=255)
                                                                                                                                    gCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGc           <  1:58942/149‑1 (MQ=255)
                                                                                                                                    gCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGc           <  5:76288/149‑1 (MQ=255)
                                                                                                                                          cGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGc     >  2:49804/1‑149 (MQ=255)
                                                                                                                                          cGAAGTCGGTCCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGc     <  5:30373/149‑1 (MQ=255)
                                                                                                                                          cGAAGTCGGTCCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGc     <  5:24946/149‑1 (MQ=255)
                                                                                                                                          cGAAGGCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGc     <  7:256657/149‑1 (MQ=255)
                                                                                                                                           gAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGCg    >  1:309035/1‑149 (MQ=255)
                                                                                                                                             aGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGCGTg  >  2:131849/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
ACCTGCAAGCGCTACGGCATGGTCAACGACGACATCCTTACCCCGGTGCAGAGCAACGGCGTGTATGTCGACTCGCTGCCATTCTTCGGCGGCCAGTTCATCTGGAAGGCCAACCCGGCCATCGTCGAGAAGCTGAGCGAAGTCGGTGCGCTGATGCACACCGAAACCATCAGCCACAGCTACATGCACTGCTGGCGCCACAAGACCCCGCTGATCTACCGCGCCACCGCGCAGTGGTTCGTGGGCATGGACAAGCAGCCGAGCACCGGCGAGCCGCTGCGTGAGCGTG  >  NC_002947/711223‑711511

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