Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,188,861 C→T V306I (GTA→ATA)  phoQ ← sensor histidine kinase PhoQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,188,8610CT95.8% 138.1 / ‑6.0 51V306I (GTA→ATA) phoQsensor histidine kinase PhoQ
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (27/19);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGACTTCTTTT  >  NC_000913/1188622‑1189102
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
gTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGATATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATGGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTTATTTGCTCCAGCATTAGCGGCTCAGCAt                                                                                                                                                                                                                                        <  2:198267/251‑1 (MQ=255)
  tGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGCGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGCGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATc                                                                                                                                                                                                                                       >  1:437338/1‑250 (MQ=255)
   gCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAc                                                                                                                                                                                                                                     >  1:92683/1‑251 (MQ=255)
   gCCCATCACATCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGCGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCAc                                                                                                                                                                                                                                     >  1:396151/1‑251 (MQ=255)
             tCGACAAAATCGTTCTGCTCACCGACAAAGCTGCTCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGCGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCAt                                                                                                                                                                                                                           >  2:93943/1‑251 (MQ=255)
                    aaTCGTTCTGCTCACCGACAAAGCTGCTCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTTGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGCTGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCCTCACTGACGCTCATCTTTTc                                                                                                                                                                                                                     >  1:545664/1‑250 (MQ=255)
                          tCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATAATGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACg                                                                                                                                                                                                               <  2:465444/250‑1 (MQ=255)
                             gCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGAAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGCCGCTCATCTTTTCACTACGc                                                                                                                                                                                                              >  2:177125/1‑248 (MQ=255)
                                 aCCGACAAAGCTGATCTCTGGAGAAATATCGAGAGAGATTTTGACCCCTTTGCGTTGATACAATTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTAGCGGCCGAGCAATGTCACGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTACAGCATTAACGGATCAGCATAACTAACGATCAGCTTTTCACTACACagtaaac                                                                                                                                                                                                       >  2:512650/1‑246 (MQ=255)
                                   agACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGc                                                                                                                                                                                                     >  2:330108/2‑251 (MQ=255)
                                      cAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGc                                                                                                                                                                                                  <  1:335257/251‑1 (MQ=255)
                                       aaaGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGAAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCg                                                                                                                                                                                                 >  2:86115/1‑251 (MQ=255)
                                        aaGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGCCCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCGGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTAAGCATCACTGCCGCTCATCTTTTCACTACGCAGAGAACGCAGCGt                                                                                                                                                                                                >  2:469142/1‑251 (MQ=255)
                                              aTCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTg                                                                                                                                                                                          >  2:469647/1‑251 (MQ=255)
                                                       gAAATATCGAGAGAGATATTGACCCCTTTGCTTTGATACACTTTGGTCATCGCTGAGGTGAGATTGTCCAGCCGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTGGCAGCACCGc                                                                                                                                                                                 <  2:404395/251‑1 (MQ=255)
                                                        aaaTATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGcc                                                                                                                                                                                >  2:443765/1‑251 (MQ=255)
                                                          atatCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAg                                                                                                                                                                              <  1:196829/251‑1 (MQ=255)
                                                          atatCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAg                                                                                                                                                                              <  1:99172/251‑1 (MQ=255)
                                                                     gaTATTGCCCCCTTGGCGTTGATCCCCTTGGTTCAGCGCGGAGGGGAGTTTGCCCAGCGGGGGGGCGCCCGGAGCCAGCTCGCGGCGGAGCAATGCCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTGGCTCCAGCATTATCGGCGCAGCATCCCTGACGCCCATCTTTTCACTACGCAGAGAACGCAGCTTACTTGGCAGCACCGCCAGTGGCGTTTTCa                                                                                                                                                                   <  1:86115/251‑1 (MQ=255)
                                                                         ttGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGAGAGATTGTCCCTCAGTGGGACGACCGGGTGCCGCTCGCGGCTGAGCAATGTCCCGCCGCGCATAATGGCACTATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTTACGCTCCTCTTTTCACTACGCCGATAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCataat                                                                                                                                                               >  2:444919/1‑247 (MQ=255)
                                                                          tGACCCCTTTGGGTGGATACACGTGGTTAAGCGTGGGGGTGAGATTGGCCAGCGGGGGGGCGCCCGGATGCAGCCCGCGGCTGAGCAAGGCCCCGCCGCGCATCCTGGCACGATGATGGTAGTAGCCAATTTTCGGGGAAAGGCGTCTGATTTGCCCCAGCATTATGGGCTCAGCATCACGGACGCTCATCTTTTCACTCCGAAGAGATCGCATCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACt                                                                                                                                                               <  2:563777/250‑1 (MQ=255)
                                                                          tGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCCGCACCGCCAGTGGCGTTTTCAGACt                                                                                                                                                               >  1:503813/1‑250 (MQ=255)
                                                                                  ttGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTc                                                                                                                                                        >  1:297831/1‑249 (MQ=255)
                                                                                  ttGCGTTGATACACCTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCCTCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCTGTGGCGTTTTCAGACTATGGGTc                                                                                                                                                        >  2:328452/1‑249 (MQ=255)
                                                                                   tGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGACCGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGCCTATGGGTCAgg                                                                                                                                                     >  2:94678/1‑251 (MQ=255)
                                                                                     cGTTGATACCCTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGTATGCAGCTCGCGGCTGAGAAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGTTTATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTc                                                                                                                                                   <  1:94678/251‑1 (MQ=255)
                                                                                      gTTGATACACTTTGTTCAGCGCTGTGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTAAGCATCACTGACGCTCATCTTTTCACTTCGCCGACAACGCAGCGTACTTTCCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCCGGTCg                                                                                                                                                  >  2:22800/1‑251 (MQ=255)
                                                                                           cacacTTTGTTCAGCGCGGAGTTGAGATTGTCCAGCAGTGGGGCGCCCGGAGGCAGCCCGCGGTTGAGCAATGCCCCGCCGCGCTTATTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAAGGCGGCTTATTTGCTCCAGCATTACCGGCGCAGCATCCCTGACGCTCATCTTTCCACTACGCAGAGACCGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAg                                                                                                                                             <  1:240654/250‑1 (MQ=255)
                                                                                            ccgcTTTGTTAAGCGCTGAGGTGAGATTGTCCGGCGGGGGGGCGACCGGATGCAGCCCGCGGCTGAGCATTGCCCCGCCGCGCATACGGGCACGATGCAGGTAGTAGCAAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATTTTTTCACTACGCAGAGAACGCAGCGTATTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGc                                                                                                                                            <  1:330108/248‑1 (MQ=255)
                                                                                            ccatttgtttCAGCGGTGGGGGGAGATGGTCCAGCAGTGGGGCGACCGGATGCAGCCCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCATTTTGCTGTGAAATGCGGCTGATTTGCCCCAGCATTATCGGCCCAGCATCACTGACGCTCATCTTTTCACTCCGCAGAGAACGCGGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGc                                                                                                                                            <  2:203245/243‑1 (MQ=255)
                                                                                            acacTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGGGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTGGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCTTACTTGGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGc                                                                                                                                            <  1:195258/251‑1 (MQ=255)
                                                                                               cTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGGAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGcgtc                                                                                                                                         >  2:341497/1‑251 (MQ=255)
                                                                                                     tCAGCGCTGAGGTGAGATTGTCCAGCGGTGGGGCGACCGGATGCAGCTCGCGCCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCACGGACCCTCATTTTTTCCCGCCGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACg                                                                                                                                    <  1:469142/250‑1 (MQ=255)
                                                                                                           gtgaggggggatgttCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGTTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACCCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTg                                                                                                                             <  1:498654/237‑1 (MQ=255)
                                                                                                           cTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGCGCGTCGTACGGTATTTg                                                                                                                             >  1:547157/1‑251 (MQ=255)
                                                                                                             ggggggagttttcccggcgggggggCGCCCGGATGCCGCCCGCGTCTGAGCAAGGCCCCGCCGCACATCCTGGCACGATGCGGGTAGTAGCCAATTTGTTGGGAAATGCGGCTGATTTGCCCCAGCATTATCGGCTCAGCATCACGGACGCTCATCTTTTCCCAACGCAGAGAACGCGGCGTTCTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTc                                                                                                                           <  2:92683/230‑1 (MQ=255)
                                                                                                                  gagaTTGTCCAGCAGTGGAGAGCCCGGATGCAGCTCGCGGCTGAGCAATGTCACGCCGCGCATACTGGCACGCTGCGGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCCCTCCGCAGAGACCGCAGCATACTTTGCCGCACCGCCCGTGGCGTTTTCAGCCTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTTTCGAAAc                                                                                                                      >  2:360233/1‑251 (MQ=255)
                                                                                                                            gccggggggccgccgggATGCAGCCGGCGGCTGAGCAATGCCCCGCCGCGCATCCTGGCCCGATGCGGGAGGTAGCCAATTTGCTGTAAAATGCGGCTTATTTGCCCCAGCATTATCGGCGCAGAATCACTGACGCTCATCTTTCCACTCCGCAGAGAACGCAGCGTTCTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTc                                                                                                            <  2:155393/236‑1 (MQ=255)
                                                                                                                                     ggggccgtaTGCAGCTCGCGGCTGAGCAATGCCCCGCCGCGCATCCTGGCCCGTTGCAGGTAGTAGCCAATTTGCTGGGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTGCGCTCCGCAGAGAACTCAGCGGTCTTTGCGGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGGCGTAACGTTCGCGTTCACttttcaa                                                                                                   <  2:190363/243‑4 (MQ=255)
                                                                                                                                     gcgctcggagcacgcgcgcgGCTGAGCAATGTCCCCCCGCGAATACTGGCTTGATGAAGGTAGTATCCAATTTGGTGTGAAATGCGGCTTATTTTCTCCAATATATGCGGCTCGGCTGCACTGCCGCGCGTCTTTTCACTCCCGAGAGAACGTAGCTTGCTTGGCAGCCCCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTaa                                                                                                   <  2:408992/235‑1 (MQ=255)
                                                                                                                                       gACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAAc                                                                                                  >  1:211016/1‑250 (MQ=255)
                                                                                                                                             gtGCAGCTCGCGGCTGAGCAATGTCCGGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAAGGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTCCTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGt                                                                                           <  1:443765/250‑1 (MQ=255)
                                                                                                                                                cAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCCCGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTc                                                                                         >  2:42321/1‑250 (MQ=255)
                                                                                                                                                            gAGCAAGGTCCCGCCGCGCATACGGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAAGGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCTTCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCTCCGCCAGTGGCGTTTTCAGACTATGTGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTAcc                                                                             <  1:93943/250‑1 (MQ=255)
                                                                                                                                                                     cccGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTcc                                                                   >  1:135379/1‑250 (MQ=255)
                                                                                                                                                                                       ccGATGGAGGTAGTAGCAAATTTGTGGTGAAATGCGGCTTATTGGCTCCAGCTTTATCGGCTCAGCATCACGGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTATTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTGTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTg                                                 <  1:22800/250‑1 (MQ=255)
                                                                                                                                                                                             cAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTg                                            >  1:472667/1‑250 (MQ=255)
                                                                                                                                                                                                                    caaTGCGGCTGAGTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTg                                                                                                                                                                                          <  1:333249/84‑1 (MQ=255)
                                                                                                                                                                                                                                tttGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCCGACTATGGGTCAGGTCGGTGAGCGACGTACGGAATTTGTCGTAACGTTCGCGTTCCCTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCCGTTCGCGGAct        >  1:300328/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGActtct     >  2:489200/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ccAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGACTTCtttt  >  2:583526/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
GTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCGCATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGAGAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGACTTCTTTT  >  NC_000913/1188622‑1189102

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