Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,996,900 A→C 5.7% M88R (ATG→AGG)  PP_3525 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,996,9000AC 5.7% 222.7 / 3.7 70M88R (ATG→AGG) PP_3525hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (37/29);  new base C (2/2);  total (39/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.93e-01

CAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGG  >  NC_002947/3996752‑3997044
                                                                                                                                                    |                                                                                                                                                
cAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCa                                                                                                                                                  >  3:339618/1‑149 (MQ=255)
     cGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATgcgc                                                                                                                                             >  6:351275/1‑149 (MQ=255)
      gCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCt                                                                                                                                            <  5:351275/149‑1 (MQ=255)
        ccTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTg                                                                                                                                           >  3:382272/1‑148 (MQ=255)
                  ttGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAGCAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCa                                                                                                                                >  5:356460/1‑149 (MQ=255)
                    gCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCt                                                                                                                                            <  1:110968/135‑1 (MQ=255)
                    gCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCt                                                                                                                                            >  2:110968/1‑135 (MQ=255)
                      aaCCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGaa                                                                                                                            >  2:90309/1‑149 (MQ=255)
                          aCTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACa                                                                                                                        <  3:226305/149‑1 (MQ=255)
                             gCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAgcg                                                                                                                     >  2:131744/1‑149 (MQ=255)
                                  tGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGcgcg                                                                                                                >  2:238804/1‑149 (MQ=255)
                                  tGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGcgcg                                                                                                                >  1:260068/1‑149 (MQ=255)
                                    gCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGt                                                                                                               <  2:260068/148‑1 (MQ=255)
                                      cctcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCt                                                                                                            >  6:228778/1‑149 (MQ=255)
                                      cctcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCt                                                                                                            >  7:266260/1‑149 (MQ=255)
                                        tcctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTc                                                                                                          >  7:118443/1‑149 (MQ=255)
                                          ctGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTcc                                                                                                         <  6:356460/148‑1 (MQ=255)
                                           tGCCCGGGGTACCGTCTGGCCACTCGCCCATCACGGTCGAGCAGCAAACGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAg                                                                                                       >  6:8626/1‑149 (MQ=255)
                                             cccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                      >  4:32158/1‑148 (MQ=255)
                                             cccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAg                                                                                                       <  4:339618/147‑1 (MQ=255)
                                              ccGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAgg                                                                                                      >  6:418686/1‑147 (MQ=255)
                                               cGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTg                                                                                                                                           >  5:192559/1‑109 (MQ=255)
                                               cGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTg                                                                                                                                           <  6:192559/109‑1 (MQ=255)
                                               cGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAg                                                                                                   <  5:418686/149‑1 (MQ=255)
                                                      ccGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTgg                                                                                            <  3:32158/149‑1 (MQ=255)
                                                         gCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCtt                                                                                              >  2:72960/1‑144 (MQ=255)
                                                         gCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCtt                                                                                              <  1:72960/144‑1 (MQ=255)
                                                         gCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGgct                                                                                         <  8:118443/149‑1 (MQ=255)
                                                             gCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCtt                                                                                              <  1:316193/140‑1 (MQ=255)
                                                             gCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCtt                                                                                              >  2:316193/1‑140 (MQ=255)
                                                               cACTCGCCCATCACGGTCGAGCAGCACCGGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCgg                                                                                   >  1:232888/1‑149 (MQ=255)
                                                               cACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTcc                                                                                                         <  4:106630/127‑1 (MQ=255)
                                                               cACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTcc                                                                                                         >  3:106630/1‑127 (MQ=255)
                                                               cACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCgg                                                                                   >  2:71719/1‑149 (MQ=255)
                                                                  tCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGc                                                                                  >  5:344586/1‑147 (MQ=255)
                                                                      ccATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATa                                                                                                                          <  4:107758/103‑1 (MQ=255)
                                                                      ccATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATa                                                                                                                          >  3:107758/1‑103 (MQ=255)
                                                                      ccATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTgc                                                                            >  7:255534/1‑149 (MQ=255)
                                                                       cATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTgcc                                                                           <  8:266260/149‑1 (MQ=255)
                                                                           aCGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTgccgcc                                                                        <  4:382272/148‑1 (MQ=255)
                                                                            cGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGTGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCtt                                                                      <  1:90309/149‑1 (MQ=255)
                                                                               tCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTTCTCGGCCAGTGCCGCCTTGAc                                                                   <  5:228778/149‑1 (MQ=255)
                                                                                      gcaCCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTg                                                            <  2:253711/149‑1 (MQ=255)
                                                                                       caCCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTgg                                                           >  8:262987/1‑149 (MQ=255)
                                                                                           cGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCtt                                                                                                           >  7:147339/1‑97 (MQ=255)
                                                                                           cGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCtt                                                                                                           <  8:147339/97‑1 (MQ=255)
                                                                                             gTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                           >  5:96824/1‑111 (MQ=255)
                                                                                             gTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTggg                                                                                           <  6:96824/111‑1 (MQ=255)
                                                                                                     aaTCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcg                                             >  3:103993/1‑149 (MQ=255)
                                                                                                      aTCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgc                                            >  2:333582/1‑149 (MQ=255)
                                                                                                         gcgcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATcgcgcg                                         >  6:318015/1‑149 (MQ=255)
                                                                                                           gcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGa                                                                    >  8:322320/1‑120 (MQ=255)
                                                                                                           gcATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGa                                                                    <  7:322320/120‑1 (MQ=255)
                                                                                                             aTGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                     >  5:122979/1‑149 (MQ=255)
                                                                                                             aTGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                     <  6:344586/149‑1 (MQ=255)
                                                                                                               ggCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                     >  6:109430/1‑147 (MQ=255)
                                                                                                                       cctcgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCTCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATc                           <  4:103993/146‑1 (MQ=255)
                                                                                                                         gcgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                     <  4:199274/137‑1 (MQ=255)
                                                                                                                         gcgAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGccac                                     >  3:199274/1‑137 (MQ=255)
                                                                                                                          cgAACAGTTTGCTGATCAGCGCATGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGc                        <  8:255534/149‑1 (MQ=255)
                                                                                                                           gAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGa                                                                    <  5:420521/104‑1 (MQ=255)
                                                                                                                           gAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGa                                                                    >  6:420521/1‑104 (MQ=255)
                                                                                                                           gAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCg                       >  2:124251/1‑149 (MQ=255)
                                                                                                                              cAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCgtgt                    <  1:71719/149‑1 (MQ=255)
                                                                                                                                    gCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAc              >  8:331133/1‑149 (MQ=255)
                                                                                                                                      cattCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTg                                                            <  6:240657/98‑1 (MQ=255)
                                                                                                                                      cattCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTg                                                            >  5:240657/4‑101 (MQ=255)
                                                                                                                                       attCAGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCAcc             >  8:52352/3‑147 (MQ=255)
                                                                                                                                           aGCGCATGCCTGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCt       <  7:52352/149‑1 (MQ=255)
                                                                                                                                             cgcAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGt     <  5:109430/149‑1 (MQ=255)
                                                                                                                                                aGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAgg  >  2:2122/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CAGCACGCCCTGCTGCAGTTGCAACCACTGCACCTGGCCCTCCTGCCCGGGGTACCGGCTGGCCACTCGCCCATCACGGTCGAGCAGCACCCGGTAGCTGTAATCGCGCATGGCCGGGACCGCGAACAGTTTGCTGATCAGCGCAGGCATGCGCTGAATGTCGGCCACGAATACAGCGGCGCGTGCTTCCAGGTAGCCCTTGGGCTGCTCGGCCAGTGCCGCCTTGACCAGCTTGGCACCGTCCATATCGCGCGCCACCAGTAAAATCTGCGTGTCGGCACCAAGGCTGTAGG  >  NC_002947/3996752‑3997044

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