Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,188,512 A→C L422* (TTA→TGA)  phoQ ← sensory histidine kinase in two‑component regulatory system with PhoP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,188,5120AC98.4% 248.0 / ‑3.5 64L422* (TTA→TGA) phoQsensory histidine kinase in two‑component regulatory system with PhoP
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (28/35);  total (28/36)
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.51e-01

CGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCG  >  NC_000913/1188365‑1188658
                                                                                                                                                   |                                                                                                                                                  
cgcgCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATg                                                                                                                                                 <  1:48756/151‑1 (MQ=255)
    cACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGaa                                                                                                                                              <  2:1006612/150‑1 (MQ=255)
    cACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGaa                                                                                                                                              >  2:805295/1‑150 (MQ=255)
        gCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCg                                                                                                                                         >  2:827013/1‑151 (MQ=255)
               aTGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCcat                                                                                                                                  >  2:859197/1‑151 (MQ=255)
               aTGCTCTCTCCGGCGACGATTTTACCCTCATATGGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCca                                                                                                                                   >  1:296903/1‑150 (MQ=255)
                  ctctctCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCcatcat                                                                                                                               >  1:422460/1‑151 (MQ=255)
                        ccGGCGACGATTTTACCCTCATATTGCTCCGTGATTTCCCGGGCTACCCCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCATCCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCAGTATCCTc                                                                                                                           <  2:429292/149‑1 (MQ=255)
                        ccGGCGACGATTTTACCCTCAGATTGCTCGGTGATTTCGCGGGCTGCCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTc                                                                                                                           <  2:32814/149‑1 (MQ=255)
                             gacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACcaca                                                                                                                    <  1:458618/151‑1 (MQ=255)
                             gacgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACcac                                                                                                                     >  2:383930/1‑150 (MQ=255)
                              acgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACcaca                                                                                                                    >  1:182805/1‑150 (MQ=255)
                               cgaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAat                                                                                                                  >  2:289301/1‑151 (MQ=255)
                                gaTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAata                                                                                                                 <  2:231001/151‑1 (MQ=255)
                                 aTTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAatat                                                                                                                >  2:611086/1‑151 (MQ=255)
                                  ttttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAatat                                                                                                                >  2:606158/1‑150 (MQ=255)
                                    ttACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATag                                                                                                              <  1:54772/150‑1 (MQ=255)
                                     tACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATaga                                                                                                             <  1:272046/150‑1 (MQ=255)
                                        ccTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATg                                                                                                         >  2:564682/1‑151 (MQ=255)
                                               ttGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCtcgtc                                                                                                   <  2:643983/150‑1 (MQ=255)
                                                 gCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGt                                                                                                >  1:1053988/1‑151 (MQ=255)
                                                 gCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGt                                                                                                <  1:785147/151‑1 (MQ=255)
                                                        gATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCTAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCtt                                                                                         >  1:594883/1‑151 (MQ=255)
                                                         aTTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCtt                                                                                         <  1:753420/150‑1 (MQ=255)
                                                         aTTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTg                                                                                        >  2:327703/1‑151 (MQ=255)
                                                         aTTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTg                                                                                        <  2:504037/151‑1 (MQ=255)
                                                                gggCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGaaa                                                                                  <  1:1025541/150‑1 (MQ=255)
                                                                gggCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAAt                                                                                 <  2:1149420/151‑1 (MQ=255)
                                                                gggCTACCGCCAGCCCTACAACTTGCCCAGGGCGTAAAGTAGCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGaa                                                                                   <  2:462329/149‑1 (MQ=255)
                                                                       cGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGAc                                                                           >  2:409270/1‑150 (MQ=255)
                                                                           aGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACt                                                                      <  1:1096512/151‑1 (MQ=255)
                                                                            gCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATAGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTc                                                                     <  1:897862/151‑1 (MQ=255)
                                                                             cccTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCg                                                                    >  2:289006/1‑151 (MQ=255)
                                                                                tACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAgg                                                                 >  2:261962/1‑151 (MQ=255)
                                                                                 acacCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGc                                                                >  2:823569/1‑151 (MQ=255)
                                                                                  cacCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCa                                                               <  1:75287/151‑1 (MQ=255)
                                                                                    ccTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCaa                                                              <  1:582316/150‑1 (MQ=255)
                                                                                          ccAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTa                                                        >  1:449844/1‑150 (MQ=255)
                                                                                          ccAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTAc                                                       >  2:305450/1‑151 (MQ=255)
                                                                                            aGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAg                                                     <  2:136677/151‑1 (MQ=255)
                                                                                               gCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCa                                                  <  1:149391/151‑1 (MQ=255)
                                                                                                cGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCAt                                                 <  1:1105950/151‑1 (MQ=255)
                                                                                                 gTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCAtt                                                >  2:298110/1‑151 (MQ=255)
                                                                                                      gTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCa                                           <  2:599115/151‑1 (MQ=255)
                                                                                                      gTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCa                                           <  2:736044/151‑1 (MQ=255)
                                                                                                        aTCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAg                                          <  1:36899/150‑1 (MQ=255)
                                                                                                         tCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGc                                         >  1:474753/1‑150 (MQ=255)
                                                                                                          cAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCa                                        <  1:758233/150‑1 (MQ=255)
                                                                                                          cAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCAc                                       <  1:946184/151‑1 (MQ=255)
                                                                                                                   gACCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCa                               <  1:338559/150‑1 (MQ=255)
                                                                                                                    aCCACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATc                             >  1:860044/1‑151 (MQ=255)
                                                                                                                     ccACGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCa                            <  1:478774/151‑1 (MQ=255)
                                                                                                                       aCGGTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCAcc                          >  1:915711/1‑151 (MQ=255)
                                                                                                                         ggTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTc                        >  2:1048502/1‑151 (MQ=255)
                                                                                                                         ggTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTc                        <  2:687839/151‑1 (MQ=255)
                                                                                                                          gTCGAAAATGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCg                       <  1:405152/151‑1 (MQ=255)
                                                                                                                                aaTGACCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACaaaa                 <  1:614478/151‑1 (MQ=255)
                                                                                                                                    aCCTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGt             <  1:841654/151‑1 (MQ=255)
                                                                                                                                     ccTCTCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGtt            >  2:915791/1‑151 (MQ=255)
                                                                                                                                     ccTCTCGCTAGCTTCATGGAATACCGGTGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTGCAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAAACgat            <  2:637936/151‑3 (MQ=255)
                                                                                                                                      ctctCGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGTTc           >  2:1033497/1‑151 (MQ=255)
                                                                                                                                        ctcGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGTTc           <  2:1148002/149‑1 (MQ=255)
                                                                                                                                          cGCTTGCTTCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGTTCTg         <  1:146147/149‑1 (MQ=255)
                                                                                                                                               gcttCATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCg  >  2:699261/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCGCACCGCCCAGCATGCTCTCTCCGGCGACGATTTTACCCTCATATTGCTCGGTGATTTCGCGGGCTACCGCCAGCCCTACACCTTGCCCAGGGCGTAAAGTATCAACCCGTTGACCACGGTCGAAAATGACCTCTCGCTTGCTTAATGGAATACCGGGGCCATCATCCTCGACCACAATATAGAGATGCTCGTCGGTTTGCCTTGCAGAAATTTCGACAAACTCGAGGCAATATTTACAGGCATTATCCAGCACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCG  >  NC_000913/1188365‑1188658

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