Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,156,141 A→T 13.7% L217* (TTA→TAA)  PP_3659 ← transcriptional regulator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,156,1410AT13.7% 144.8 / 12.6 64L217* (TTA→TAA) PP_3659transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base T (5/4);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGAC  >  NC_002947/4155993‑4156288
                                                                                                                                                    |                                                                                                                                                   
ggTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTa                                                                                                                                                     <  5:284042/149‑1 (MQ=255)
    gAACTGGACCTGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCACTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                 <  6:219025/149‑1 (MQ=255)
    gAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGGGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                 <  6:223882/149‑1 (MQ=255)
       cTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGc                                                                                                                                              <  8:217014/149‑1 (MQ=255)
            cATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTTAAAAGGCacgcc                                                                                                                                         <  6:320597/149‑5 (MQ=255)
                cTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAagag                                                                                                                                      <  5:20371/148‑1 (MQ=255)
                          gggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGtt                                                                                                                           <  2:144307/149‑1 (MQ=255)
                           ggCGAGCTCATAGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                          <  7:246947/149‑1 (MQ=255)
                            gCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTg                                                                                                                         <  5:307388/149‑1 (MQ=255)
                              gAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACAGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGc                                                                                                                       <  2:284805/149‑1 (MQ=255)
                                gCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCa                                                                                                                      >  7:87480/1‑148 (MQ=255)
                                  tCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCTTCTTTACGGTATTTCACAAGGAAGTTAATTGGCTTGCATTTTAAAAAGGAAGGAAGAGAAACGGGCGTTTGGaaaa                                                                                                                    >  3:211068/1‑148 (MQ=255)
                                  tCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAAc                                                                                                                   >  1:216956/1‑149 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCt                                                                                                                 >  2:75640/1‑149 (MQ=255)
                                     tCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTAGTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTGTGGCAAACCtt                                                                                                                <  1:140240/149‑1 (MQ=255)
                                             aGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGc                                                                                                          >  8:207049/1‑147 (MQ=255)
                                               gCCACGTGTGTGGAGAGAGCACTTTTAATCGCCGCCGGAACATAATGCACGTACTGTCGTCGAAACGGTATTACCCAAGGAACTTAATTGACGTGCCTTTTAAAAAGGAAGGAAGACAAACGGGCGTTTGGCAAACATTGAACGCCAg                                                                                                       >  6:343610/1‑148 (MQ=255)
                                                    gtgtgtGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACTGGCGTTTGGCAAACCTTGTACGCCAGAAAGct                                                                                                 <  5:151684/149‑1 (MQ=255)
                                                         tgGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCAGGGTTTAGGTATTTGAAAATGAAGTTAATTTGCGTGGCTTTTTAAAAGGCAGGAAGAGAAAAGGGCGTTTTGCAAAACGTGTAAGCCAGAAAGCtgttct                                                                                            >  1:249894/1‑149 (MQ=255)
                                                           gagagaGCAGTTTTAATCGCCGCCGTATAATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGAAGTACGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTtctc                                                                                           <  1:75640/148‑1 (MQ=255)
                                                                agCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                             >  6:161771/1‑93 (MQ=255)
                                                                agCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                             <  5:161771/93‑1 (MQ=255)
                                                                agCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAaca                                                                                      >  3:68533/1‑148 (MQ=255)
                                                                        tAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAATAAGGCAGGAAGAGAAACGGGAGTTTGTCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACAGAGCGGAg                                                                              <  2:249894/148‑1 (MQ=255)
                                                                        tAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAg                                                                              <  2:216956/148‑1 (MQ=255)
                                                                          aTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGtt                                                                           <  4:68533/149‑1 (MQ=255)
                                                                              ccgccgGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTAt                                                                       <  8:87480/149‑1 (MQ=255)
                                                                                    gATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGccc                                                                 >  2:228446/1‑149 (MQ=255)
                                                                                     aTCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCt                                                                <  7:207049/149‑1 (MQ=255)
                                                                                     aTCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCt                                                                >  8:80939/1‑149 (MQ=255)
                                                                                      tCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTc                                                               >  3:289247/1‑149 (MQ=255)
                                                                                        aTAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCa                                                             >  5:105938/1‑149 (MQ=255)
                                                                                          aaTGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACt                                                           >  3:192300/1‑149 (MQ=255)
                                                                                          aaTGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACt                                                           >  5:224440/1‑149 (MQ=255)
                                                                                           aTGGTCGTACTGTCGGCGTTACGGTATTTCACAAGGAAGTTACTTGGCGGTCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTGTCGAACACAGCGGAGGTTCTATCGGCCCCCCACTc                                                          <  5:343610/149‑1 (MQ=255)
                                                                                               tCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAtt                                                      >  8:72570/1‑149 (MQ=255)
                                                                                                  tACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAACAGGCAGGAAGAGAAACGGGCGTGTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTcc                                                   <  2:64782/149‑1 (MQ=255)
                                                                                                  tACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTcc                                                   <  1:228446/149‑1 (MQ=255)
                                                                                                       tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTcac                                               >  3:266001/1‑148 (MQ=255)
                                                                                                       tcgtcgTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTcac                                               >  5:278934/1‑148 (MQ=255)
                                                                                                          tcgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACaaa                                            >  8:117093/1‑148 (MQ=255)
                                                                                                            gtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAgc                                         >  2:101015/1‑149 (MQ=255)
                                                                                                                     tttCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGa                                 >  1:202569/1‑148 (MQ=255)
                                                                                                                     tttCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGAt                                >  1:197610/1‑149 (MQ=255)
                                                                                                                        cacaAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCgg                                                                    >  2:225912/1‑110 (MQ=255)
                                                                                                                        cacaAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCgg                                                                    <  1:225912/110‑1 (MQ=255)
                                                                                                                          caAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTTGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCg                           >  5:7982/1‑149 (MQ=255)
                                                                                                                             ggAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTg                        <  6:105938/149‑1 (MQ=255)
                                                                                                                             ggAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTAAGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTATATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGATg                        >  1:124090/1‑149 (MQ=255)
                                                                                                                                 gTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtcgt                    >  1:265605/1‑149 (MQ=255)
                                                                                                                                      ttGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCAcc               >  1:226890/1‑149 (MQ=255)
                                                                                                                                        ttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtcgt                    >  1:17965/5‑142 (MQ=255)
                                                                                                                                        ttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtcgt                    <  2:17965/138‑1 (MQ=255)
                                                                                                                                        ggCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCg              <  6:224440/148‑1 (MQ=255)
                                                                                                                                         tcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGc            >  5:218703/4‑149 (MQ=255)
                                                                                                                                         gCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGc            <  7:72570/149‑1 (MQ=255)
                                                                                                                                            tGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCt           >  5:321260/1‑147 (MQ=255)
                                                                                                                                            tGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga         <  2:197610/149‑1 (MQ=255)
                                                                                                                                            tGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga         <  2:202569/149‑1 (MQ=255)
                                                                                                                                             gCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATa                               >  3:212342/1‑126 (MQ=255)
                                                                                                                                             gCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATa                               <  4:212342/126‑1 (MQ=255)
                                                                                                                                             gCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgag        >  6:341513/1‑149 (MQ=255)
                                                                                                                                               cTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgagga      <  6:321260/149‑1 (MQ=255)
                                                                                                                                               cTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgag        >  8:282316/1‑147 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTAGATAGTCGCTGTCGTCCACCGGCTgaggag     >  5:186521/1‑149 (MQ=255)
                                                                                                                                                  ttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg    <  1:101015/148‑1 (MQ=255)
                                                                                                                                                   tAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAggagga   <  7:117093/148‑1 (MQ=255)
                                                                                                                                                    aaaaaGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGAc  <  4:289247/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGAC  >  NC_002947/4155993‑4156288

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