Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 689,978 A→G 6.5% V220A (GTG→GCG)  PP_0591 ← adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947689,9780AG 6.5% 178.3 / 2.6 61V220A (GTG→GCG) PP_0591adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (31/26);  new base G (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

TTTCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTGGGTGGGTGGC  >  NC_002947/689830‑690125
                                                                                                                                                    |                                                                                                                                                   
tttCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCa                                                                                                                                                     >  5:178705/1‑149 (MQ=255)
      accacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCAt                                                                                                                                               >  6:90654/1‑149 (MQ=255)
        cacGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATgg                                                                                                                                             >  5:215359/1‑149 (MQ=255)
         acGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTTGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGTGCGCGCCAtcga                                                                                                                                            <  1:258378/149‑4 (MQ=255)
                 gCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGt                                                                                                                                    >  3:218290/1‑149 (MQ=255)
                    tGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTg                                                                                                                                 >  3:55689/1‑149 (MQ=255)
                        gAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCt                                                                                                                             >  1:136703/1‑149 (MQ=255)
                        gAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCt                                                                                                                             >  8:293954/1‑149 (MQ=255)
                           gATGTTGTGCTGGCTCATATGGTCGAACACGCATAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTc                                                                                                                           <  3:199906/148‑1 (MQ=255)
                                         tCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGttt                                                                                                                                  >  2:238539/1‑127 (MQ=255)
                                         tCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGttt                                                                                                                                  <  1:238539/127‑1 (MQ=255)
                                          cATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGt                                                                                                                                    >  2:251076/1‑124 (MQ=255)
                                          cATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGt                                                                                                                                    <  1:251076/124‑1 (MQ=255)
                                                 tCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGAt                                                                                                    >  5:66741/1‑149 (MQ=255)
                                                 tCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGAt                                                                                                    <  4:55689/149‑1 (MQ=255)
                                                   gAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATg                                                                                                   <  4:91186/148‑1 (MQ=255)
                                                    aaCACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTa                                                                                                 >  3:25551/1‑149 (MQ=255)
                                                           agagCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGttt                                                                                                                                  >  1:181284/1‑109 (MQ=255)
                                                           agagCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGttt                                                                                                                                  <  2:181284/109‑1 (MQ=255)
                                                            gagCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGTGCGCGCCAtcga                                                                                                                                            >  7:63350/1‑95 (MQ=255)
                                                            gagCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGTGCGCGCCAtcga                                                                                                                                            <  8:63350/98‑4 (MQ=255)
                                                              gCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTcgggcggg                                                                                       <  7:293954/149‑1 (MQ=255)
                                                                 tGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTc                                                                                              >  8:51959/1‑139 (MQ=255)
                                                                 tGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGCCGAGGGCTTCCCAGATGTACTc                                                                                              <  7:51959/139‑1 (MQ=255)
                                                                      ttGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCcttctt                                                                                <  4:218290/148‑1 (MQ=255)
                                                                         gACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGcc                                                                            >  4:278349/1‑149 (MQ=255)
                                                                           cAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAg                                                                          <  6:178705/149‑1 (MQ=255)
                                                                             gCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCa                                                                        >  6:169345/1‑149 (MQ=255)
                                                                              cGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCAt                                                                       >  6:141096/1‑149 (MQ=255)
                                                                               gggCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATg                                                                      <  6:66741/149‑1 (MQ=255)
                                                                                ggCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATgg                                                                     <  6:215359/149‑1 (MQ=255)
                                                                                   acacCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGc                                                                   >  7:157191/1‑148 (MQ=255)
                                                                                      ccGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACg                                                               <  4:25551/149‑1 (MQ=255)
                                                                                          gAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGcga                                                           >  5:179296/1‑149 (MQ=255)
                                                                                                 aTCTGCTCGTCAATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGcc                                                    >  1:156239/1‑149 (MQ=255)
                                                                                                  tCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATggg                                                                    <  5:16401/132‑1 (MQ=255)
                                                                                                  tCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATggg                                                                    >  6:16401/1‑132 (MQ=255)
                                                                                                      cTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCg                                                <  5:90654/148‑1 (MQ=255)
                                                                                                       tcgtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGc                                               >  4:91342/1‑148 (MQ=255)
                                                                                                       tcgtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGc                                               >  1:79860/1‑148 (MQ=255)
                                                                                                         gtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTAc                                            >  7:42131/1‑149 (MQ=255)
                                                                                                         gtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTAc                                            <  1:51860/149‑1 (MQ=255)
                                                                                                         gtcgATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCCTCCCACATGTACACGGGCGGGCCTTCTTCGCCAGCAAGGGCAACGGCGACGAAGCCTTCGCTAc                                            >  7:234286/1‑149 (MQ=255)
                                                                                                             atgatgCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACATTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGgc                                        >  5:306481/1‑149 (MQ=255)
                                                                                                                       gCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGaa                               >  1:251188/1‑148 (MQ=255)
                                                                                                                        cATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACTTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAaca                             >  1:70975/1‑149 (MQ=255)
                                                                                                                               cgcTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTg                      <  8:42131/149‑1 (MQ=255)
                                                                                                                                   cGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAAc                  <  3:278349/149‑1 (MQ=255)
                                                                                                                                    gTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGaa                   <  2:79860/147‑1 (MQ=255)
                                                                                                                                    gTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGaa                   <  2:70975/147‑1 (MQ=255)
                                                                                                                                    gTCCTCGATGGCGCCCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACt                 >  6:214852/1‑149 (MQ=255)
                                                                                                                                      ccTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTg               >  3:125480/1‑149 (MQ=255)
                                                                                                                                       cTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTg               <  5:141096/148‑1 (MQ=255)
                                                                                                                                       cTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGc              >  6:222914/1‑149 (MQ=255)
                                                                                                                                        tCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCt             >  1:205338/1‑149 (MQ=255)
                                                                                                                                        accaTGGCGTGCGCGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCt             >  4:128297/4‑149 (MQ=255)
                                                                                                                                          gATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtg            <  2:156239/148‑1 (MQ=255)
                                                                                                                                          gATGGCGCGCACGCCATGGTCGATACATTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtgg           <  6:306481/149‑1 (MQ=255)
                                                                                                                                           aTGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCtggg          <  6:179296/149‑1 (MQ=255)
                                                                                                                                                gcgcACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTgggtgggt     <  4:125480/149‑1 (MQ=255)
                                                                                                                                                  gcACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTGGgtgggtgg   >  7:16651/1‑149 (MQ=255)
                                                                                                                                                   cACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTGGGTGGGTGGc  <  1:86934/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTTCACACCACGCTCGAGCATGTCGAGGATGTTGTGCTGGCTCATATGGTCGAACACGCAGAGCTTGGTGTTGGACAGCGGGCACACCGTGAGCGGGATCTGCTCGTCGATGATGCGCTGCATCAGGCGCTCGTCCTCGATGGCGCGCACGCCATGGTCGATACGTTTGATCTTCAGCAGGTCGAGGGCTTCCCAGATGTACTCGGGCGGGCCTTCTTCGCCAGCATGGGCAACGGCGACGAAGCCTTCGCTACGGGCGCGGTCGAACACCCGCTGGAACTTGCTGGGTGGGTGGC  >  NC_002947/689830‑690125

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: