Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,032 T→C 8.2% Q63R (CAG→CGG)  PP_4580 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,0320TC 8.2% 167.3 / 5.2 62Q63R (CAG→CGG) PP_4580hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (21/35);  new base C (3/2);  total (25/37)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

GGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199884‑5200169
                                                                                                                                                    |                                                                                                                                          
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)
                                                                                ggATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑Ac                                                           <  1:147393/149‑1 (MQ=255)
                                                                                  atcatcAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAg                                                         <  7:81886/149‑1 (MQ=255)
                                                                                            tttCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCCTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTc                                                >  6:310150/1‑148 (MQ=255)
                                                                                                gATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGAccc                                           >  3:238345/1‑149 (MQ=255)
                                                                                                 aTGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCg                                          >  2:55560/1‑149 (MQ=255)
                                                                                                         ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCta                                   <  2:59013/148‑1 (MQ=255)
                                                                                                         ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCta                                   <  7:16138/148‑1 (MQ=255)
                                                                                                         ggCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAACGCTCGGCCGTTTTCCGGAACCCATAGAGCGCCTCGCGCGAATCGA‑ACAAGTGCGCCTCGACCCCCTGGCta                                   >  8:16138/1‑148 (MQ=255)
                                                                                                            aTCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGt                               >  3:184716/1‑149 (MQ=255)
                                                                                                                cTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGATACAGGTGCGCCTCGAACCGCTGGCTATAGTCCt                            <  1:197782/149‑1 (MQ=255)
                                                                                                                   ttGAACACCGGGTCGTCGCTGCTGCGGTATTGCTGCAGCAGCGACTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        <  8:218232/149‑1 (MQ=255)
                                                                                                                   ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  4:25755/148‑1 (MQ=255)
                                                                                                                   ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  6:182712/148‑1 (MQ=255)
                                                                                                                     gAACACCGGGTCGTCGCTTGTGCGGTAGTGCTGCAGCACCGCCATCAAGTACCCCGTTAAAAAACGCTCGGCCGTTTGCCGGCACCCAACGAGCGCCACGCGCGAATCGA‑AAAGGTGCGCCTCGACCCGCTGGCTATAGTcaatgaac                       >  8:15386/1‑141 (MQ=255)
                                                                                                                      aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          <  4:197823/144‑1 (MQ=255)
                                                                                                                      aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          >  3:197823/1‑144 (MQ=255)
                                                                                                                       acacCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGc                    <  4:184716/149‑1 (MQ=255)
                                                                                                                               tcgtcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc            <  6:245020/149‑1 (MQ=255)
                                                                                                                               caggcgCTGCTGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGAc            >  1:291587/5‑149 (MQ=255)
                                                                                                                                  tcgCTGCTGCGGTAGTGCTGCAGCAGGGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg         <  5:25864/149‑1 (MQ=255)
                                                                                                                                  tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          >  4:42223/1‑148 (MQ=255)
                                                                                                                                  tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          <  3:42223/148‑1 (MQ=255)
                                                                                                                                  tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTg         >  3:129517/1‑149 (MQ=255)
                                                                                                                                    gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg               <  7:152340/141‑1 (MQ=255)
                                                                                                                                    gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCgg               >  8:152340/1‑141 (MQ=255)
                                                                                                                                    gctgctGCAGCACTACCGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGa       >  7:185555/1‑149 (MQ=255)
                                                                                                                                        ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  7:15386/148‑1 (MQ=255)
                                                                                                                                        ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  3:23979/148‑1 (MQ=255)
                                                                                                                                           cGGTAGTGCTGCAGCATCGCCTGCAAGTCGCAGTTGAAATAGCGCTCGGCCGTGTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAcc   <  7:206389/146‑1 (MQ=255)
                                                                                                                                           cGGTAGTGCTGCAGCAGCGCCTGCAAGTCTCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCTCCTCGCGCGATTCGA‑ACAGGTTCGCCTCGACCCGCTGTCTATAGTCATTGACCAGCCCCGGGGCCTGGATACCg  <  4:289070/147‑1 (MQ=255)
                                                                                                                                           cGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCg  <  6:19569/147‑1 (MQ=255)
                                                                                                                                               aGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        <  3:100450/121‑1 (MQ=255)
                                                                                                                                               aGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        >  4:100450/1‑121 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
GGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199884‑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: ‑