Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,969 A→G 8.0% A303A (GCT→GCC) ‡ PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,9690AG 8.0% 170.6 / 5.2 63A303A (GCT→GCC) ‡PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base G (3/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.00e-01

CTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTC  >  NC_002947/3865821‑3866112
                                                                                                                                                    |                                                                                                                                                
cTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCa                                                                                                                                                  >  6:87854/1‑149 (MQ=255)
      gAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGccgcccgc                                                                                                                                            >  7:219264/1‑149 (MQ=255)
           gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                      >  1:224074/1‑149 (MQ=255)
           gcgcgAGAAGGCCTTGGCGGGCACTCACAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                      >  1:9123/1‑149 (MQ=255)
            cgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAACGCCCAGATCCTGCTGCGCATTCGGCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtggg                                                                                                                                     <  2:344876/149‑3 (MQ=255)
              cgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCg                                                                                                                                   <  6:140427/149‑1 (MQ=255)
               gagaAGGCCTTGGCGGGCACTCTCAGGTTCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGa                                                                                                                                  <  1:304858/149‑1 (MQ=255)
               gagaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGtgtggcg                                                                                                                                  <  2:9123/149‑4 (MQ=255)
               gagaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAAATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCACACCGCGGGCTGCGgtggggcg                                                                                                                                  <  7:183593/149‑6 (MQ=255)
                 gaAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAAt                                                                                                                                >  7:222159/1‑149 (MQ=255)
                  aaGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATg                                                                                                                               <  2:345036/149‑1 (MQ=255)
                    ggCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGa                                                                                                                             <  7:102458/149‑1 (MQ=255)
                     gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                                                                                                                            <  7:314665/149‑1 (MQ=255)
                     gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                                                                                                                            <  2:165311/149‑1 (MQ=255)
                     gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                                                                                                                            <  5:293034/149‑1 (MQ=255)
                             gggCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCg                                                                                                                    >  6:148108/1‑149 (MQ=255)
                                   tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              >  8:317968/1‑149 (MQ=255)
                                   tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              >  5:65349/1‑149 (MQ=255)
                                   tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              >  4:294815/1‑149 (MQ=255)
                                   tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGCGGTGCGAATGGAGGTACGCCGTGGCa                                                                                                               >  8:181296/1‑148 (MQ=255)
                                     tcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              <  4:249053/147‑1 (MQ=255)
                                             tctctctcAACATCAGATCAAAGCCAAGATCCTTCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCAAGCAGCGCCAGTGCCGCGCCCCGCCCCACCGCACCCAGCCGGCTGCGAATGGACGACAGACGTGGCACCATAAGCGcc                                                                                                    >  8:315683/1‑149 (MQ=255)
                                               tcGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCCGCCCGCGGGGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGcc                                                                                                    >  7:297178/1‑147 (MQ=255)
                                               tcGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACACGCTTCCCATCCATCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGcc                                                                                                    >  7:291904/1‑147 (MQ=255)
                                                 gcgcACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGa                                                                                                  <  8:275220/147‑1 (MQ=255)
                                                     aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            <  3:294815/149‑1 (MQ=255)
                                                     aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            <  6:146301/149‑1 (MQ=255)
                                                       cATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGga                                                                                          <  5:25424/149‑1 (MQ=255)
                                                         tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            >  3:203690/1‑145 (MQ=255)
                                                         tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            >  5:148687/1‑145 (MQ=255)
                                                             cTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGt                                                                                    >  2:297785/1‑149 (MQ=255)
                                                                aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                   <  3:349444/147‑1 (MQ=255)
                                                                aaaGCCCAGATCCCGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATTGtt                                                                                   <  5:355634/147‑1 (MQ=255)
                                                                    cccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAAcc                                                                             >  6:258235/1‑149 (MQ=255)
                                                                     ccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAAcc                                                                             >  3:322503/1‑148 (MQ=255)
                                                                     ccAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAAcc                                                                             >  7:305081/1‑148 (MQ=255)
                                                                       aGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGAc                                                                          >  3:21120/1‑149 (MQ=255)
                                                                              gctgcGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGc                                                                    >  5:11969/1‑148 (MQ=255)
                                                                                     aTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGc                                                            <  6:148687/149‑1 (MQ=255)
                                                                                        cGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGt                                                          <  5:148108/148‑1 (MQ=255)
                                                                                                gCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTcc                                                 <  4:203690/149‑1 (MQ=255)
                                                                                                    cccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAccc                                             >  3:217855/1‑149 (MQ=255)
                                                                                                    cccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAc                                               <  2:224074/147‑1 (MQ=255)
                                                                                                    cccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAc                                               <  7:181296/147‑1 (MQ=255)
                                                                                                     ccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                            <  4:217855/149‑1 (MQ=255)
                                                                                                       aTCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGTACCTCCACCCCtt                                          <  7:315683/149‑1 (MQ=255)
                                                                                                       aTCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCtt                                          >  4:228229/1‑149 (MQ=255)
                                                                                                          cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGa                                        <  6:11969/148‑1 (MQ=255)
                                                                                                           aGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACg                                      >  7:169639/1‑149 (MQ=255)
                                                                                                            gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACg                                      >  8:5268/1‑148 (MQ=255)
                                                                                                              aaCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCaa                                   <  5:87854/149‑1 (MQ=255)
                                                                                                              aaCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCaa                                   <  8:169639/149‑1 (MQ=255)
                                                                                                               aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                  <  4:21120/149‑1 (MQ=255)
                                                                                                                        gcgcCTTTTCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATTGAGGTCAGCCGTGGCACCATATGCGCAGAGGCCGGGAGATCGTTGAAACCTACCATCGCCACCTGTCGTGGCACCTCCACCCCTTGACGCAAGGCCTgcagc                         <  7:270533/149‑1 (MQ=255)
                                                                                                                         cgcCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgcagca                        >  6:389189/1‑149 (MQ=255)
                                                                                                                           ccTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCAcc                      <  7:317968/149‑1 (MQ=255)
                                                                                                                                 ccGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGcc                 <  4:322503/148‑1 (MQ=255)
                                                                                                                                    cgccccgcccCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGa             <  8:305081/149‑1 (MQ=255)
                                                                                                                                          cgcacACCGCGGGCTGCGGTGGGGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAg        >  2:366253/5‑149 (MQ=255)
                                                                                                                                            cccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtc     <  8:222159/149‑1 (MQ=255)
                                                                                                                                            cacACCGCGGGCTGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtc     >  6:161852/3‑149 (MQ=255)
                                                                                                                                            cacACCGCGGGCTGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtc     >  5:8736/3‑149 (MQ=255)
                                                                                                                                             ccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCTGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtcg    >  2:324334/1‑149 (MQ=255)
                                                                                                                                              cACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtcgt   >  5:220286/1‑149 (MQ=255)
                                                                                                                                               aCCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtcgtc  <  5:389189/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTC  >  NC_002947/3865821‑3866112

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: