Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 772,182 C→A 7.3% intergenic (+43/‑156) PP_0663 → / → PP_0664 AsnC family transcriptional regulator/homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947772,1820CA 7.3% 200.5 / 4.1 63intergenic (+43/‑156)PP_0663/PP_0664AsnC family transcriptional regulator/homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base A (3/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAGGACA  >  NC_002947/772039‑772330
                                                                                                                                               |                                                                                                                                                    
ccTCCAGCCCGGAGCGGGTGCAGCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGGGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGGCGGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCa                                                                                                                                                 <  7:266508/149‑1 (MQ=255)
 cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAg                                                                                                                                                >  1:125664/1‑149 (MQ=255)
 cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAg                                                                                                                                                >  7:142133/1‑149 (MQ=255)
 cTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAg                                                                                                                                                >  3:89092/1‑149 (MQ=255)
  tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAgg                                                                                                                                               >  3:121209/1‑149 (MQ=255)
  tCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTCGCCGCGGTGAACTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAgg                                                                                                                                               >  1:119064/1‑149 (MQ=255)
    cAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAggcg                                                                                                                                             >  2:90611/1‑149 (MQ=255)
      gCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAggcggc                                                                                                                                           <  1:90611/149‑1 (MQ=255)
      gCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGGCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAggcggc                                                                                                                                           <  5:222701/149‑1 (MQ=255)
         cGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGtt                                                                                                                                        <  4:89092/149‑1 (MQ=255)
          ggAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTtg                                                                                                                                       >  1:136487/1‑149 (MQ=255)
                  tGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGtttt                                                                                                                               <  6:330405/149‑1 (MQ=255)
                        aTGCGTCGAAGGAACTCGCGGCCCTGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTg                                                                                                                         <  8:274774/149‑1 (MQ=255)
                         tgtgtgGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTgg                                                                                                                        <  8:304018/149‑1 (MQ=255)
                                 aggaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTg                                                                                                                >  3:295357/1‑149 (MQ=255)
                                   gaaATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgcg                                                                                                              >  5:329823/1‑149 (MQ=255)
                                    aaaTCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTgcgc                                                                                                             <  6:98996/149‑1 (MQ=255)
                                         gCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCttt                                                                                                        <  5:212976/149‑1 (MQ=255)
                                             ccATGCCAGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGgcggcg                                                                                                                                          <  6:265747/111‑1 (MQ=38)
                                             ccATGCCAGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGCGCGAACCCAGgcggcg                                                                                                                                          >  5:265747/1‑111 (MQ=38)
                                                    cGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTCAATTGGATAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTGGGGGGTTTTGCCATGGCGCAGGTGCGCGCTTTGCAATGGCCa                                                                                              >  8:173550/1‑148 (MQ=255)
                                                           gagaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGc                                                                                       >  1:332796/1‑148 (MQ=255)
                                                             gaACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCag                                                                                    <  4:295357/149‑2 (MQ=255)
                                                              cacacGGGGACTTCCGTCGGCCTGGCCGCGGGGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCagt                                                                                   <  7:173550/148‑3 (MQ=255)
                                                              aaCACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAATCTGGCCagt                                                                                   <  2:125664/149‑3 (MQ=255)
                                                                     tGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAAAAAATACa                                                                            <  2:298330/149‑1 (MQ=255)
                                                                         ttcgttcgTCCTGGCCGCGCTGGTCTGAACCGCGTCTTTCCATTGGATATCACAACGCCGCCTGGATTCGCGCCCACGCGGCGTTGTGCGTTTTTCCCTGGCGAAGGTGCGCGCTTTGCACTGGCCACATTGGCCAATAAATACAGCg                                                                         >  3:213809/1‑148 (MQ=255)
                                                                            gttcgtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGAtt                                                                     >  3:15805/1‑149 (MQ=255)
                                                                                gtCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGcc                                                                  <  3:15003/148‑1 (MQ=255)
                                                                                 tCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCa                                                                 >  8:116709/1‑148 (MQ=255)
                                                                                  ccTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCata                                                               >  7:27383/1‑149 (MQ=255)
                                                                                    tGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCatat                                                              >  6:107615/1‑148 (MQ=255)
                                                                                      gCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGcc                                                                                      <  2:336979/122‑1 (MQ=255)
                                                                                      gCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGcc                                                                                      >  1:336979/1‑122 (MQ=255)
                                                                                      gCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATAtt                                                             <  8:142133/147‑1 (MQ=255)
                                                                                       ccGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCaa                                                          <  6:329823/149‑1 (MQ=255)
                                                                                        cgcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAAc                                                         <  7:201137/149‑1 (MQ=255)
                                                                                        cgcgGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAAc                                                         >  8:201137/1‑149 (MQ=255)
                                                                                             tggtCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTgg                                                      >  3:239360/1‑147 (MQ=255)
                                                                                                     aCCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACgg                                            >  8:171188/1‑149 (MQ=255)
                                                                                                       cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACt                                                        >  5:244971/1‑135 (MQ=255)
                                                                                                       cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACt                                                        <  6:244971/135‑1 (MQ=255)
                                                                                                       cgcgTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGCATTGCCATATTCAACTGGCAGTGCACGGc                                           >  4:26559/1‑148 (MQ=255)
                                                                                                          gTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACt                                                        >  8:6388/1‑132 (MQ=255)
                                                                                                          gTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACt                                                        <  7:6388/132‑1 (MQ=255)
                                                                                                               tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTAc                                     <  4:15805/146‑1 (MQ=255)
                                                                                                               tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGc                                   <  2:136487/148‑1 (MQ=255)
                                                                                                               tGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGc                                   <  4:121209/148‑1 (MQ=255)
                                                                                                                  aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGct                                          >  8:40507/1‑138 (MQ=255)
                                                                                                                  aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGct                                          <  7:40507/138‑1 (MQ=255)
                                                                                                                  aTTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGt                               >  1:196910/1‑149 (MQ=255)
                                                                                                                      gCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCt                            >  2:181707/1‑148 (MQ=255)
                                                                                                                          gCACTACGCGGAGTGGGTTCGCGGCCAGGCGGCGTTGTGGGTTTTGCCCTGTGGCAGGTGCGCGCTTTGCTCTGGCAAGATTTTCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTTCACGGCTCTTAGGCTTGTCCTCATc                        <  4:213809/148‑1 (MQ=255)
                                                                                                                          gCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCTATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCAt                         <  2:119064/147‑1 (MQ=255)
                                                                                                                             caACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                    >  7:40223/1‑149 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                    >  5:10975/1‑148 (MQ=255)
                                                                                                                              aaCGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGAcc                    >  8:42208/1‑148 (MQ=255)
                                                                                                                                cgccgcCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCATACGCTAGTCCTCATCgccaca                  >  8:346957/1‑143 (MQ=255)
                                                                                                                                  ccgccTGGGCGCGAACCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGAGCACGTCTCTTACGCTAGTCCTCATCGACCCAcag               >  8:152246/1‑147 (MQ=255)
                                                                                                                                   cgccTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGc              <  8:27383/149‑1 (MQ=255)
                                                                                                                                      cTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTg           <  1:181707/149‑1 (MQ=255)
                                                                                                                                         ggTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTgga         <  7:171188/148‑1 (MQ=255)
                                                                                                                                           ttCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggacag      <  7:116709/149‑1 (MQ=255)
                                                                                                                                               cgcCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCTCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTggacaggaca  <  3:26559/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCTCCAGCCCGGAGCGGGTGCACCATGTGTGGAAGGAAATCGCTGCCATGCCCGGGGTGGAGAACACGGTGACTTCGTTCGTCCTGGCCGCGGTGGTCTGAACCGCGTCTTTGCATTGGCTAGCACAACGCCGCCTGGGTTCGCGCCCAGGCGGCGTTGTGCGTTTTGCCCTGGCGCAGGTGCGCGCTTTGCACTGGCCAGATTGGCCAATAAATACAGCGGATTGCCATATTCAACTGGCAGTGCACGGCTCTTACGCTAGTCCTCATCGACCCACTGCCTGGACAGGACA  >  NC_002947/772039‑772330

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