Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 3,943,262 G→A G302D (GGC→GAC)  LY180_18795 → bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065843,943,2620GA93.8% 86.5 / NA 32G302D (GGC→GAC) LY180_18795bifunctional (p)ppGpp synthetase II/ guanosine‑3',5'‑bis pyrophosphate 3'‑pyrophosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (18/12);  total (20/12)
Fisher's exact test for biased strand distribution p-value = 5.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGAGATGATCCAAAAAATCCTCTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGGCTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCT  >  CP006584/3942985‑3943556
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
aaaGAGATGATCCAAAAAATCCTCTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGCTGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGGCACAGAC‑ATAACAGCCATTCCAAACCCGAACGGTTATCatgttttccaccccc                                                                                                                                                                                                                                                                                       >  2:321967/1‑284 (MQ=255)
               aaaaTCCTCTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCACTCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGc                                                                                                                                                                                                                                                                      >  1:115017/1‑298 (MQ=255)
                      tcgttGAATTCGAAGGGGGTTGGCAGGAAGCGGGATAACCGTGCCGGGTAAGTGGTCGCAAAAAGCAGTTTTATTCGCATTACTGAAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCGATCATGGACATCTACGCTGTCCGCGTGATCGTCAATGATTCGGACACCTGTCATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGt                                                                                                                                                                                                                                                               <  2:295325/294‑1 (MQ=255)
                         ctGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGcc                                                                                                                                                                                                                                                            <  1:45295/298‑1 (MQ=255)
                             aaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAACTGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCAGCGTGACCGCCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGACTGTACAAGCCCCCTCCGGGCCGCGTGCAAGCCTATATC‑GCCATTCCCAAAGCAACCAACTAACAGCCTTTGACACACCGC‑ATGATCGGCC‑CGccaggccgcccgggt                                                                                                                                                                                                                                                        >  2:107906/1‑284 (MQ=255)
                             aaTCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAACTGGTGCTCAAAGAACAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTTCTGGGCCAGATGCACAGCCTGTACACACCGCGTCCGGCCCGCGTGAAAGACTATATC‑GCACTTCCAAAACCGAAAGGCCATCACGCTTTG‑CCCAACTCGCTGATTGACC‑CCCACGGCGGGCTGGta                                                                                                                                                                                                                                                        >  2:301448/1‑297 (MQ=255)
                                 gAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGCTTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTCCAAGCCGCGTCCGGGCCGCGTGACCGACTATATC‑GCCCTTCCAAAAGCGCACCACTACCCGTATTTG‑CACCCCTCCAGCATAGGCCACGC‑CCGCGTCCCGGTTGCgg                                                                                                                                                                                                                                                    >  2:295679/1‑298 (MQ=255)
                                                                 ccGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCCGCTGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAACGCCTATATC‑GCCATTCCAAACGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CCCCCCGAGGGCCGGTTGCGGGGCCGAACCGTCCCGAAGCTGTGGAGCAGagg                                                                                                                                                                                                                    >  2:186331/1‑296 (MQ=255)
                                                                           ggTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATgg                                                                                                                                                                                                          >  1:30629/1‑298 (MQ=255)
                                                                                 gAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTCTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTg                                                                                                                                                                                                    >  1:176991/1‑298 (MQ=255)
                                                                                  aaaaGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGc                                                                                                                                                                                                   <  1:321967/298‑1 (MQ=255)
                                                                                                                                  gTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGcaca                                                                                                                                                   >  1:274620/1‑298 (MQ=255)
                                                                                                                                                                                  aTTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACagag                                                                                                   >  1:280006/1‑298 (MQ=255)
                                                                                                                                                                                        acacCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCgg                                                                                             >  1:186139/1‑298 (MQ=255)
                                                                                                                                                                                          cccTGTTAGCAGGGGCGGGGGCAGGAGCACAGCCTGTACAAGGAGCGTCCGGGCCGGGTGAAGGACTAGATC‑GCCAGTGCAAAAGCGAACGAGTATCAGTCGTTG‑GACGGCTCGATGATGGGCC‑CGCACGGCGTGCCGGTTGTGTTCCAGGTCCGTAGCGAGGATATGGATGAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATCAGGAGCACGGCGAAACAAGTACTACCGCACAAATCCGCGCGCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTa                                                                                           <  2:15186/297‑1 (MQ=255)
                                                                                                                                                                                          acCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTa                                                                                           <  1:283682/298‑1 (MQ=255)
                                                                                                                                                                                                       tGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAAtt                                                                              >  1:15258/1‑298 (MQ=255)
                                                                                                                                                                                                                aGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGagag                                                                     <  1:100436/298‑1 (MQ=255)
                                                                                                                                                                                                                       caGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGGGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCCAAGCCTGCCGGAACTGCAACACAGCGCCGGTGGTTCGGTTGGATTTCTCGAGAGGGTTaaa                                                              >  2:102009/1‑298 (MQ=255)
                                                                                                                                                                                                                                     cgcgTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGCAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCCAAACCAGTACTACCGCACAAAGCCGCGCCCAGCGCTGGATGCAAACCCTGCTGGAGCCCCCAACGAGCGCCGGCAGTCCGCTTGAAGTTATAGAAGGCGTTAAATCCCATCTCTTccc                                                >  2:58789/1‑298 (MQ=255)
                                                                                                                                                                                                                                          ccGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGGTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCCAAACCAGTACTACCGCCCAAATCCGACCCCAGCCCCGGAAGGCAAGCCTTCTGGAACTGCAACAGAGCGCCGGTACTTCGTTTAACTTTATTGGACTCGTTACATCCGATCTCTTCCCgtgag                                           >  2:141512/1‑294 (MQ=255)
                                                                                                                                                                                                                                                 gcgTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTa                                   <  1:319167/298‑1 (MQ=255)
                                                                                                                                                                                                                                                 gcgTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTa                                   <  1:171459/298‑1 (MQ=255)
                                                                                                                                                                                                                                                      aaaGACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGttt                              <  1:295679/298‑1 (MQ=255)
                                                                                                                                                                                                                                                         gACTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGACGCTATGGATCAGATTCCGGCGATGGGTGCTGCCGCGCACTGGGCTTCTAACGACCACGGCCAAACCAATACTACCGCCCAAATCCGCACCCCGCGCTGGATGCCAAGCCTGCTGGCGCTGCACCATCGCGCCGGTAGTCCCTTCACCTTTACTGAGAACGTTAACTCCGACTGCTTCCCGGACGA‑GACTTACGTTTTca                           >  2:195240/1‑298 (MQ=255)
                                                                                                                                                                                                                                                           cTATATC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGTTTTcaca                         <  1:107906/298‑1 (MQ=255)
                                                                                                                                                                                                                                                            tataTC‑GCCATTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGTTGA‑GATTTACGTTTTCacac                        >  1:274661/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                      aTTCCAAAAGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGATGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCCAAAATCCGCGCCCCGCGCTGAATGCAAAGCCTGCCGGAGCTCCAACCGAGCGCGCGGAGTTCTTTCGACTCTACCGCGAGCGTCAAACCCGACCCCTTCCCCG‑TGAGGACTTACGTTTGCACCCCGGacgcgc               >  2:147127/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                           aaaaGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGTTTTCACACCGGAAGGGCGCAtt          <  1:186331/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                              aGCGAACGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGACGATATGGATCAGATGGCGGCGATGGGTGTTGCCGCGCACTGGGCTTCTCAACACCACCGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCCGGCTGCAAAGCCACCTGGAGCCGCACCCGCCCGCCGCTACTCCGCTCCCACTAATCAAGATCGTTAACTCCGCTCCCCCCCCGCATGA‑CATTCACGCGTTCACACCGCACGGGCCCaccttc       >  2:288952/1‑293 (MQ=255)
                                                                                                                                                                                                                                                                                  aaCGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCAAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGc   >  1:271138/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                   aCGACTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGct  <  1:301448/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                                                                         
AAAGAGATGATCCAAAAAATCCTCTCTGAAATCGAAGGGCGTTTGCAGGAAGCGGGAATACCGTGCCGCGTCAGTGGTCGCGAAAAGCATCTTTATTCGATTTACTGCAAAATGGTGCTCAAAGAGCAGCGTTTTCACTCAATCATGGACATCTACGCTTTCCGCGTGATCGTCAATGATTCTGACACCTGTTATCGCGTGCTGGGCCAGATGCACAGCCTGTACAAGCCGCGTCCGGGCCGCGTGAAAGACTATATC‑GCCATTCCAAAAGCGAACGGCTATCAGTCTTTG‑CACACCTCGATGATTGGCC‑CGCACGGCGTGCCGGTTGAGGTCCAGATCCGTACCGAAGATATGGATCAGATGGCGGAGATGGGTGTTGCCGCGCACTGGGCTTATAAAGAGCACGGCGAAACCAGTACTACCGCACAAATCCGCGCCCAGCGCTGGATGCAAAGCCTGCTGGAGCTGCAACAGAGCGCCGGTAGTTCGTTTGAATTTATCGAGAGCGTTAAATCCGATCTCTTCCCGGATGA‑GATTTACGTTTTCACACCGGAAGGGCGCATTGTCGAGCT  >  CP006584/3942985‑3943556

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