Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,105,450 A→T 9.8% intergenic (+211/‑76) hisG → / → hisD ATP phosphoribosyltransferase/histidinol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,105,4500AT 9.8% 612.6 / 8.9 222intergenic (+211/‑76)hisG/hisDATP phosphoribosyltransferase/histidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (126/68);  new base T (2/26);  total (128/94)
Fisher's exact test for biased strand distribution p-value = 3.33e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

CCCGCGCGGATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCT  >  NC_002947/1105303‑1105599
                                                                                                                                                    |                                                                                                                                                     
   gcgcgGATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGAtc                                                                                                                                                  <  2:113472/150‑1 (MQ=255)
     gcgGATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGAtcgt                                                                                                                                                <  2:2779946/150‑1 (MQ=255)
      cgGATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGAtcgtc                                                                                                                                               <  2:7994456/150‑1 (MQ=255)
       ggATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAGATGAGCATACCGAGCCTGTTTTTATCGTAGCTTGAtcgtcg                                                                                                                                              >  1:4869883/1‑150 (MQ=255)
       ggATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGAtcgtcg                                                                                                                                              >  1:8217089/1‑150 (MQ=255)
        gATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGACAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTcgc                                                                                                                                             <  2:8163099/150‑1 (MQ=255)
        gATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTcgc                                                                                                                                             <  2:8324306/150‑1 (MQ=255)
        gATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTcgc                                                                                                                                             <  1:7320351/150‑1 (MQ=255)
            gACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAg                                                                                                                                         <  2:5935891/150‑1 (MQ=255)
              cTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGc                                                                                                                                       >  1:8526576/1‑150 (MQ=255)
               tGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCa                                                                                                                                      <  1:6233127/150‑1 (MQ=255)
                 gTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACt                                                                                                                                    >  1:1179321/1‑150 (MQ=255)
                 gTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACt                                                                                                                                    >  2:5511112/1‑150 (MQ=255)
                 gTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACt                                                                                                                                    >  1:4718449/1‑150 (MQ=255)
                      ttGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCg                                                                                                                               >  2:9089383/1‑150 (MQ=255)
                        gggATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGca                                                                                                                             >  2:632513/1‑150 (MQ=255)
                        gggATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACGTTTCGca                                                                                                                             >  2:5209718/1‑150 (MQ=255)
                                  gAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTa                                                                                                                   >  1:4622602/1‑150 (MQ=255)
                                  gAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTa                                                                                                                   >  1:3568699/1‑150 (MQ=255)
                                   aaacGTCGCCGCGGGGCTATAATGAAACAACCAAAACGGGCGTGTACGGGCGGGGGGCGCGCCCGACTTTGCTGGTTGAAATTGAGAATACGGAGCCTTTTTTTAACGTAGATTGATCGTCGCGCAGCCACTTTTCGCACAATCCCCTAg                                                                                                                  <  2:2676792/150‑1 (MQ=255)
                                    atgatCGCCGCTCGGCGATCATGCAAGAACGAAAACAGGCGAGGACCGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAgg                                                                                                                 >  1:8059939/5‑150 (MQ=255)
                                          gCCGCGCGGCGTACATGAAAGACCGAAACCCGGCGTGGAACGGCCGGGGTCGCGCCCGACTTTGCGGGTTGCAAATGAGCATACCGAGCCTGTTTTTACCGTCGCTTGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCt                                                                                                           <  2:4283752/150‑1 (MQ=255)
                                            cGAGCGGCGCTCATGCAAGACCGAAAACCGGCGAGGAACGGCGGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATCCCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGAAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGa                                                                                                         <  2:6630884/150‑1 (MQ=255)
                                             gAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAg                                                                                                        >  1:3105907/1‑150 (MQ=255)
                                                cggcgAACATGCAACAACGAAAACCGGCTTGGACCGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCAt                                                                                                     <  2:5646346/150‑1 (MQ=255)
                                                 ggctATCTTGCAAGAACTAAAACTGGCGAGTAACGGCGGGGGTCGCGCCCGACTTTGCTGGTTGAAAATGAGAATACCGAGCTTTTTTTTAACTTCCCTTGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCAtt                                                                                                    <  2:1013697/150‑1 (MQ=255)
                                                    gATCATGCAAGAACGAAAACAGGCGAGTAACGGCCGGGGTCGCGCCCGACTTTGCTGTTTGCAATTGAGCATACGGAGCCTTTTTTTAACTTAGCTTGATCGGCGCGCTGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGc                                                                                                 <  2:3656142/150‑1 (MQ=255)
                                                     aTCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCt                                                                                                >  1:8925917/1‑150 (MQ=255)
                                                      tcttGCAATAACTAAAACCGGCGCGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGAATGTTCGACGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTa                                                                                               <  2:3214564/147‑1 (MQ=255)
                                                        aTGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAat                                                                                             >  2:7244192/1‑150 (MQ=255)
                                                          gCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAatat                                                                                           >  1:5004560/1‑150 (MQ=255)
                                                          gCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAatat                                                                                           >  2:3293794/1‑150 (MQ=255)
                                                          gCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGAAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAatat                                                                                           >  2:8446233/1‑150 (MQ=255)
                                                           cAATAACTAAAACCGGCTTGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGAAAATGAGCATACCGCGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATAtt                                                                                          <  2:7660347/150‑1 (MQ=255)
                                                           cAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATAtt                                                                                          >  2:6825231/1‑150 (MQ=255)
                                                            aaGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTa                                                                                         >  2:4501259/1‑150 (MQ=255)
                                                            aaGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTa                                                                                         >  1:2507371/1‑150 (MQ=255)
                                                             agacctaaaaCCGGCGCGGCCCGGCGGGGGTCGCGCCCGACTTTGCTGGTTGCAATTGAGCATCCCTAGCCTTTTTTTACCGTAGCATGACCGCCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTAc                                                                                        <  1:6003343/144‑1 (MQ=255)
                                                               aaCCAAAACCGGCTTGGAACGGCGGGGGTCGCGCCCGACTTTGCTGTTTGAAAATGAGCATACCGAGCCTTTTTTTAACGTAGCATGATCTTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGa                                                                                      <  2:5845924/150‑1 (MQ=255)
                                                                  gAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGc                                                                                   >  1:3924843/1‑150 (MQ=255)
                                                                  gAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGc                                                                                   >  1:1558278/1‑150 (MQ=255)
                                                                  gAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGc                                                                                   >  2:4058019/1‑150 (MQ=255)
                                                                        cggCGAGGACCGGCGGGGGTCGCGCCCGACTTGGTTGGTTGAAAATGAGAATACCGAGCCTGTTTTTAAGGTAGCATGATGGTCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGc                                                                             <  2:4803466/149‑1 (MQ=255)
                                                                        aGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGc                                                                             >  1:6849410/1‑150 (MQ=255)
                                                                        aGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGc                                                                             >  2:7739143/1‑150 (MQ=255)
                                                                        aGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGc                                                                             >  2:8128980/1‑150 (MQ=255)
                                                                          gCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCa                                                                           >  1:914832/1‑150 (MQ=255)
                                                                          gCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCa                                                                           >  2:3258374/1‑150 (MQ=255)
                                                                            gAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATg                                                                         >  1:1967387/1‑150 (MQ=255)
                                                                            gAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATg                                                                         >  2:1522230/1‑150 (MQ=255)
                                                                              tgcaCGGGCGGGGTGCGCGCCCTAATTTGGTGTTTCAAATTCAGAATACCGTGCCTTTTTTTAACTTAGATTTATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGa                                                                        <  1:2756205/147‑1 (MQ=255)
                                                                               gAACGGCGGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGCGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAAGCCCTCGCCATGAcc                                                                      <  2:565419/150‑1 (MQ=255)
                                                                                accggcggGGGTCGCGCCCGACTTTGCTGGTTGCAATTGAGCAACCCGCGCCTTTTTTTACCGTCGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCg                                                                     <  2:1036071/148‑1 (MQ=255)
                                                                                 aCGGCCGGGTGCGCGCCCGACTTTGCTGTTTGAAAATGAGCATACCGAGCATGTTTTTAACGTAGAATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCgt                                                                    <  2:8856275/150‑1 (MQ=255)
                                                                                 aCGGCCGGGGTCGCGCCCTAATTTGCTTTTTTAAAATTACAATACATACCATTTTTTAAACTTACAATTATCTTCGCCCAGCCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCgt                                                                    <  1:5315492/150‑1 (MQ=255)
                                                                                 aCGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCgt                                                                    >  2:6141841/1‑150 (MQ=255)
                                                                                 aCGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCgt                                                                    >  2:6414378/1‑150 (MQ=255)
                                                                                    gCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGACCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTc                                                                 <  2:8699022/150‑1 (MQ=255)
                                                                                     cgggggTCGCGCCCGACTTTGCTGGTTGCAAATTAGCATACCGCGCCTGTTTTTAACGTAGCTTGATCGTGGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTcc                                                                <  2:8653574/148‑1 (MQ=255)
                                                                                      cGGGGACGCCCCCAAATTTGATGTTTCAAAATCACCAAACCCAGCATTTTTTTAACGTACATTTATCGTCGCGCAGCCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCa                                                               <  1:9055524/150‑1 (MQ=255)
                                                                                        gggTGGCGCCCGATTTTGCTGTTTTCAATTGTGAATCCGGCGCGTTTTTTTAACGTCGCTTGATGGGCGCGCCGCCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACt                                                             <  1:149320/150‑1 (MQ=255)
                                                                                        gggTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACt                                                             >  1:6863614/1‑150 (MQ=255)
                                                                                         ggtccgccCCCGACTTTGCTGTTTGAAATTGAGCATACCGCGCCTTTTTTTACCTTAGCTTTATCGTCGCGCCGCCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACt                                                             <  1:4243530/146‑1 (MQ=255)
                                                                                             cgcgCCCGACTTTGCTGGTTGCAAATGAGAATACCGAGCCTGTTTTTAACGTAGCTTGATGGTCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCaa                                                         <  2:1909704/150‑1 (MQ=255)
                                                                                                gcCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTg                                                      >  1:8036442/1‑150 (MQ=255)
                                                                                                 ccctattttGGTGGTTGAAATTTAGAATACCGAGCCTTTTTTTAAGGTAGCATGTTCGTCGCGCTGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGc                                                     <  1:8190978/144‑1 (MQ=255)
                                                                                                  ccGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGcc                                                    >  2:6538888/1‑150 (MQ=255)
                                                                                                   ccaatgggCTGGTTGCAAAGGCCCATCCCGTGCCTGTTTTTACCGTCGCAGGATCGTCGGCCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGccc                                                   <  1:1113947/143‑1 (MQ=255)
                                                                                                          gtgggTTGAAATTCAGCATCCCGCGCCTTTTTTTACCTTAGCTTTATGGGCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCaa                                            <  1:7772820/147‑1 (MQ=255)
                                                                                                            tGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAAcg                                          >  2:6027945/1‑150 (MQ=255)
                                                                                                              gTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAAcgcc                                        >  2:231618/1‑150 (MQ=255)
                                                                                                                tGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGAAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAAcgccgc                                      <  2:128818/150‑1 (MQ=255)
                                                                                                                 gCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCt                                     >  1:2347394/1‑150 (MQ=255)
                                                                                                                   aaaTGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGa                                   >  1:8433184/1‑150 (MQ=255)
                                                                                                                   aaaTGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGa                                   >  2:184867/1‑150 (MQ=255)
                                                                                                                   aaaTGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGa                                   >  2:1627086/1‑150 (MQ=255)
                                                                                                                   aaaTGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGa                                   >  2:1379129/1‑150 (MQ=255)
                                                                                                                      tGAGAATCCGGAGCGTTTTTTTAAGGTAGCTTTATCGTCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATcc                                <  1:5608688/150‑1 (MQ=255)
                                                                                                                       gAGCATACCGAGCCTTTTTTTAACGTAGAATTATCGTGGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCg                               <  1:6908234/150‑1 (MQ=255)
                                                                                                                            tagcGAGCGTGTTTTTAAGGTAGGTTGATCGTCGCGCAGGCAGTTTTCGCTCAATCCCCTAGTGCGCGTGAGCATTCGCTAATATTAGGAGGCCCTCGCCATGACCGTGTCCTCTGCAATTGCCCGTCTGAACGCCGCTGATCCGGAttt                          >  2:6468097/4‑150 (MQ=255)
                                                                                                                            tACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGAttt                          >  2:5160990/1‑150 (MQ=255)
                                                                                                                            tACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGAttt                          >  1:8763960/1‑150 (MQ=255)
                                                                                                                              ccGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCg                        <  2:3505663/150‑1 (MQ=255)
                                                                                                                                gAGCCTTTTTTTAACGTAGCTTGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGcc                      <  2:2598923/150‑1 (MQ=255)
                                                                                                                                gAGCCTTATTTTAACTAACCATGATCGTCGCCCAGCCACTATTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGcc                      <  2:7216955/150‑1 (MQ=255)
                                                                                                                                 aGCCTTTTTTTAACGTAGCATTATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGccc                     <  2:4385528/150‑1 (MQ=255)
                                                                                                                                 aGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGccc                     >  1:4760237/1‑150 (MQ=255)
                                                                                                                                 aGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGccc                     >  2:7415198/1‑150 (MQ=255)
                                                                                                                                   ccTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGa                   >  1:2376193/1‑150 (MQ=255)
                                                                                                                                    ctttttttACCTTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGAc                  <  1:4034927/147‑1 (MQ=255)
                                                                                                                                    cTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGAc                  >  1:3222133/1‑150 (MQ=255)
                                                                                                                                    cTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGAc                  >  1:4623414/1‑150 (MQ=255)
                                                                                                                                         tttAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTg             >  2:1788431/1‑150 (MQ=255)
                                                                                                                                         tttAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTg             >  1:972468/1‑150 (MQ=255)
                                                                                                                                           taagttAGCATGTTGGGCGCGCCGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGa           <  1:7855083/145‑1 (MQ=255)
                                                                                                                                           tAACGTTGCATGATCGTCGCGACGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGa           <  1:5921561/150‑1 (MQ=255)
                                                                                                                                               gTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGatcat       >  2:7647988/1‑150 (MQ=255)
                                                                                                                                               gTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGatcat       >  2:8902443/1‑150 (MQ=255)
                                                                                                                                               gTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGatcat       >  2:8996227/1‑150 (MQ=255)
                                                                                                                                                    aTGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATctgct  >  1:4413605/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCGCGCGGATGGACTGGTAGCTTGGGATAGTGCGAAAAGTCGCCGAGCGGCGATCATGCAAGAACGAAAACAGGCGAGGAACGGCCGGGGTCGCGCCCGACTTTGCTGGTTGCAAATGAGCATACCGAGCCTGTTTTTAACGTAGCATGATCGTCGCGCAGGCACTTTTCGCACAATCCCCTAGGGCGCCTGAGCATTCGCTAATATTACGAGGCCCTCGCCATGACCGTGTCCACTGCAATTGCCCGTCTCAACGCCGCTGATCCGGATTTCGCCCGACATCTGGATCATCTGCT  >  NC_002947/1105303‑1105599

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