Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,971 C→G 7.9% A303P (GCT→CCT) ‡ 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,9710CG 7.9% 197.1 / 7.0 63A303P (GCT→CCT) ‡PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base G (3/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

GAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAG  >  NC_002947/3865827‑3866119
                                                                                                                                                |                                                                                                                                                     
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                                                                                                                                   <  5:293034/149‑1 (MQ=255)
               gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                                                                                                                                   <  7:314665/149‑1 (MQ=255)
               gCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAg                                                                                                                                   <  2:165311/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                                                                                                   <  6:146301/149‑1 (MQ=255)
                                               aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                   <  3:294815/149‑1 (MQ=255)
                                                 cATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGga                                                                                                 <  5:25424/149‑1 (MQ=255)
                                                   tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                   >  5:148687/1‑145 (MQ=255)
                                                   tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                                   >  3:203690/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                                          <  8:169639/149‑1 (MQ=255)
                                                                                                        aaCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCaa                                          <  5:87854/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)
                                                                                                                                                cccGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTa     >  5:155872/1‑146 (MQ=255)
                                                                                                                                                cccGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACag  <  1:297785/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTCGTTACAG  >  NC_002947/3865827‑3866119

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: