Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,903,890 C→T 6.8% P285P (CCG→CCA)  spuH ← spermidine/putrescine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,903,8900CT 6.8% 175.2 / 3.0 59P285P (CCG→CCA) spuHspermidine/putrescine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (26/29);  new base T (2/2);  total (28/31)
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

CTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGTTGA  >  NC_002947/5903752‑5904035
                                                                                                                                          |                                                                                                                                                 
cgggAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGGGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGaa                                                                                                                                          <  1:205185/146‑1 (MQ=255)
 tGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCGGAAGGCTGTGAGGTCCCTGTGGGCGCGGGCGTGCCCGCGCAAGGGCCAGGTCAGGCCGCCGCCGCCACCTGCGCCGGGCGCTTGCCCTGCGGCTGCGAAGc                                                                                                                                        >  6:109499/1‑149 (MQ=39)
   gAGCCTTCGCGGGCAGGCCCGCTCCCACAGGGCGCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGGGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGt                                                                                                                                      <  5:142791/149‑1 (MQ=255)
    aGCCTTCGCGGGCATGCCCGCTCCCACCGGGCTCACCACAACCCTGAAGGCTGTGAGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGGCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                                                                                                     <  5:200906/149‑1 (MQ=255)
    aGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGATGGCTGTGAGGTCCCTGTGGGAGGGGGGGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGtt                                                                                                                                     <  2:276319/149‑1 (MQ=255)
        ttCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCGTGACCTGCGGCTGCGAAGCGTTGGCt                                                                                                                                 >  8:301130/1‑149 (MQ=255)
           gcgGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgc                                                                                                                               <  5:109499/148‑1 (MQ=255)
           gcgGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCGGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTgc                                                                                                                               <  5:130147/148‑1 (MQ=255)
                     ccGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCctct                                                                                                                     <  2:298931/148‑1 (MQ=255)
                     ccGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCctct                                                                                                                     >  1:298931/1‑148 (MQ=255)
                     ccGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                                                                                                                    >  5:125865/1‑149 (MQ=255)
                     ccGCGCCCACCGGGCTCACCACACCCCTGAAGGCTGTGAGGGCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                                                                                                                    <  3:5093/149‑1 (MQ=255)
                         tCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGAt                                                                                                                >  4:239565/1‑149 (MQ=255)
                           ccACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATc                                                                                                               <  2:234989/148‑1 (MQ=255)
                            cacaGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGc                                                                                                             <  6:125865/149‑1 (MQ=255)
                             acaGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                                                                                   >  6:359487/1‑142 (MQ=255)
                             acaGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTc                                                                                                                   <  5:359487/142‑1 (MQ=255)
                                gggCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgc                                                                                                         <  3:205512/149‑1 (MQ=255)
                                    tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                                                                                                                    <  4:172305/134‑1 (MQ=255)
                                    tCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCtt                                                                                                                    >  3:172305/1‑134 (MQ=255)
                                       ccacAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCg                                                                                                  >  3:320354/1‑149 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGTGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgct                                                                                                        >  5:267987/1‑142 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGTGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgct                                                                                                        <  6:138869/142‑1 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGTGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgct                                                                                                        >  5:138869/1‑142 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGTGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCctgct                                                                                                        <  6:267987/142‑1 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                                                                                        <  7:204436/110‑1 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGc                                                                                                                                        >  8:204436/1‑110 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                                                                                                 >  3:304841/1‑149 (MQ=255)
                                        cacaACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGc                                                                                                 <  4:320354/149‑1 (MQ=255)
                                         acaACCCTGAAGGATGTGAGGTCCCTGTGGGAGCGGGAGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGcc                                                                                                >  6:92688/1‑149 (MQ=255)
                                           aaCCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCtt                                                                                              <  4:259363/149‑1 (MQ=255)
                                                 gAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac                                                                                         >  7:177143/1‑148 (MQ=255)
                                                 gAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGcgac                                                                                         <  8:177143/148‑1 (MQ=255)
                                                             gTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGc                                                                            >  5:362617/1‑149 (MQ=255)
                                                                ccTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTg                                                                         >  8:106614/1‑149 (MQ=255)
                                                                ccTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTg                                                                         <  6:362617/149‑1 (MQ=255)
                                                                    tgGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTggcgg                                                                     <  7:259147/149‑1 (MQ=255)
                                                                    tgGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTggc                                                                       <  7:106614/147‑1 (MQ=255)
                                                                      ggAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCt                                                                   >  8:151706/1‑149 (MQ=255)
                                                                                   ccGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGaa                                                                <  5:1860/139‑1 (MQ=255)
                                                                                   ccGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGaa                                                                >  6:1860/1‑139 (MQ=255)
                                                                                   ccGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGAc                                                       >  3:234187/1‑148 (MQ=255)
                                                                                   ccGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACg                                                      <  1:325241/149‑1 (MQ=255)
                                                                                   ccGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCCGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGAc                                                       >  2:218636/1‑148 (MQ=255)
                                                                                   ccGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCCGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGAc                                                       >  2:213103/1‑148 (MQ=255)
                                                                                         aaGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCTCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGaagaag                                                <  4:304841/149‑1 (MQ=255)
                                                                                            ggCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTc                                             <  5:92688/149‑1 (MQ=255)
                                                                                                aGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCAc                                           >  3:49952/1‑147 (MQ=255)
                                                                                                       gCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATAc                                  >  7:61626/1‑149 (MQ=255)
                                                                                                          gACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc                                                                                            <  6:330460/88‑1 (MQ=255)
                                                                                                          gACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGc                                                                                            >  5:330460/1‑88 (MQ=255)
                                                                                                                          ccGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTg               <  4:234187/149‑1 (MQ=255)
                                                                                                                                  ttGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGc       >  6:232058/1‑149 (MQ=255)
                                                                                                                                     aCCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGtt    >  3:163929/1‑149 (MQ=255)
                                                                                                                                     aCCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGt     <  4:49952/148‑1 (MQ=255)
                                                                                                                                     aCCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCCGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGtt    <  1:218636/149‑1 (MQ=255)
                                                                                                                                      ccTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGtt    >  1:85310/1‑148 (MQ=255)
                                                                                                                                      ccTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCCGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGtt    <  1:213103/148‑1 (MQ=255)
                                                                                                                                       cTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGTTGa  >  3:364591/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                 
CTGGAGCCTTCGCGGGCATGCCCGCTCCCACAGGGCTCACCACAACCCTGAAGGCTGTGAGGTCCCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGGTCAGGCCGACGCAGCCACCTGCGCCGGGCGCTTGACCTGCGGCTGCGAAGCGTTGGCTGCAGCACCCTCTTCGATCGCCTGCTGAATCGCCTTGCGACGACGCTCTTCGGCCTGGCGGCTGAAGTACCAGACGAAGAAGGTCACCAGCGATACCGACAGCAAGATCAGGCTGGCCACGGCGTTGA  >  NC_002947/5903752‑5904035

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