Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,364,175 C→A 8.3% A238D (GCC→GAC)  PP_2960 → transcriptional regulator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,364,1750CA 8.3% 180.0 / 6.9 60A238D (GCC→GAC) PP_2960transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.35e-01

TGGCGGAACTGGCGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAAT  >  NC_002947/3364028‑3364316
                                                                                                                                                   |                                                                                                                                             
tGGCGGAACTGGCGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGcc                                                                                                                                              >  4:173508/1‑149 (MQ=255)
tGGCGGAACTGGCGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGcc                                                                                                                                              >  8:253521/1‑149 (MQ=255)
tGGCGGAACTGGCGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGcc                                                                                                                                              <  8:192164/149‑1 (MQ=255)
          ggCGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGcc                                                                                                                                     >  2:199281/1‑148 (MQ=255)
            cGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGAt                                                                                                                                  <  8:211177/149‑1 (MQ=255)
                  tttCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCGATTGCCCTGATCGCCGATGAGGTc                                                                                                                            >  4:172197/1‑149 (MQ=255)
                   ttCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCg                                                                                                                           >  6:135388/1‑149 (MQ=255)
                    tCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTc                                                                                                                            <  7:271480/147‑1 (MQ=255)
                    tCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCg                                                                                                                           <  3:173508/148‑1 (MQ=255)
                     cGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGACCTGATCGCCGATGAGGTCgg                                                                                                                          <  1:213150/148‑1 (MQ=255)
                       ggCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTa                                                                                                                       >  7:43930/1‑149 (MQ=255)
                             tCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGaaa                                                                                                                  <  4:160820/148‑1 (MQ=255)
                             tCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAg                                                                                                                 >  3:326558/1‑149 (MQ=255)
                               gcgGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTc                                                                                                               <  4:326558/149‑1 (MQ=255)
                                   gCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTc                                                                                                           <  5:135388/149‑1 (MQ=255)
                                    cAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCt                                                                                                          <  7:253521/149‑1 (MQ=255)
                                    cAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGACTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGAGGCAAAAACGCCTTCGCGTAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCt                                                                                                          <  7:61128/149‑1 (MQ=255)
                                       gcttcgcttcGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTgcg                                                                                                       <  1:199281/149‑1 (MQ=255)
                                       gcttcgcttcGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACACAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGAt                                                                                                                                  >  7:154035/1‑122 (MQ=255)
                                       gcttcgcttcGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACACAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGAt                                                                                                                                  <  8:154035/122‑1 (MQ=255)
                                              ttcgcACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGccc                                                                                                >  1:17007/1‑149 (MQ=255)
                                              ttcgcACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGACCTGATCGCCGATGAGGTCGGCTACGAAAGTCCtt                                                                                                            >  6:145631/1‑137 (MQ=255)
                                              ttcgcACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGACCTGATCGCCGATGAGGTCGGCTACGAAAGTCCtt                                                                                                            <  5:145631/137‑1 (MQ=255)
                                                 gcACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGaa                                                                                             >  1:46051/1‑149 (MQ=255)
                                                 gcACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGaa                                                                                             >  3:72201/1‑149 (MQ=255)
                                                        cgcgAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTc                                                                                       <  2:17007/148‑1 (MQ=255)
                                                           gagaTCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTcgccg                                                                                   >  6:260273/1‑149 (MQ=255)
                                                           gagaTCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGACCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTcgccg                                                                                   >  4:142265/1‑149 (MQ=255)
                                                                 gTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAggt                                                                             <  2:46051/149‑1 (MQ=255)
                                                                     gTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAggtgggt                                                                         <  8:43930/149‑1 (MQ=255)
                                                                        ttACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGc                                                                       >  6:307848/1‑148 (MQ=255)
                                                                          aCGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGa                                                                     <  4:72201/148‑1 (MQ=255)
                                                                              cTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGccc                                                                 >  6:153643/1‑148 (MQ=255)
                                                                               tGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAGCACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGCAAGTCCTTCTGCGCTAGCCCGAACCGTGCGCCGAAAGGTGCGTACGAGCCCg                                                                >  6:252334/1‑148 (MQ=255)
                                                                                 cGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGaa                                                              <  5:153643/148‑1 (MQ=255)
                                                                                  ggACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCGATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGc                                                            <  3:172197/149‑1 (MQ=255)
                                                                                        ccTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATgg                                                       <  5:307848/148‑1 (MQ=255)
                                                                                          tGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTg                                                    >  8:150354/1‑149 (MQ=255)
                                                                                                aCTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGcagcag                                              >  6:230948/1‑149 (MQ=255)
                                                                                                    gcgcATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGcagcag                                              >  2:2586/1‑145 (MQ=255)
                                                                                                    gcgcATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGcagcag                                              <  1:2586/145‑1 (MQ=255)
                                                                                                    gcgcATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGcagcag                                              >  2:8413/1‑145 (MQ=255)
                                                                                                    gcgcATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGcagcag                                              <  1:8413/145‑1 (MQ=255)
                                                                                                    gcgcATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGaa                                          >  4:171243/1‑149 (MQ=255)
                                                                                                           aGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGACCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGgtgt                                    <  3:142265/148‑1 (MQ=255)
                                                                                                             tCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGAc                                 >  3:58060/1‑149 (MQ=255)
                                                                                                                 aCGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGcc                              <  5:230948/148‑1 (MQ=255)
                                                                                                                   gCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTgg                           >  1:46843/1‑149 (MQ=255)
                                                                                                                        aaCGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTcgac                      >  7:269617/1‑149 (MQ=255)
                                                                                                                         aCGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTc                         >  1:105975/1‑145 (MQ=255)
                                                                                                                         aCGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTc                         <  2:105975/145‑1 (MQ=255)
                                                                                                                            ccTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTcgaccgac                  <  3:171243/149‑1 (MQ=255)
                                                                                                                              ttCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGaa                                          <  6:305460/123‑1 (MQ=255)
                                                                                                                              ttCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGaa                                          >  5:305460/1‑123 (MQ=255)
                                                                                                                                cgcgAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTgca                                                  <  2:305220/113‑1 (MQ=255)
                                                                                                                                 gcgAAGGGCGCCCCATTGCCCTGATCGCCGATTAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCATGGTCGACCGACTTGcc             >  5:285536/1‑148 (MQ=255)
                                                                                                                                     aGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGcc         <  8:269617/149‑1 (MQ=255)
                                                                                                                                           ccccATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCa                                                           >  6:167990/1‑93 (MQ=255)
                                                                                                                                           ccccATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCa                                                           <  5:167990/93‑1 (MQ=255)
                                                                                                                                           ccccATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGaa   <  7:150354/149‑1 (MQ=255)
                                                                                                                                             ccATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAAt  >  1:14228/1‑148 (MQ=255)
                                                                                                                                              cATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGaa   >  2:239937/1‑146 (MQ=255)
                                                                                                                                              cATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGaa   <  1:239937/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
TGGCGGAACTGGCGAATCTTTCGGGCATGTCGCGGGCAAGCTTCGCTTCGCACTTTCGCGAGATCGTGGGTATTACGCCTGCGGACTACCTGGCGAACTGGCGCATCAGTCTGACGCAAAAACGCCTTCGCGAAGGGCGCCCCATTGCCCTGATCGCCGATGAGGTCGGCTACGAAAGTCCTTCTGCGCTAGCCCGAACCTTTCGCCGCAAGGTGGGTGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAAT  >  NC_002947/3364028‑3364316

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