Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,319,052 C→T 5.9% A8V (GCG→GTG)  PP_3790 → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,319,0520CT 5.9% 202.7 / 2.1 68A8V (GCG→GTG) PP_3790diaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base C (32/32);  new base T (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

GCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACGGTTGTCG  >  NC_002947/4318915‑4319197
                                                                                                                                         |                                                                                                                                                 
gcgcCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCgg                                                                                                                                                 >  1:359664/1‑140 (MQ=255)
gcgcCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCgg                                                                                                                                                 <  2:359664/140‑1 (MQ=255)
  gcCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGAc                                                                                                                                      >  8:231227/1‑149 (MQ=255)
   cccTCGACTGACCCCGGTCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACt                                                                                                                                     <  7:305680/149‑1 (MQ=255)
          cTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGt                                                                                                                              <  3:185472/149‑1 (MQ=255)
           tGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGta                                                                                                                             >  2:246827/1‑149 (MQ=255)
             aCCCCGGGCAACTTTGCGCGTGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGtat                                                                                                                            <  1:240171/148‑1 (MQ=255)
               cccGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATAcc                                                                                                                         <  1:137898/149‑1 (MQ=255)
                 cGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGta                                                                                                                             <  4:97173/143‑1 (MQ=255)
                 cGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGta                                                                                                                             >  3:97173/1‑143 (MQ=255)
                  gggCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCgg                                                                                                                       >  2:54780/1‑148 (MQ=255)
                        cTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGTGGCTGGAAATGACTACCTTGTATACCGGGACTGCg                                                                                                                >  3:168270/1‑149 (MQ=255)
                           tGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGTGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTa                                                                                                              <  4:168270/148‑1 (MQ=255)
                                  ggAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGa                                                                                                      <  4:123313/149‑1 (MQ=255)
                                      aTCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGc                                                                                                  <  8:52587/149‑1 (MQ=255)
                                         tttAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTCATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTc                                                                                                >  5:283360/1‑148 (MQ=255)
                                          ttAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTcc                                                                                               >  6:87143/1‑148 (MQ=255)
                                           tAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGt                                                                                                                              <  3:122361/116‑1 (MQ=255)
                                           tAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGt                                                                                                                              >  4:122361/1‑116 (MQ=255)
                                              tACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAAcc                                                                                          >  6:38476/1‑149 (MQ=255)
                                               aCTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAAcc                                                                                          >  4:158460/1‑148 (MQ=255)
                                                      tACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGAttt                                                                                  >  8:272814/1‑149 (MQ=255)
                                                        cgcgCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGGAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTc                                                                                 <  6:284155/148‑1 (MQ=255)
                                                         tcgcATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTAGTATACCGTCACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCAc                                                                               >  6:181317/2‑149 (MQ=255)
                                                         gcgcATTGGTGGTCAATGTCTCACCTTGACTTCTCTGATCATTCTGGTGCCTTAGGTAATGCAGTTTTATAAATACCTTGCGGCTGGAAATGACTACCTTGTATTCCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCAc                                                                               <  7:268962/149‑1 (MQ=255)
                                                         gcgcATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCAc                                                                               >  6:82834/1‑149 (MQ=255)
                                                           gcATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCg                                                                                                                        >  7:194447/1‑106 (MQ=255)
                                                           gcATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCg                                                                                                                        <  8:194447/106‑1 (MQ=255)
                                                           gcATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACg                                                                              >  2:51823/1‑148 (MQ=255)
                                                             aTTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTAt                                                                           >  5:344873/1‑149 (MQ=255)
                                                                 gtcgtcAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGc                                                                       <  1:246827/149‑1 (MQ=255)
                                                                 gtcgtcAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGc                                                                       >  5:353445/1‑149 (MQ=255)
                                                                  tcgtcAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGc                                                                       <  1:51823/148‑1 (MQ=255)
                                                                    gtcAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGcgac                                                                    >  5:191825/1‑149 (MQ=255)
                                                                       aaTGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGcgaccga                                                                 >  5:102955/1‑149 (MQ=255)
                                                                       aaTGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGcgaccga                                                                 <  7:272814/149‑1 (MQ=255)
                                                                        aTGTCTCACCCTCACTTCTCTGATCATTCTCGTGCCTTAGCTAATGAAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCTGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGcgaacga                                                                 >  1:244536/1‑148 (MQ=255)
                                                                           tctcACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACAcc                                                             >  7:187378/1‑149 (MQ=255)
                                                                             tcACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACcgc                                                           >  6:323041/1‑149 (MQ=255)
                                                                               aCCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCTGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCgg                                                         >  8:127194/1‑149 (MQ=255)
                                                                               aCCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGTGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGcgaccga                                                                 <  8:108127/141‑1 (MQ=255)
                                                                               aCCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGTGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGcgaccga                                                                 >  7:108127/1‑141 (MQ=255)
                                                                                cccTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTAGCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCggg                                                        <  2:75594/149‑1 (MQ=255)
                                                                                 ccTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGt                                                       >  6:371634/1‑149 (MQ=255)
                                                                                 ccTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGt                                                       >  1:206954/1‑149 (MQ=255)
                                                                                     aCTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTggg                                                   >  6:35413/1‑149 (MQ=255)
                                                                                            tGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTAt                                                                           <  1:35166/118‑1 (MQ=255)
                                                                                            tGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTAt                                                                           >  2:35166/1‑118 (MQ=255)
                                                                                              aTCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGa                                             <  5:38476/146‑1 (MQ=255)
                                                                                               tCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAATCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGc                                         <  1:54780/149‑1 (MQ=255)
                                                                                               tCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGc                                         <  5:87143/149‑1 (MQ=255)
                                                                                                 aTTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCa                                        <  5:371634/148‑1 (MQ=255)
                                                                                                 aTTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCa                                        <  6:353445/148‑1 (MQ=255)
                                                                                                       gTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGt                                 >  3:37374/1‑149 (MQ=255)
                                                                                                            ttAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCgg                            >  1:120956/1‑149 (MQ=255)
                                                                                                            ttAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCgg                            >  1:43663/1‑149 (MQ=255)
                                                                                                               gCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTg                          <  7:231227/148‑1 (MQ=255)
                                                                                                                      cAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAaaga                  <  4:37374/149‑1 (MQ=255)
                                                                                                                      cAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAaaga                  <  2:43663/149‑1 (MQ=255)
                                                                                                                        gTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATTCCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAaagaag                <  5:323041/149‑1 (MQ=255)
                                                                                                                        gTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGAATGAGTACATTTTGATTGCTCCAAACCGATGATTTCACGTATTTGAGACCGACAAAGCGGGTTGGGCAGCGACGTAATAATGGTGACTGTGATCaaaaaat                >  5:80609/1‑148 (MQ=255)
                                                                                                                           ttATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGATTGATTGCTCAAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTTGGCAGCGACGGCATACTGGTGACGGTGATCAAAGAAGGt              >  5:267173/1‑148 (MQ=255)
                                                                                                                              taAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAAc          <  5:181317/149‑1 (MQ=255)
                                                                                                                                aaTACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACgg        <  6:344873/149‑1 (MQ=255)
                                                                                                                                aaTACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACgg        <  6:191825/149‑1 (MQ=255)
                                                                                                                                  tACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACGGtt      <  2:120956/149‑1 (MQ=255)
                                                                                                                                    ccATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACGGTTGt    <  3:158460/149‑1 (MQ=255)
                                                                                                                                      aTGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACGGTTGTCg  <  6:102955/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                 
GCGCCCTCGACTGACCCCGGGCAACTTTGCGCGGGGAGATCTTTAATACTTCAGTACGCGCATTGGTCGTCAATGTCTCACCCTCACTTCTCTGATCATTCTGGTGCCTTAGCTAATGCAGTTTTATAAATACCATGCGGCTGGAAATGACTACCTTGTATACCGGGACTGCGTACCGTTTGATTGCTCCAAACCGATGATTTCACGTATTTGCGACCGACACCGCGGGTTGGGCAGCGACGGCATACTGGTGCCGGTGATCAAAGAAGGTGAACGGTTGTCG  >  NC_002947/4318915‑4319197

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: