Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 632,082 C→A 5.6% intergenic (+83/‑100) cstA → / → ybdD carbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913632,0820CA 5.6% 308.5 / 3.4 88intergenic (+83/‑100)cstA/ybdDcarbon starvation protein involved in peptide utilization; APC peptide transporter family protein/DUF466 family protein
Reads supporting (aligned to +/- strand):  ref base C (39/44);  new base A (2/3);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

CCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGC  >  NC_000913/631937‑632228
                                                                                                                                                 |                                                                                                                                                  
ccGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                                                                                                                                                 <  5:56620/149‑1 (MQ=255)
 cGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAgg                                                                                                                                                >  1:387855/1‑149 (MQ=255)
  gTATGACCCAATGCCGGAAAATGTCGAGGCGAGCGGGGCGCAGGCAAAAGGCGCACACTAAAGGCAGAGGGAGGGGGGGATGGTGGCGAATGTGGGCTTAGTGCCCAGGGTGCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                                                                                                                                               <  1:252403/149‑1 (MQ=255)
  gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                                                                                                                                               <  5:8322/149‑1 (MQ=255)
  gTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGGTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAggg                                                                                                                                               <  1:246105/149‑1 (MQ=255)
     tGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                                                                                                                                            <  1:250338/149‑1 (MQ=255)
          cAATGCCGGAAAATGTCGAGGAGATCGTGGCGCCGGCAAAAGGCGCACACTAAAGTCAGAGTCAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGCGTTCCCTCTCAACCTAACCCTATCCACGGTGGGGCGAGGGGACTGAcc                                                                                                                                       >  5:181377/1‑149 (MQ=255)
             tGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAg                                                                                                                                    >  7:9381/1‑149 (MQ=255)
              gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                                                                                                                                            >  3:285578/1‑140 (MQ=255)
              gCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGAc                                                                                                                                            <  4:285578/140‑1 (MQ=255)
               ccGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgc                                                                                                                                   <  6:321318/148‑1 (MQ=255)
                cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcgc                                                                                                                                 >  2:123767/1‑149 (MQ=255)
                cGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAgcgc                                                                                                                                 >  6:149131/1‑149 (MQ=255)
                  gAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCCCCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGCGCGCGt                                                                                                                               >  1:10711/1‑149 (MQ=255)
                   aaaaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGtt                                                                                                                              <  4:351152/149‑1 (MQ=255)
                     aaTGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGa                                                                                                                            <  8:381007/149‑1 (MQ=255)
                       tggcGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATa                                                                                                                          <  4:81335/146‑1 (MQ=255)
                           gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCa                                                                                                                       <  3:43728/148‑1 (MQ=255)
                           gaggagATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCAt                                                                                                                      <  2:394963/149‑1 (MQ=255)
                              gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                                   <  7:202011/149‑1 (MQ=255)
                              gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                                   <  3:165441/149‑1 (MQ=255)
                              gagaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGGGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                                   <  3:160634/149‑1 (MQ=255)
                               agaTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                                  <  4:153511/149‑1 (MQ=255)
                                 aTCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAg                                                                                                                <  2:387855/149‑1 (MQ=255)
                                          cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGc                                                                                                                        >  7:39248/1‑132 (MQ=255)
                                          cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGc                                                                                                                        <  8:39248/132‑1 (MQ=255)
                                          cAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATaa                                                                                                       <  8:9381/149‑1 (MQ=255)
                                               aaaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                                                                                                                                     >  6:414726/1‑114 (MQ=255)
                                               aaaaGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGa                                                                                                                                     <  5:414726/114‑1 (MQ=255)
                                                  aGGCGCACACTAAAGTCAGAGGGAGTGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCGCGGGGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCa                                                                                               <  6:181377/149‑1 (MQ=255)
                                                         cacTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                                  <  4:267934/123‑1 (MQ=255)
                                                         cacTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGt                                                                                                                  >  3:267934/1‑123 (MQ=255)
                                                         cacTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCa                                                                                        >  7:129129/1‑149 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                                                                                                                                             >  3:196674/1‑95 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGa                                                                                                                                             <  4:196674/95‑1 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTc                                                                                                <  1:266460/140‑1 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTc                                                                                                >  2:266460/1‑140 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTc                                                                                                <  1:365170/140‑1 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTc                                                                                                >  2:365170/1‑140 (MQ=255)
                                                          acTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACACTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTAGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                       >  2:247524/1‑149 (MQ=255)
                                                            tAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                                   <  3:311501/119‑1 (MQ=255)
                                                            tAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTg                                                                                                                   >  4:311501/1‑119 (MQ=255)
                                                            tAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                     >  2:149142/1‑149 (MQ=255)
                                                            tAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAggg                                                                                                                                               >  4:232296/1‑91 (MQ=255)
                                                            tAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCGCCCCACCGGGGAGAggg                                                                                                                                               <  3:232296/91‑1 (MQ=255)
                                                                gTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGa                                                                                                          <  8:366902/124‑1 (MQ=255)
                                                                gTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGa                                                                                                          >  7:366902/1‑124 (MQ=255)
                                                                gTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                 <  5:232049/149‑1 (MQ=255)
                                                                gTCAGAGTGAGGGGGGGATGGTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                 <  1:149142/149‑1 (MQ=255)
                                                                   actcggcgggggggATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCActc                                                                              <  4:368658/142‑1 (MQ=255)
                                                                     agTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                                                                                                                                                 <  4:283109/80‑1 (MQ=255)
                                                                     agTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAg                                                                                                                                                 >  3:283109/1‑80 (MQ=255)
                                                                       tGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTc                                                                          <  7:219776/149‑1 (MQ=255)
                                                                                gATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATg                                                                 >  2:428419/1‑149 (MQ=255)
                                                                                  tGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGt                                                               >  1:358370/1‑149 (MQ=255)
                                                                                    ttGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGtt                                                              <  8:129129/148‑1 (MQ=255)
                                                                                       gCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACa                                                                   <  1:98031/140‑1 (MQ=255)
                                                                                       gCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACa                                                                   >  2:98031/1‑140 (MQ=255)
                                                                                         gAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAg                                                                        >  7:4437/1‑133 (MQ=255)
                                                                                         gAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAg                                                                        <  8:4437/133‑1 (MQ=255)
                                                                                                   ttAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATg                                              >  6:340581/1‑149 (MQ=255)
                                                                                                   ttAGTGCCCAGGGTTCACTCTAACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGTATAAGGCGTTACAGCAGAATCCGGAACTCTTTCAGCAACATGTTTAGCTGAGGCCAAGATg                                              >  2:157592/1‑149 (MQ=255)
                                                                                                      gTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                                           <  6:41952/149‑1 (MQ=255)
                                                                                                       tGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                                           >  2:305962/1‑148 (MQ=255)
                                                                                                         cccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                        >  7:320692/1‑149 (MQ=255)
                                                                                                         cccAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGAt                                        >  7:368810/1‑149 (MQ=255)
                                                                                                           cAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTc                                      <  5:149131/149‑1 (MQ=255)
                                                                                                            aGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCa                                     >  7:103702/1‑149 (MQ=255)
                                                                                                            aGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCa                                     <  1:305962/149‑1 (MQ=255)
                                                                                                               gTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTg                                  <  5:340581/149‑1 (MQ=255)
                                                                                                                ttCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTgg                                 >  8:27940/1‑149 (MQ=255)
                                                                                                                ttCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTgg                                 <  1:428419/149‑1 (MQ=255)
                                                                                                                     tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGtt                                            <  5:61377/133‑1 (MQ=255)
                                                                                                                     tctcACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGtt                                            >  6:61377/1‑133 (MQ=255)
                                                                                                                          cccTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCg                        >  8:304984/1‑148 (MQ=255)
                                                                                                                             tAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaa                     >  5:33900/1‑148 (MQ=255)
                                                                                                                              tccccTCGCCCCACCGGGGAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaaa                   >  4:211443/3‑149 (MQ=39)
                                                                                                                              aaCCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGaaaa                   >  5:63639/1‑149 (MQ=255)
                                                                                                                                 cccctcCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAatat                >  1:413924/4‑149 (MQ=255)
                                                                                                                                 ccTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                                    <  2:109495/65‑1 (MQ=255)
                                                                                                                                 ccTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGc                                                                                                    >  1:109495/1‑65 (MQ=255)
                                                                                                                                      ccccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAgg           >  1:201646/1‑149 (MQ=255)
                                                                                                                                       cccaccggggAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGt          <  3:211443/143‑1 (MQ=255)
                                                                                                                                       cccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGt          <  7:304984/149‑1 (MQ=255)
                                                                                                                                        ccaccggggAGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTc         <  7:158706/144‑1 (MQ=255)
                                                                                                                                        ccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                                           >  6:309233/1‑115 (MQ=255)
                                                                                                                                        ccGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGttt                                           <  5:309233/115‑1 (MQ=255)
                                                                                                                                                gCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAggcggc  >  4:119377/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCGTATGAACCAATGCCGGAAAATGTCGAGGAGATCGTGGCGCAGGCAAAAGGCGCACACTAAAGTCAGAGTGAGGGGGGGATGTTGGCGAATGTTGGCTTAGTGCCCAGGGTTCCCTCTCACCCTAACCCTCTCCCCGGTGGGGCGAGGGGACTGACCGAGCGCGTTGATAGCATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACTCTTTCAGCAACATGGTTAGCGGAGGCCAAGATGTTTGATTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGC  >  NC_000913/631937‑632228

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: