Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,409,921 G→C 6.7% intergenic (‑487/‑49) 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,9210GC 6.7% 191.2 / 3.6 60intergenic (‑487/‑49)PP_3022/PP_3024AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/24);  new base C (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 = 8.90e-01

CCTATTCCTGGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATG  >  NC_002947/3409774‑3410069
                                                                                                                                                   |                                                                                                                                                     
ccTATTCCTGGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTggcgg                                                                                                                                                      >  1:225991/1‑149 (MQ=255)
 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                                                                                                                                           <  2:308011/113‑2 (MQ=255)
                                               tCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAccca                                                                                                                                           >  1:308011/1‑112 (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                                                                            >  4:175796/1‑149 (MQ=255)
                                                                         gTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGAc                                                                            <  2:249728/149‑1 (MQ=255)
                                                                          tGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACt                                                                           <  8:147669/149‑1 (MQ=255)
                                                                          tGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACt                                                                           >  2:130749/1‑149 (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                                                                >  5:7852/1‑145 (MQ=255)
                                                                                         tggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                <  6:7852/145‑1 (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)
                                                                                                                                                   |                                                                                                                                                     
CCTATTCCTGGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCTTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATG  >  NC_002947/3409774‑3410069

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: