Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,409,923 C→G 6.7% intergenic (‑489/‑47) PP_3022 ← / → PP_3024 AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,409,9230CG 6.7% 192.5 / 3.4 60intergenic (‑489/‑47)PP_3022/PP_3024AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base G (2/2);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

CTATTCCTGGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATGTG  >  NC_002947/3409775‑3410071
                                                                                                                                                    |                                                                                                                                                     
cTATTCCTGGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTggcggc                                                                                                                                                       >  1:346915/1‑149 (MQ=255)
        ggCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAAcc                                                                                                                                               >  3:239325/1‑149 (MQ=255)
          ccATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAAccc                                                                                                                                              >  3:315850/1‑148 (MQ=255)
             tCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGc                                                                                                                                          <  1:263311/149‑1 (MQ=255)
               aaGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGa                                                                                                                                        >  4:338285/1‑149 (MQ=255)
                           gcgGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAcccag                                                                                                                                            <  6:98537/133‑3 (MQ=255)
                           gcgGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAcccag                                                                                                                                            >  5:98537/1‑131 (MQ=255)
                             ggCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAAc                                                                                                                          <  4:239325/149‑1 (MQ=255)
                                 tACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTa                                                                                                                      >  3:230928/1‑149 (MQ=255)
                                   ccTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAg                                                                                                                     <  7:114372/148‑1 (MQ=255)
                                         cTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAg                                                                                                              <  3:338285/148‑1 (MQ=255)
                                          tACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAggg                                                                                                            >  3:204486/1‑149 (MQ=255)
                                           aCATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAggg                                                                                                            >  4:181512/1‑148 (MQ=255)
                                            cATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCACCCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGt                                                                                                           >  7:217750/1‑148 (MQ=255)
                                             aTCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTa                                                                                                          >  4:114417/1‑148 (MQ=255)
                                              tCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAccca                                                                                                                                             >  1:308011/1‑112 (MQ=255)
                                              tCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAccca                                                                                                                                             <  2:308011/113‑2 (MQ=255)
                                                 ccTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTg                                                                                                     <  2:225991/149‑1 (MQ=255)
                                                 ccTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTg                                                                                                     >  1:249728/1‑149 (MQ=255)
                                                              gcgcgcACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTc                                                                                         <  5:135664/148‑1 (MQ=255)
                                                                      ccGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGc                                                                                    <  7:314388/145‑1 (MQ=255)
                                                                      ccGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGc                                                                                    >  8:314388/1‑145 (MQ=255)
                                                                      ccGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCAtt                                                                                 >  1:155054/1‑148 (MQ=255)
                                                                        gTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGAc                                                                              <  2:249728/149‑1 (MQ=255)
                                                                        gTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGAc                                                                              >  4:175796/1‑149 (MQ=255)
                                                                         tGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACt                                                                             >  2:130749/1‑149 (MQ=255)
                                                                         tGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACt                                                                             <  8:147669/149‑1 (MQ=255)
                                                                           aGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTc                                                                            <  2:346915/148‑1 (MQ=255)
                                                                             cGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTcagc                                                                         <  4:204486/149‑1 (MQ=255)
                                                                               agaCCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTcagcag                                                                       <  3:175796/149‑1 (MQ=255)
                                                                                 aCCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAgg                                                                      <  4:230928/148‑1 (MQ=255)
                                                                                 aCCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGc                                                                     >  8:268422/1‑149 (MQ=255)
                                                                                    gtggtgGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                  >  6:322749/1‑149 (MQ=255)
                                                                                        tggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                  <  6:7852/145‑1 (MQ=255)
                                                                                        tggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                  >  5:7852/1‑145 (MQ=255)
                                                                                         ggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACagag                                                             >  2:22937/1‑149 (MQ=255)
                                                                                          gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCCACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                            >  5:393064/1‑149 (MQ=255)
                                                                                          gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                            >  2:275088/1‑149 (MQ=255)
                                                                                          gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                            >  1:208116/1‑149 (MQ=255)
                                                                                          gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                            >  5:64767/1‑149 (MQ=255)
                                                                                            aaGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCCTCCATTGACTCAGCATGCGAACAGAGGc                                                           >  8:384243/1‑148 (MQ=255)
                                                                                               cGATCACCCAGGCGTCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAAt                                                       <  6:64767/149‑1 (MQ=255)
                                                                                                     gccAGGCGCCAAGTGCCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCAc                                                 >  7:76303/2‑149 (MQ=255)
                                                                                                      ccAGGCGGCACTTGGCGCCTGGCGCGACAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACg                                                >  2:118693/1‑149 (MQ=255)
                                                                                                       ctggcggcCCTTTGCGCCTGGCGCGCAAGGAGTTCTGGGTTGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGc                                               <  4:315850/147‑1 (MQ=255)
                                                                                                            ggcACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTc                                          >  7:104516/1‑149 (MQ=255)
                                                                                                                   ggCGCCTGGCTCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTcgtacgta                                   >  4:21279/1‑149 (MQ=255)
                                                                                                                   ggCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACt                                     <  3:114417/147‑1 (MQ=255)
                                                                                                                    gcgcCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTa                                   >  1:194223/1‑148 (MQ=255)
                                                                                                                           gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCt                                                                                          >  7:6337/1‑86 (MQ=255)
                                                                                                                           gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCt                                                                                          <  8:6337/86‑1 (MQ=255)
                                                                                                                           gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATc                             <  2:155054/147‑1 (MQ=255)
                                                                                                                           gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCa                            <  6:393064/148‑1 (MQ=255)
                                                                                                                                    gAGTTCTGGGTGGCGGCCGGCACCCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCCCGCTGGTCGTACTTTCAAATCACCTGCGCagt                  >  4:283578/1‑147 (MQ=255)
                                                                                                                                    gAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCAtt                  <  1:22937/149‑1 (MQ=255)
                                                                                                                                          tGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCTTGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATgg            <  7:384243/149‑1 (MQ=255)
                                                                                                                                          tGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATg             >  4:350340/1‑148 (MQ=255)
                                                                                                                                            ggTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCATTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCACATTGAAt              >  1:82905/1‑146 (MQ=255)
                                                                                                                                            ggTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCACATTGAAt              <  2:82905/145‑1 (MQ=255)
                                                                                                                                            ggTGGCGGCAGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGat          <  8:104516/149‑1 (MQ=255)
                                                                                                                                              tGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGata         >  2:133317/1‑148 (MQ=255)
                                                                                                                                                  ggccggcAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAAAAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCAtg    >  6:123467/1‑149 (MQ=255)
                                                                                                                                                    ccggcAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCAtgtg  <  1:133317/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTATTCCTGGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATGTG  >  NC_002947/3409775‑3410071

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: