Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,200,026 T→C 7.2% Y65C (TAC→TGC)  PP_4580 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,200,0260TC 7.2% 198.1 / 4.9 70Y65C (TAC→TGC) PP_4580hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/41);  new base C (3/2);  total (27/43)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

GGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199879‑5200169
                                                                                                                                                   |                                                                                                                                                
gggTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTa                                                                                                                                                 <  5:73380/149‑1 (MQ=255)
gggTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTa                                                                                                                                                 <  7:42926/149‑1 (MQ=255)
  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)
      gCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATAATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTAGCATCTTCGGTaattaat                                                                                                                                           >  6:261584/1‑143 (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                                                                                                                                   >  5:139834/1‑149 (MQ=255)
              cAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAgcg                                                                                                                                   >  4:283209/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:153133/147‑1 (MQ=255)
                                                           ccATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGaa                                                                                        <  6:139834/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          >  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    <  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)
                                                                                                                                                   |                                                                                                                                                
GGGTCGGCCTGTACCAGCACATGCCAGGTACGCCCGGCCCATGACCCTTGCAGCACCTGCCATTCGTTTTCCAGCAGTTCGTTGCGGATCATCAGGCTTTCGATGCTGTTGGCATCACTGTTGAACACCGGGTCGTCGCTGCTGCGGTAGTGCTGCAGCAGCGCCTGCAAGTCGCCCTTGAAAAAGCGCTCGGCCGTTTGCCGGAACCCATCGAGCGCCTCGCGCGATTCGA‑ACAGGTGCGCCTCGACCCGCTGGCTATAGTCCTTGACCAGCCCCGGGACCTGGATACCG  >  NC_002947/5199879‑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: ‑