Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,663,119 A→T 9.6% intergenic (+23/+15) purT → / ← yhjE phosphoribosylglycinamide formyltransferase/metabolite transport protein YhjE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,663,1190AT 9.6% 201.3 / 12.6 73intergenic (+23/+15)purT/yhjEphosphoribosylglycinamide formyltransferase/metabolite transport protein YhjE
Reads supporting (aligned to +/- strand):  ref base A (36/30);  new base T (4/3);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

GGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGTGC  >  NC_002947/1662977‑1663265
                                                                                                                                              |                                                                                                                                                  
ggCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGccc                                                                                                                                             >  1:95732/1‑150 (MQ=255)
ggCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGccc                                                                                                                                             <  2:1107630/150‑1 (MQ=255)
  cAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCa                                                                                                                                           >  2:289387/1‑150 (MQ=255)
    aGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCCAAGCGGCCCCAGt                                                                                                                                         <  2:1051875/150‑1 (MQ=255)
    aGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGt                                                                                                                                         <  1:289387/150‑1 (MQ=255)
      ccGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGAGGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGCGTGCTGGGGGCGGAGAGCGGAGCCAGTgg                                                                                                                                       >  2:913723/1‑150 (MQ=255)
       cGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGCCCCAg                                                                                                                                          >  2:684556/1‑146 (MQ=255)
       cGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGCCCCAg                                                                                                                                          <  1:684556/146‑1 (MQ=255)
       cGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGc                                                                                                                                      >  2:1343672/1‑150 (MQ=255)
               actGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGCCGAACTGGCGCGGACCACGGCGGCCCGCGCTTCGCAGGCGGGCAAGGGCTAGGTCTGATCCGAGAGTGCGGGGGCCTCAAACCGGCCCCAGTGGCTTACAGCg                                                                                                                              >  2:329213/3‑150 (MQ=255)
               aaTGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCg                                                                                                                              >  2:822178/1‑150 (MQ=255)
                  ggCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTgg                                                                                                                           <  1:912314/150‑1 (MQ=255)
                   gCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCCGTGGCTTACAGCGTGGt                                                                                                                          >  2:805864/1‑150 (MQ=255)
                     aCCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGCCGAGTTCTGATCCGCGAGAGCTCGGCCCGCAAAGCGGCCACAGTGGCTTACAGCGCGGtct                                                                                                                        >  2:1151524/1‑148 (MQ=255)
                       ccgccgGATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGGCAAGGTCGAGTTCTGCTCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTg                                                                                                                      <  1:817591/150‑1 (MQ=255)
                       ccgccgCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCAGCAAAGCGGCCCCAGTGGGTTACAGCGTGGTATTg                                                                                                                      >  1:330184/1‑150 (MQ=255)
                          ccgcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGc                                                                                                                                               <  2:1538779/122‑1 (MQ=255)
                          ccgcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGc                                                                                                                                               >  1:1538779/1‑122 (MQ=255)
                            gcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGt                                                                                                                 <  1:113063/150‑1 (MQ=255)
                                  gtgtgtGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGc                                                                                                                                               <  2:599249/114‑1 (MQ=255)
                                  gtgtgtGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGc                                                                                                                                               >  1:599249/1‑114 (MQ=255)
                                   tgtgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTc                                                                                                          <  1:946959/150‑1 (MQ=255)
                                     tgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAg                                                                                                        <  1:39908/150‑1 (MQ=255)
                                       tgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGc                                                                                                      >  1:1327702/1‑150 (MQ=255)
                                                 gTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGc                                                                                            >  2:726083/1‑150 (MQ=255)
                                                    aCGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGCGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAAt                                                                                         >  1:1341231/1‑150 (MQ=255)
                                                      gactCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCa                                                                                       <  2:497301/147‑1 (MQ=255)
                                                       aaTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAg                                                                                      <  2:859213/150‑1 (MQ=255)
                                                        aTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                     <  1:726083/150‑1 (MQ=255)
                                                            gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTa                                                                                                                         <  2:1189362/110‑1 (MQ=255)
                                                            gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTa                                                                                                                         >  1:1189362/1‑110 (MQ=255)
                                                            gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCa                                                                                                         >  2:165011/1‑126 (MQ=255)
                                                            gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCa                                                                                                         <  1:165011/126‑1 (MQ=255)
                                                                aaCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAc                                                                             <  1:822178/150‑1 (MQ=255)
                                                                 aCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAcc                                                                            >  1:403353/1‑150 (MQ=255)
                                                                         gggCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCAc                                                                    >  2:1083810/1‑150 (MQ=255)
                                                                             cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCg                                                                                                                                                 >  2:78714/1‑69 (MQ=255)
                                                                             cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCg                                                                                                                                                 <  1:78714/69‑1 (MQ=255)
                                                                             cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAg                                                                                      >  1:152652/1‑128 (MQ=255)
                                                                             cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAg                                                                                      <  2:152652/128‑1 (MQ=255)
                                                                             cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGac                                                                >  2:802895/1‑150 (MQ=255)
                                                                             cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGac                                                                >  1:1146795/1‑150 (MQ=255)
                                                                              aaGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGaca                                                               >  1:383083/1‑150 (MQ=255)
                                                                                   gACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATaa                                                          >  1:1466490/1‑150 (MQ=255)
                                                                                      ccGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAAccc                                                       >  2:206372/1‑150 (MQ=255)
                                                                                       cgcgCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAAcccc                                                      >  2:1225419/1‑150 (MQ=255)
                                                                                       cgcgCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAAcccc                                                      <  1:802895/150‑1 (MQ=255)
                                                                                        gcgcTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACAGCAATCAAGCTGATCACCGCCGCCATCGCCAACTAAccccc                                                     >  1:1411838/1‑150 (MQ=255)
                                                                                                cAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGc                                             >  2:1300773/1‑150 (MQ=255)
                                                                                                cAGGAGGTTAAAGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGc                                             <  2:383083/150‑1 (MQ=255)
                                                                                                          aggcggcgTTCTTACCAGATAGTGCTGTGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAAccccc                                   <  2:1327702/143‑1 (MQ=255)
                                                                                                          aGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAAccccc                                   >  1:734659/1‑150 (MQ=255)
                                                                                                           ggTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAAcccccc                                  <  2:95732/150‑1 (MQ=255)
                                                                                                              cGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATg                               <  1:206372/150‑1 (MQ=255)
                                                                                                              cGAGTTCTGATCCGAGAGGGCTGGGTCCGCAAAGCGGCCCCACTCGCTTACACCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATNANCNCCGCCACCGACACANNACCCCCCACCCAGCTCAAACCCCCCATg                               <  1:913723/150‑1 (MQ=255)
                                                                                                                  ttCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTATAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCa                           <  1:1343672/150‑1 (MQ=255)
                                                                                                                  ttCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCa                           >  1:1151132/1‑150 (MQ=255)
                                                                                                                   tCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCAc                          <  2:1341231/150‑1 (MQ=255)
                                                                                                                          cGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACg                                                                                              <  2:1365564/75‑1 (MQ=255)
                                                                                                                          cGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCGTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACg                                                                                              >  1:1365564/1‑75 (MQ=255)
                                                                                                                           gagaGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCACCTTCt                  >  2:934277/1‑150 (MQ=255)
                                                                                                                               ccacTGGGGCCGCTTTGCGGCCCCAGCGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCAc                          >  1:1308314/4‑138 (MQ=255)
                                                                                                                               ccacTGGGGCCGCTTTGCGGCCCCAGCGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCAc                          <  2:1308314/135‑1 (MQ=255)
                                                                                                                                tGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACGCGTAAGCCCCCACCCAGCTCAAACCCCGCATGCTAACCAGCTCCTGAGCa             >  1:596326/1‑150 (MQ=255)
                                                                                                                                cacTGGGGCCGCTTTGCGGCCCCAGCGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCa             >  1:1017100/3‑150 (MQ=255)
                                                                                                                                 acTGGGGCCGCTTTGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGc                                             >  2:157537/2‑117 (MQ=255)
                                                                                                                                 acTGGGGCCGCTTTGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGc                                             <  1:157537/116‑1 (MQ=255)
                                                                                                                                   tGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCCCGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTGAAGCCCCCCACGCTCACCAGCGTCCACGCAAAg          >  2:524511/1‑150 (MQ=255)
                                                                                                                                    ggggCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCttt                                                                                                           >  1:1554962/1‑52 (MQ=255)
                                                                                                                                    ggggCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCttt                                                                                                           <  2:1554962/52‑1 (MQ=255)
                                                                                                                                       gCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACg      <  2:1146795/150‑1 (MQ=255)
                                                                                                                                         cGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGt    <  2:1151132/150‑1 (MQ=255)
                                                                                                                                           cAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCATGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGtgc  >  2:1294851/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGTGC  >  NC_002947/1662977‑1663265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑