Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,028 G→C 7.4% H64Q (CAC→CAG) ‡ PP_4580 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,0280GC 7.4% 195.0 / 6.6 68H64Q (CAC→CAG) ‡PP_4580hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/39);  new base C (3/2);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

GTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199881‑5200169
                                                                                                                                                   |                                                                                                                                              
gTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGt                                                                                                                                               <  7:173661/149‑1 (MQ=255)
 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)
                                                                                   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                         <  6:182712/148‑1 (MQ=255)
                                                                                                                      ttGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGa                         <  4:25755/148‑1 (MQ=255)
                                                                                                                        gAACACCGGGTCGTCGCTTGTGCGGTAGTGCTGCAGCACCGCCATCAAGTACCCCGTTAAAAAACGCTCGGCCGTTTGCCGGCACCCAACGAGCGCCACGCGCGAATCGA‑AAAGGTGCGCCTCGACCCGCTGGCTATAGTcaatgaac                       >  8:15386/1‑141 (MQ=255)
                                                                                                                         aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          >  3:197823/1‑144 (MQ=255)
                                                                                                                         aaCACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTg                          <  4:197823/144‑1 (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          <  3:42223/148‑1 (MQ=255)
                                                                                                                                     tcgCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCt          >  4:42223/1‑148 (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    <  3:23979/148‑1 (MQ=255)
                                                                                                                                           ctgcGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATAc    <  7:15386/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                        >  4:100450/1‑121 (MQ=255)
                                                                                                                                                  aGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGAc                        <  3:100450/121‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199881‑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: ‑