Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,577,639 T→C 6.0% G18G (GGT→GGC)  PP_3158 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,577,6390TC 6.0% 199.8 / 2.6 68G18G (GGT→GGC) PP_3158hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base C (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CATTCTCTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTG  >  NC_002947/3577492‑3577783
                                                                                                                                                      |                                                                                                                                                
cATTCTCTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTg                                                                                                                                                 >  8:72046/1‑149 (MQ=255)
  ttCTCTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGcc                                                                                                                                               <  8:219706/149‑1 (MQ=255)
  ttCTCTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCACCGCCTGCTACTG‑‑‑ACAGGTGcc                                                                                                                                               >  7:376579/1‑149 (MQ=255)
      ctcCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTgg                                                                                                                                           >  7:326427/1‑149 (MQ=255)
       tcCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTgg                                                                                                                                           >  5:61465/1‑148 (MQ=255)
         cAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                        <  6:169346/149‑1 (MQ=255)
          aGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                        >  4:278657/1‑148 (MQ=255)
          aGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                        <  3:278657/148‑1 (MQ=255)
           gCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGc                                                                                                                                      <  8:326427/149‑1 (MQ=255)
           gCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGc                                                                                                                                      <  8:109335/149‑1 (MQ=255)
           gCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGc                                                                                                                                      <  2:166208/149‑1 (MQ=255)
            ccACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                        <  1:388858/146‑1 (MQ=255)
            ccACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCAt                                                                                                                                        >  2:388858/1‑146 (MQ=255)
             cACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCt                                                                                                                                    <  1:186146/149‑1 (MQ=255)
                cGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCtgtg                                                                                                                                 >  1:109937/1‑149 (MQ=255)
                   aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                              >  7:179686/1‑149 (MQ=255)
                   aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                              >  6:20488/1‑149 (MQ=255)
                       tttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTg                                                                                                                           >  7:105352/1‑148 (MQ=255)
                        ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGt                                                                                                                          <  8:179686/148‑1 (MQ=255)
                        ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTa                                                                                                                         >  6:160556/1‑149 (MQ=255)
                         tCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTAc                                                                                                                        >  3:244080/1‑149 (MQ=255)
                           ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCTACCCTGTTCCGATGGTTGCAAGGTACTCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                       <  8:295010/148‑1 (MQ=255)
                           ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                       >  8:270317/1‑148 (MQ=255)
                            cATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGCGCACAGGCATGCCAGGGCACCTGtcagta                                                                                                                     <  6:219781/149‑6 (MQ=38)
                             aTTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACgcag                                                                                                                    <  1:266494/149‑1 (MQ=255)
                                   gggCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGcagcagcag                                                                                                              <  2:109937/149‑1 (MQ=255)
                                      ctctATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGc                                                                                                           >  4:114196/1‑149 (MQ=255)
                                         tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGcc                                                                                                          <  6:256792/147‑1 (MQ=255)
                                         tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAAGTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                        >  1:87417/1‑149 (MQ=255)
                                           tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCg                                                                                                         <  2:271396/146‑1 (MQ=255)
                                           tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                        <  7:39925/147‑1 (MQ=255)
                                            gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                             <  4:137270/144‑5 (MQ=18)
                                            gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                             >  3:137270/1‑140 (MQ=18)
                                                 gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATgg                                                                                                 <  4:244080/148‑1 (MQ=255)
                                                 gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                <  6:61465/149‑1 (MQ=255)
                                                      gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                <  4:185055/144‑1 (MQ=255)
                                                      gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                >  3:185055/1‑144 (MQ=255)
                                                            cTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcg                                                                                     <  7:270317/149‑1 (MQ=255)
                                                                 tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCGTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                >  7:218030/1‑149 (MQ=255)
                                                                 tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                >  4:183299/1‑149 (MQ=255)
                                                                   ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                              >  6:242441/1‑149 (MQ=255)
                                                                   ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacg                                                                               <  5:160556/148‑1 (MQ=255)
                                                                      gcgcAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacggca                                                                            >  1:379792/1‑148 (MQ=255)
                                                                          aGCCGCCCAATATCGACCCGTAGCCGAAGGTAGCAAGCTACGCGAATTTCCTACCCCGCCTGCTACTG‑‑‑ACAGTTGCCCTGGCATGCCTGTGCGCCTTTACGCATCAGCAGGTCCGTGACATGGTCAACCAGCTCGGCAACGTCAAAgcc                                                                       <  7:166119/149‑1 (MQ=255)
                                                                                  aaTATCGACCCGGTGCCGATGGTTGCAAGTTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCt                                                               >  6:215466/1‑149 (MQ=255)
                                                                                     atCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGtt                                                            >  6:374055/1‑149 (MQ=255)
                                                                                        gACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCa                                                         <  8:105352/149‑1 (MQ=255)
                                                                                          cccGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGa                                                       <  2:379792/149‑1 (MQ=255)
                                                                                                gCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCg                                                 <  7:72046/149‑1 (MQ=255)
                                                                                                         tGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                              <  2:281141/111‑1 (MQ=255)
                                                                                                           cAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTg                                      <  3:114196/149‑1 (MQ=255)
                                                                                                            aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                     <  5:215466/149‑1 (MQ=255)
                                                                                                            aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                     >  2:102131/1‑149 (MQ=255)
                                                                                                                  cgcgAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGc                               >  1:1461/1‑149 (MQ=255)
                                                                                                                      aaTTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTc                           <  7:120340/149‑1 (MQ=255)
                                                                                                                             tccCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                    >  6:15830/1‑149 (MQ=255)
                                                                                                                               ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCTCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGTCCGCCTGGCCACGCTGTCGATgcgc                    <  8:138593/147‑1 (MQ=255)
                                                                                                                               ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                    <  8:218030/147‑1 (MQ=255)
                                                                                                                               ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGcgcg                   >  4:365268/1‑148 (MQ=255)
                                                                                                                                     cTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCt            >  1:93426/1‑149 (MQ=255)
                                                                                                                                        ccaCTG‑‑‑ACCCGGGCCCTTGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGGGACATGGGCAACCACCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa         <  5:144422/147‑1 (MQ=255)
                                                                                                                                         tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa         >  2:128631/1‑148 (MQ=255)
                                                                                                                                         tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGACTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACAGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa         >  4:93652/1‑148 (MQ=255)
                                                                                                                                              cgtaCAGGCGCACAGGCATGCCAGGGCACCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAg      >  7:43058/4‑149 (MQ=255)
                                                                                                                                                 aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                    <  4:69522/84‑1 (MQ=255)
                                                                                                                                                 aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                    >  3:69522/1‑84 (MQ=255)
                                                                                                                                                 aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGc     <  2:1461/147‑1 (MQ=255)
                                                                                                                                                  cAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGGCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTg  <  8:376579/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
CATTCTCTCCAGCCACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTG  >  NC_002947/3577492‑3577783

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