Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,352,175 G→C 6.4% intergenic (+65/+77) PP_2947 → / ← PP_2948 transcriptional regulator MvaT/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,352,1750GC 6.4% 251.9 / 5.8 78intergenic (+65/+77)PP_2947/PP_2948transcriptional regulator MvaT/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (39/34);  new base C (2/3);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

AGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGA  >  NC_002947/3352028‑3352322
                                                                                                                                                   |                                                                                                                                                   
agaCCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGc                                                                                                                                                    <  2:68579/149‑1 (MQ=255)
  aCCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGCCAGCGt                                                                                                                                                  >  4:134381/1‑149 (MQ=255)
    cAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCg                                                                                                                                                >  8:33260/1‑149 (MQ=255)
    cAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCg                                                                                                                                                >  1:241743/1‑149 (MQ=255)
          tGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGcttcct                                                                                                                                          <  8:319765/149‑1 (MQ=255)
          tGGCAACCACAAGCAGGTCAAGGAATGGAATGCCCAGTACGGGGCTGAAGAAGTGGAATCCTGGGGTTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTTTCTGAAACGACTTGTTCTTTTTCGAGGAGGCGACGCTGTCGcttcct                                                                                                                                          <  2:298823/149‑1 (MQ=255)
            gCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCctcctc                                                                                                                                         <  6:141004/148‑1 (MQ=255)
            gCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGAAGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCctcctc                                                                                                                                         >  6:234619/1‑148 (MQ=255)
              aaCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACCGCGTCGCCTCCTCGa                                                                                                                                       >  8:90030/1‑148 (MQ=255)
               aCCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGaaa                                                                                                                                     <  4:85155/149‑1 (MQ=255)
                      gctgctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCg                                                                                                                              <  2:1016/149‑1 (MQ=255)
                          ctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAgg                                                                                                                                 >  2:199208/1‑142 (MQ=255)
                          ctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAgg                                                                                                                                 <  1:199208/142‑1 (MQ=255)
                          ctCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCa                                                                                                                          >  7:47604/1‑149 (MQ=255)
                           tCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCAt                                                                                                                         >  7:230789/1‑149 (MQ=255)
                           tCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGGGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCAt                                                                                                                         >  3:315715/1‑149 (MQ=255)
                             aaGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcc                                                                                                                                              >  1:105781/1‑126 (MQ=255)
                             aaGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGcc                                                                                                                                              <  2:105781/126‑1 (MQ=255)
                               ggAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtg                                                                                                                     >  4:86820/1‑149 (MQ=255)
                                         gCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGa                                                                                                           <  7:33260/149‑1 (MQ=255)
                                          cccAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGaa                                                                                                          >  6:21554/1‑149 (MQ=255)
                                              gTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAg                                                                                                      <  8:230789/149‑1 (MQ=255)
                                                  ggggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGa                                                                                                  >  6:90738/1‑149 (MQ=255)
                                                  ggggCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGa                                                                                                  >  5:79277/1‑149 (MQ=255)
                                                         aagaagTGGAATCCTGGGTGTCCTGATTTGTGCTGAAGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGtt                                                                                           <  5:234619/149‑1 (MQ=255)
                                                             agTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGTCGGAATGAGCTGACGTGGTTAAAc                                                                                       <  3:62968/149‑1 (MQ=255)
                                                              gTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACa                                                                                      <  6:257682/149‑1 (MQ=255)
                                                                         ggTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTc                                                                             <  4:315715/147‑1 (MQ=255)
                                                                               cTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgggtgg                                                                                                                >  3:290027/1‑106 (MQ=255)
                                                                               cTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgtgggtgg                                                                                                                <  4:290027/106‑1 (MQ=255)
                                                                               cTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAg                                                                                     <  1:120213/133‑1 (MQ=255)
                                                                               cTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAg                                                                                     >  2:120213/1‑133 (MQ=255)
                                                                               cTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTa                                                                      <  5:21554/148‑1 (MQ=255)
                                                                                       tgCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAg                                                             >  5:324723/1‑149 (MQ=255)
                                                                                         cTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCg                                                                                                             >  3:30518/1‑99 (MQ=255)
                                                                                         cTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCg                                                                                                             <  4:30518/99‑1 (MQ=255)
                                                                                         cTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAgt                                                            <  3:86820/148‑1 (MQ=255)
                                                                                                 gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCTTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                    >  2:266889/1‑148 (MQ=255)
                                                                                                 gCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGttt                                                    <  8:11618/148‑1 (MQ=255)
                                                                                                  cccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgcg                                                                  >  2:265191/1‑133 (MQ=255)
                                                                                                  cccTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAgcg                                                                  <  1:265191/133‑1 (MQ=255)
                                                                                                    cTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCt                                                                                                    >  8:170802/1‑97 (MQ=255)
                                                                                                    cTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCt                                                                                                    <  7:170802/97‑1 (MQ=255)
                                                                                                       aTTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCttcctc                                                                                                                                         >  7:107179/1‑57 (MQ=18)
                                                                                                       aTTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAGGCGACGCTGTCGCttcctc                                                                                                                                         <  8:107179/57‑1 (MQ=18)
                                                                                                       aTTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATa                                               <  4:165415/147‑1 (MQ=255)
                                                                                                               aTCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAt                                     >  8:45170/1‑149 (MQ=255)
                                                                                                                tCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAtt                                    >  5:53747/1‑149 (MQ=255)
                                                                                                                  tGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAg                                                                                                                                  >  4:43772/1‑53 (MQ=255)
                                                                                                                  tGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAg                                                                                                                                  <  3:43772/53‑1 (MQ=255)
                                                                                                                   gAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              <  2:145224/120‑1 (MQ=255)
                                                                                                                   gAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              >  1:145224/1‑120 (MQ=255)
                                                                                                                   gAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                 >  8:320614/1‑149 (MQ=255)
                                                                                                                    aaaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAAtttt                                  <  6:67446/147‑1 (MQ=255)
                                                                                                                    aaaCGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTg                                 <  6:324723/148‑1 (MQ=255)
                                                                                                                         ccTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCAt                             >  7:79108/1‑147 (MQ=255)
                                                                                                                          cTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAg                          <  8:47604/149‑1 (MQ=255)
                                                                                                                             gTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGcc                        >  6:186976/1‑148 (MQ=255)
                                                                                                                             gTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGc                         >  7:150150/1‑147 (MQ=255)
                                                                                                                             gTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGc                         <  7:90030/147‑1 (MQ=255)
                                                                                                                             gTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGc                         <  8:150150/147‑1 (MQ=255)
                                                                                                                               cctgtttCGAGGAGGCGACGCTGTCGCTTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGAt                     >  6:109114/5‑149 (MQ=39)
                                                                                                                                cTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGAtt                    >  1:346941/1‑149 (MQ=255)
                                                                                                                                 tttttCGAGGAAGCGACAGCGTCGCCTCGTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTg                   >  1:214542/1‑149 (MQ=255)
                                                                                                                                 tttttCGAGGAAGCGACAGCGTCGCCTCCTCGAACCAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAGCAGGATGTTTCGGCACTAAGCGCTCAGTGTCCGTTTCTATAGTTCGTAAATTTTGGCATCAGCCGATTg                   >  7:274069/1‑149 (MQ=255)
                                                                                                                                 tttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgt                                                                                                                      <  3:244578/50‑1 (MQ=255)
                                                                                                                                 tttttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCgt                                                                                                                      >  4:244578/1‑50 (MQ=255)
                                                                                                                                   tttCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTgg                  <  5:90738/148‑1 (MQ=255)
                                                                                                                                     tCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGc                <  2:214542/148‑1 (MQ=255)
                                                                                                                                      cGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCg               <  2:241743/148‑1 (MQ=255)
                                                                                                                                        aGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              >  5:248338/1‑99 (MQ=255)
                                                                                                                                        aGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCa                                                              <  6:248338/99‑1 (MQ=255)
                                                                                                                                        aGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACa            >  1:272144/1‑149 (MQ=255)
                                                                                                                                         ggAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACAt           >  4:251821/1‑149 (MQ=255)
                                                                                                                                          gAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACAtt          >  5:97218/1‑149 (MQ=255)
                                                                                                                                              cGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTAGGCGACATTGaaa      >  1:293636/1‑149 (MQ=255)
                                                                                                                                                aCAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACg    >  1:200833/1‑149 (MQ=255)
                                                                                                                                                  aGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACgga  <  1:266889/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGACCAAGGGTGGCAACCACAAGCTGCTCAAGGAATGGAAAGCCCAGTACGGGGCAGAAGAAGTGGAATCCTGGGTGTCCTGATTTGTGCTGATGTCGCCCTCATTGCGCTATCTGAAACGCCTTGTTCTTTTTCGAGGAAGCGACAGCGTCGCCTCCTCGAAACAGGGCGTTCATCGTGGGTGGGCGGAATGAGCTGACGTGGTTAAACAGGAAGTTTCGCCACTAAGCGCTCAGTGTCCGTTTCAATACTTCGTAAATTTTGGCATCAGCCGATTGGGCGACATTGAAACGGA  >  NC_002947/3352028‑3352322

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: ‑