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,8610CT97.4% 104.6 / ‑6.0 39V306I (GTA→ATA) phoQsensor histidine kinase PhoQ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/11);  total (28/11)
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.17e-01

CACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGACTTCTTTTGC  >  NC_000913/1188619‑1189104
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
cACGTTGCCCATCCCCTCTACAAAATCGTCATGCTCCCCGACAAAGGTGATCTCGGGCGAAATATCGAGAGAGTTATTGTCCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTTTCCAGCAGTGGGGCGACCGGAGGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAg                                                                                                                                                                                                                                              <  2:118801/251‑1 (MQ=255)
   gTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCAt                                                                                                                                                                                                                                           >  1:371221/1‑251 (MQ=255)
                tCGACAAAATCGTACTGGTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTACAGCAGTGCGGCGACCGTATGCAGCTCACGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCCGTGACACGCGGCTGATTTGCTCCAGCATTATCGGCTCCGCCTCACTTACGCTCac                                                                                                                                                                                                                              >  2:311930/1‑250 (MQ=255)
                 cGACAAAATAGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTTATACACTTTGTTCAACGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGTTGCCGGAAGTAGCCACTTTGCTGGGAAATGCGGCTGATTTGTTCCAGCATTATCGGCTCCGCATCAATTACGCTCATc                                                                                                                                                                                                                             >  2:556399/1‑251 (MQ=255)
                  gACACAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGCAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGAGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGCTGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCGTTATCGGCTCAGCATCACTGACGCTCATCt                                                                                                                                                                                                                            >  2:86491/1‑251 (MQ=255)
                    cAAAATAGTTTTGCTCACCGAAAAAGCTGATCTCTGGCGAAATATCGAGAGAGATTTTGTCCCGTGTGTGTTGATACTCTTTGTGCAGCGCTGAGGTGAGATTGTCCATCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCttt                                                                                                                                                                                                                          <  2:515114/251‑1 (MQ=255)
                      aaaTCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTc                                                                                                                                                                                                                        <  1:311930/251‑1 (MQ=255)
                        aTCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCAc                                                                                                                                                                                                                      >  1:19186/1‑251 (MQ=255)
                         gcGTTCTGCTCACCGACAAAGCTGATCTCTGGCTAAATAGCGAGAGAGATATTGACCCCTTTGCGTTGATACACGTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACt                                                                                                                                                                                                                     <  2:281661/250‑1 (MQ=255)
                           gTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCCTCACTGACGCTCATCTTTTCACTAc                                                                                                                                                                                                                   >  2:399824/1‑251 (MQ=255)
                                           aaGCTGATATCTGGCGAAATATGGAGAGAGATATTGACCCCTTTGCGTTGCTACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGAGAGACCGGATGCAGCGCGCGGCTGAGCAATGTCCCGCCGGG‑CATAATGGCACGATGCAGGAAGGAGCCAATTTGCTGTGAAATGAGACTGATTTGCTCCAGCATTATCGGCTCGGCATCAATAAAGCACATATTTTCACTACG‑AGAAGAACGCAGCg                                                                                                                                                                                                   >  2:183710/1‑250 (MQ=255)
                                                    tctGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCCGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTgc                                                                                                                                                                                           >  2:301110/1‑249 (MQ=255)
                                                     ctGGCGAAATATCGAGAGAGTTCTTGACCCCTTTGCGGTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATCCTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGCAATGCGCCTGATTTGCTCCCGCATTATCGGCTCAGCATCACTGACGCTACTCTTTTCACTACGCAGA‑GAACGCAGCGTACtctgaagc                                                                                                                                                                                        >  2:347206/1‑244 (MQ=255)
                                                        gCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGACCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCCCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCAcc                                                                                                                                                                                     <  2:513999/251‑1 (MQ=255)
                                                          gAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGc                                                                                                                                                                                   <  1:399824/251‑1 (MQ=255)
                                                          gAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGc                                                                                                                                                                                   >  1:109930/1‑251 (MQ=255)
                                                             atatCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGAGAGATTGTCCAGCAGTGGTGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCCGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCCATGACGCTCAACTTTTCACTACGCAGA‑GACCGCAGCGTACTTTGCAGCACCGccc                                                                                                                                                                                 >  2:335958/1‑249 (MQ=255)
                                                              tatCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGAGTACTTTGCAGCACCGCCAGt                                                                                                                                                                               >  1:329212/1‑251 (MQ=255)
                                                                                caccTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGGACTTTGCAGCACCGCCAGTGGAGTTTTCAGACTATgg                                                                                                                                                             >  2:364877/3‑251 (MQ=255)
                                                                                   cTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTc                                                                                                                                                          >  1:336286/1‑251 (MQ=255)
                                                                                             aTACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTgg                                                                                                                                                >  1:1887/1‑250 (MQ=255)
                                                                                             aTACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTg                                                                                                                                                 >  1:353045/1‑250 (MQ=255)
                                                                                                    tttttctggcgtgtggTGATATTGCCCCGCGGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAAGTTCCCGCCGCT‑CATAGTGGCACGATGAGGGTAGTCGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCTTTATCGGCTCAGCACCACTGACCCTAATCTTTTCACTACGCAGA‑GACCGCAGCGTTCTTTGCACCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGcgtcg                                                                                                                                          <  1:570433/236‑1 (MQ=255)
                                                                                                                  ggTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGAATTATCGGCTCAGCCTCACTGACGCTCATCTTTTCACTACGCCGC‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGAGTCGTACGGGATTTGTCGt                                                                                                                           >  2:328861/1‑251 (MQ=255)
                                                                                                                               agcagTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGAATTTGTCGTAACGTTCGCGTTc                                                                                                              >  1:400750/1‑251 (MQ=255)
                                                                                                                                          gACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGCTTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GTACGCCGCGTACTTTGCAGCACCGCCAGTGGCGTTGTCAGACTATGGGTCTGGTCGGTGAGCGTCGTACGGTATTTGTTGTAACGTTCGCGTTCACTTTTTAAc                                                                                                    >  2:90129/1‑250 (MQ=255)
                                                                                                                                               gATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGGAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCCGCGTACTTTGCAGCACCGCCAGTGGCGTTATCAGCCTATGGGTCAGATCGGTGAGCGTCGTTCGGGATTTGTGGTAACGTTTGCGGTCACTTTTTAACAAtagg                                                                                              >  2:427616/1‑248 (MQ=255)
                                                                                                                                                        cgcgGCTGAGCAAGTTCCCGCCGCG‑CATCCTGGCACGATGCAGGTAGTAGCCAATTTGCTGGAAAATGGGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTCCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGttt                                                                                     <  1:335958/251‑1 (MQ=255)
                                                                                                                                                               gAGCAATGTCCCGCGGCG‑CATACGGGCACGAGGCAGGTGGGAGCCATTTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCCTCCCTGACGCTCATCTTTTCACTACGAAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCa                                                                              <  1:216917/251‑1 (MQ=255)
                                                                                                                                                                    aTGTCCCGC‑GCGGCATCCTGGCACGTGGCAGTTATTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTATCGCCCCAGCATCACGGACGCTCATCTTTTCCCTCCGCAGA‑GAACGCAGCGTATTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTg                                                                         <  1:86491/251‑1 (MQ=255)
                                                                                                                                                                                                 cgttttgtagcaaattgtcTGTGAAAGGGGGCTGATTTGCTCCAGATTTATCGGCTCATCATCACGGACGCTCATCTTTTCCCTACGCAGA‑GAACGCGGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGa                                             <  1:20574/233‑1 (MQ=255)
                                                                                                                                                                                                       gtagCCAATTTGCTGTGAAATTCGGCTGATTTGCTCCAGCATTATCGGCTCAGCCTCACTGACGCTCATCTTTTCACTCCGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGATCGTCGTACGGAATTTGTCGTAACGATCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGGCTGACCAGTTCTCACGTTGTGGCTGGATTggagct                                         >  2:20574/1‑244 (MQ=255)
                                                                                                                                                                                                           tcactttGCTGTGAAATGCGGCTGATTTGCTCCAGAATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGAACATTGCAGCCCCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTAGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCAGTTGTTAACACTCGGGACAGGTTTCTTACCAGACTGGGCAGTTTTCGAGTTGTGGCCGGAATGAgaacttacgg                                   >  2:201757/5‑242 (MQ=255)
                                                                                                                                                                                                                   cTGTGAAAAGCGGCTGATTTGATCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑AAACGCAGAGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTAAGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTAGTAACGTTCGCGTTCAATTTTTAACAATCGGTTCAGGTTTCGTACAAGACTGGTCAGTTCACGCCTTGTAGCTGGGTTGAGCCATTCGCGGttgttgt                           >  2:280398/1‑246 (MQ=255)
                                                                                                                                                                                                                      tgAAATGCGGCTGATTTGCTCCAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGttctt                        >  1:281719/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ccAGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGATATGTTCTTCCAGTTCGCGGACTTCtttt    >  1:43173/1‑251 (MQ=255)
                                                                                                                                                                                                                                          ccAGCATTACAGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTAATTTGCAGCACAGCACGTGGCGTTTTCAGACTATGGGTCAGCTCGGTGAGCGTCGTACGGAATTTGTCGAAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGATTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGACTTCtttt    >  1:64826/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cAGCATTATAGGCTCAGCATCACTGACGCTCCTCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCATCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTAACAGACTGGTCAGTTCTCGCCTTGTGGCTGGTTTGAGGAATTAGCGGATAGGTTCTTCCacatcccggcattcttttg   >  2:312292/1‑233 (MQ=255)
                                                                                                                                                                                                                                            aGCATTATCGGCTCAGCATCACTGACGCTCATCTTTTCACTCCGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGCCTGGTCAGTTCTCGCGTTGTGGCTTGATTGAGCAATTCGCGGTTATGTTCTTCCATTTCGCGGCCTTCTTTTGc  >  2:5004/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                    
CACGTTGCCCATCACCTCGACAAAATCGTTCTGCTCACCGACAAAGCTGATCTCTGGCGAAATATCGAGAGAGATATTGACCCCTTTGCGTTGATACACTTTGTTCAGCGCTGAGGTGAGATTGTCCAGCAGTGGGGCGACCGGATGCAGCTCGCGGCTGAGCAATGTCCCGCCGCG‑CATACTGGCACGATGCAGGTAGTAGCCAATTTGCTGTGAAATGCGGCTGATTTGCTCCAGCATTACCGGCTCAGCATCACTGACGCTCATCTTTTCACTACGCAGA‑GAACGCAGCGTACTTTGCAGCACCGCCAGTGGCGTTTTCAGACTATGGGTCAGGTCGGTGAGCGTCGTACGGTATTTGTCGTAACGTTCGCGTTCACTTTTTAACAATCGGTTCAGGTTTCGTACCAGACTGGTCAGTTCTCGCGTTGTGGCTGGATTGAGCAATTCGCGGTTATGTTCTTCCAGTTCGCGGACTTCTTTTGC  >  NC_000913/1188619‑1189104

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