Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,030 G→A 8.0% H64Y (CAC→TAC) ‡ PP_4580 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,0300GA 8.0% 164.3 / 4.6 62H64Y (CAC→TAC) ‡PP_4580hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/35);  new base A (3/2);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

TCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199882‑5200169
                                                                                                                                                    |                                                                                                                                            
tCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtg                                                                                                                                              <  4:219475/149‑1 (MQ=255)
  ggCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgct                                                                                                                                            <  3:57198/149‑1 (MQ=255)
   gCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctg                                                                                                                                           >  3:239522/1‑149 (MQ=255)
   gCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctg                                                                                                                                           <  3:146394/149‑1 (MQ=255)
    ccTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGtgctgc                                                                                                                                          <  6:4929/149‑1 (MQ=255)
     cTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTTGGGTAgtggtg                                                                                                                                           >  5:245020/1‑147 (MQ=255)
          ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCagcagc                                                                                                                                    <  5:96859/149‑1 (MQ=255)
          ccAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCAGCACTACCGcagcag                                                                                                                                     <  5:183497/148‑1 (MQ=255)
           cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcg                                                                                                                                   >  4:283209/1‑149 (MQ=255)
           cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcg                                                                                                                                   >  5:139834/1‑149 (MQ=255)
                    ccAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGGGGTTGTGCTGCAGCAGCGCCTGCaact                                                                                                                          >  5:19569/1‑147 (MQ=255)
                      aGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCg                                                                                                                        <  2:222168/149‑1 (MQ=255)
                      aGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCAGTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGGTGCGGGTGTGCTGGAGCGGCGGCTggacttgg                                                                                                                        >  8:81886/1‑142 (MQ=255)
                        gTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGcc                                                                                                                      <  3:283209/149‑1 (MQ=255)
                         tACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGccc                                                                                                                     <  4:239522/149‑1 (MQ=255)
                             cccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTg                                                                                                                  >  7:91747/1‑148 (MQ=255)
                              ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTg                                                                                                                               <  3:155767/134‑1 (MQ=255)
                              ccGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGa                                                                                                                 >  1:187936/1‑148 (MQ=255)
                                 gCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGaaaaa                                                                                                             <  2:187936/149‑1 (MQ=255)
                                          ccTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGc                                                                                                    >  2:219219/1‑149 (MQ=255)
                                              gcagcaCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATAATCAGGCTTTCGATGCTGTTGGCATAACTGTTGAACACCGGGTCGTCGCTGCTGGGGTATTGCTGCAGACCCGCCTTCAATTCTCCCATGAAAAAACGGTCGGACggt                                                                                                >  3:289070/1‑147 (MQ=255)
                                                agcaCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTg                                                                                              <  8:13830/149‑1 (MQ=255)
                                                  caCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGcc                                                                                            <  5:73218/149‑1 (MQ=255)
                                                   aCCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCg                                                                                           <  8:100836/149‑1 (MQ=255)
                                                      tGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        >  1:51829/1‑149 (MQ=255)
                                                        ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  6:139834/147‑1 (MQ=255)
                                                        ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  6:153133/147‑1 (MQ=255)
                                                        ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAAc                                                                                       <  8:91747/148‑1 (MQ=255)
                                                        ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAAc                                                                                       <  1:219219/148‑1 (MQ=255)
                                                        ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTATCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  5:261584/147‑1 (MQ=255)
                                                               ttttCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAg                                                                               <  2:51829/149‑1 (MQ=255)
                                                                  tCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAgcg                                                                             <  6:37080/148‑1 (MQ=255)
                                                                   ccAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGAGCCTGCAAGTAGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCAACGAccgc                                                                            >  6:25864/1‑148 (MQ=255)
                                                                    cagcagTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCt                                                                          >  2:197782/1‑149 (MQ=255)
                                                                      gcagTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCTCTGTTTAACACCGGGTCGTCGCTGGTTCGGTAGTGGAGCAGCAGCGGCCGGAAGTCGCCCTTTAAAAACCGCTAGGCCGTTTTCCGGAACCCAAGCAGCGCCTcg                                                                        >  1:59013/1‑149 (MQ=255)
                                                                                  ggATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAAc                                                           <  1:147393/149‑1 (MQ=255)
                                                                                    atcatcAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAg                                                         <  7:81886/149‑1 (MQ=255)
                                                                                              tttCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCCTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTc                                                >  6:310150/1‑148 (MQ=255)
                                                                                                  gATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGAccc                                           >  3:238345/1‑149 (MQ=255)
                                                                                                   aTGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCg                                          >  2:55560/1‑149 (MQ=255)
                                                                                                           ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCta                                   <  7:16138/148‑1 (MQ=255)
                                                                                                           ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCta                                   <  2:59013/148‑1 (MQ=255)
                                                                                                           ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAACGCTCGGCCGTTTTCCGGAACCCATAGAGCGCCTCGCGCGAATCGAACAAGTGCGCCTCGACCCCCTGGCta                                   >  8:16138/1‑148 (MQ=255)
                                                                                                              aTCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGt                               >  3:184716/1‑149 (MQ=255)
                                                                                                                  cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGATACAGGTGCGCCTCGAACCGCTGGCTATAGTCCt                            <  1:197782/149‑1 (MQ=255)
                                                                                                                     ttGAACACCGGGTCGTCGCTGCTGCGGTATTGCTGCAGCAGCGACTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        <  8:218232/149‑1 (MQ=255)
                                                                                                                     ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  4:25755/148‑1 (MQ=255)
                                                                                                                     ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  6:182712/148‑1 (MQ=255)
                                                                                                                       gAACACCGGGTCGTCGCTTGTGCGGTAGTGCTGCAGCACCGCCATCAAGTACCCCGTTAAAAAACGCTCGGCCGTTTGCCGGCACCCAACGAGCGCCACGCGCGAATCGAAAAGGTGCGCCTCGACCCGCTGGCTATAGTcaatgaac                       >  8:15386/1‑141 (MQ=255)
                                                                                                                        aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          <  4:197823/144‑1 (MQ=255)
                                                                                                                        aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          >  3:197823/1‑144 (MQ=255)
                                                                                                                         acacCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGc                    <  4:184716/149‑1 (MQ=255)
                                                                                                                                 tcgtcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc            <  6:245020/149‑1 (MQ=255)
                                                                                                                                 caggcgCTGCTGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc            >  1:291587/5‑149 (MQ=255)
                                                                                                                                    tcgCTGCTGCGGTAGTGCTGCAGCAGGGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg         <  5:25864/149‑1 (MQ=255)
                                                                                                                                    tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          >  4:42223/1‑148 (MQ=255)
                                                                                                                                    tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          <  3:42223/148‑1 (MQ=255)
                                                                                                                                    tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg         >  3:129517/1‑149 (MQ=255)
                                                                                                                                      gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg               <  7:152340/141‑1 (MQ=255)
                                                                                                                                      gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg               >  8:152340/1‑141 (MQ=255)
                                                                                                                                      gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa       >  7:185555/1‑149 (MQ=255)
                                                                                                                                          ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  7:15386/148‑1 (MQ=255)
                                                                                                                                          ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  3:23979/148‑1 (MQ=255)
                                                                                                                                             cGGTAGTGCTGCAGCATCGCCTGCAAGTCGCAGTTGAAATAGCGCTCGGCCGTGTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAcc   <  7:206389/146‑1 (MQ=255)
                                                                                                                                             cGGTAGTGCTGCAGCAGCGCCTGCAAGTCTCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCTCCTCGCGCGATTCGAACAGGTTCGCCTCGACCCGCTGTCTATAGTCATTGACCAGCCCCGGGGCCTGGATACCg  <  4:289070/147‑1 (MQ=255)
                                                                                                                                             cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg  <  6:19569/147‑1 (MQ=255)
                                                                                                                                                 aGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        <  3:100450/121‑1 (MQ=255)
                                                                                                                                                 aGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        >  4:100450/1‑121 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGAACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199882‑5200169

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: