Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,370,777 C→T 13.0% M1M (ATG→ATA) † sspA ← stringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,370,7770CT13.0% 183.1 / 8.4 73M1M (ATG→ATA) †sspAstringent starvation protein A, phage P1 late gene activator, RNAP‑associated acid‑resistance protein, inactive glutathione S‑transferase homolog
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base T (6/4);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCAATACCCGACATTTGGGCAGAAAATTG  >  CP009273/3370631‑3370925
                                                                                                                                                  |                                                                                                                                                    
aGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACTTGATGTCTATGGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCAt                                                                                                                                                    <  1:544978/149‑1 (MQ=255)
   tCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCCCAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTCCCGACCGTTTGTTGGCAGCGGCAGCTATGaa                                                                                                                                                 >  1:460997/1‑149 (MQ=255)
   tCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGaa                                                                                                                                                 >  1:1015535/1‑149 (MQ=255)
   tCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGaa                                                                                                                                                 <  1:102865/149‑1 (MQ=255)
   tCCTGAGGCGCATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGaa                                                                                                                                                 <  1:1082800/149‑1 (MQ=255)
   acgtGAGGCGGACTGTCCTTTTCCACGTGGTCGATCTAGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTTATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGaa                                                                                                                                                 <  1:403988/146‑1 (MQ=255)
    ccTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGaaa                                                                                                                                                >  1:318916/1‑149 (MQ=255)
      tGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGAAAAc                                                                                                                                              <  1:197432/149‑1 (MQ=255)
         ggCGGATTGTCCTTTTCCACGTGTTCGATCTCGCAACTTCCACCGTTCTTAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTc                                                                                                                                           <  2:1182956/149‑1 (MQ=255)
         ggCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGCTGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTc                                                                                                                                           <  2:815180/149‑1 (MQ=255)
          gCGGATTGCCCTGTTCCACGTCTTCCATCTCGAAACTCACACCTTCCTCCGCCAGCACAATGCGGACCGGATGGCTATAGATGTCATTAGAACCGTAAATCAACGTCATTACCGACCGTTTGTTGGCAGCGACAGCCATGAAAACCTcc                                                                                                                                          <  2:934111/149‑1 (MQ=255)
           cGGATTGTCCTTTTACTCATGTTCCACCCCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCa                                                                                                                                         <  1:532358/149‑1 (MQ=255)
           cGGATTGTCCGTTTCCACGTGTGCGCTCTCGAAACTTTCACCTTTCTAAGCGCACAAATGGTGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCa                                                                                                                                         <  1:287074/149‑1 (MQ=255)
            ggATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTTGCAGCGACAGCTATGAAAACCTCCAg                                                                                                                                        >  1:498905/1‑149 (MQ=255)
            ggATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGAAAACCTCCAg                                                                                                                                        >  1:669679/1‑149 (MQ=255)
             gATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCGGATGGCTATAGATGTCAGTAGGACCGGAAACCAGCGTCATTACCGAACGTATGTGGGCAGCGACAGCTATGGAAACCTCCAgt                                                                                                                                       >  1:115806/1‑148 (MQ=255)
                 gTCCTTTTCCACGTGTTCGATCTCGAAACTTTCACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGATAGCCATGAAAACCTCCAGGtata                                                                                                                                   <  2:846140/149‑1 (MQ=255)
                  tCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCGTTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGAATAg                                                                                                                                  >  2:812335/1‑149 (MQ=255)
                  tCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAg                                                                                                                                  >  2:547792/1‑149 (MQ=255)
                     ttttCCACGTGTTCGATTTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCa                                                                                                                               >  2:50242/1‑149 (MQ=255)
                     ttttCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGAAAACCTCCAGGTATAGTCa                                                                                                                               >  1:264612/1‑149 (MQ=255)
                     ttttCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCa                                                                                                                               >  1:1200769/1‑149 (MQ=255)
                     ttttCCACGTGTTCGATCTCGAAACTTACACCATTCTCAGCCAGCCCAGTGAGCACCTGTTGGACAAAGAAATCACTAGCCCCGGAAAACTGCGACATTCCCCAAGATTGTTTTTCAGCTGTTGCCATGAAAACCACCAgaactagcaa                                                                                                                               >  2:1032064/1‑140 (MQ=255)
                       ttCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATCGATGTCGGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTCTAGTCAGa                                                                                                                             >  1:1094807/1‑149 (MQ=255)
                       ttCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGa                                                                                                                             >  1:447063/1‑149 (MQ=255)
                       ttCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGa                                                                                                                             >  1:765893/1‑149 (MQ=255)
                           tcGTGTTCGATTCCGTACCTTAGCCCTTTCTCAGCCAGCACAATGGGGACCTGATGGCTATAGATTTCAGTAAGACCGGAAAACAGCGTCATTACCGAACGTTTTTTGGCAGCGACAGCCGTGAAAACCTCCAGTTATAGTCAGAAttt                                                                                                                         <  2:699834/148‑1 (MQ=255)
                            cGTGTTCGATCTCGTAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCGTGAAAACCTCCAGGTATAGTCAGAAtttt                                                                                                                        <  1:777887/149‑1 (MQ=255)
                                       tcGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCa                                                                                                             <  1:718194/149‑1 (MQ=255)
                                        cGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAg                                                                                                            >  2:479041/1‑149 (MQ=255)
                                         gAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCGTTGCCGAACGTAGGTTGGCAGCGACAGCCATGAAAGCCGCCAGGGCTCGTCAGAATTTTTTTGGCTACCAGc                                                                                                           >  1:38408/1‑149 (MQ=255)
                                           aaCTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGcca                                                                                                         <  2:782029/149‑1 (MQ=255)
                                            aCTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCTTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGccac                                                                                                        <  1:793479/149‑1 (MQ=255)
                                            aCTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTACAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTGGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGccac                                                                                                        <  1:793710/149‑1 (MQ=255)
                                                 gcccTTTCTCAGCGAGCAGACTGCGGCCCTGGTGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGGTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGt                                                                                                   <  1:613436/147‑1 (MQ=255)
                                                  gccctactcAGCCAGCCCCATGAGGACCTGAGAGCTATAGATGTCAGTAGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGCATAGTCAGAATTTGTACTGCTACCAGCCACCAGGTgg                                                                                                 <  2:842720/143‑1 (MQ=255)
                                                   cgctccgcgggcgtcacaATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTTCTACCAGCCACCAGGTgg                                                                                                 <  2:798521/135‑1 (MQ=255)
                                                    cTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACGGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGGCAGAATTTTTACTCCTACCAGCCACCAGGTGGc                                                                                                >  2:363263/1‑149 (MQ=255)
                                                    cTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGc                                                                                                >  1:743814/1‑149 (MQ=255)
                                                      ttCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAGCGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCa                                                                                              <  1:208908/149‑1 (MQ=255)
                                                          cagccatcACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCa                                                                                          <  1:168756/149‑1 (MQ=255)
                                                             ccagcACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGaa                                                                                       <  1:191459/149‑1 (MQ=255)
                                                             ccagcACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGaa                                                                                       <  1:394283/149‑1 (MQ=255)
                                                                       gCGGACCTGATGGCTATCGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTAc                                                                             <  2:394685/149‑1 (MQ=255)
                                                                           aCCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCaa                                                                         <  2:351541/149‑1 (MQ=255)
                                                                           aCCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCaa                                                                         <  2:351382/149‑1 (MQ=255)
                                                                              tGATAGATATAGATGTCAGTAGGAACATAAAAACGAGTCATTAGCGCACGTTTGTTGGCAGCGACAACCATGGAAACCTCCAGGTGTAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCTAGTCAGAAGTTGCTTTACCCAGTaa                                                                      >  2:532105/1‑149 (MQ=255)
                                                                                aTGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAgg                                                                    <  1:271746/149‑1 (MQ=255)
                                                                                 tGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGACCGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGa                                                                   <  2:57962/149‑1 (MQ=255)
                                                                                  ggCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGaa                                                                  >  1:919888/1‑149 (MQ=255)
                                                                                      ataGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGAc                                                              >  1:154251/1‑149 (MQ=255)
                                                                                           ttgcggTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGAAAACCTCCAGGTATATTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGActctct                                                         <  2:498905/144‑1 (MQ=255)
                                                                                           tGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGActctct                                                         <  2:669679/149‑1 (MQ=255)
                                                                                                      ccGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGaaaa                                              <  2:318916/149‑1 (MQ=255)
                                                                                                      ccGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGaaaa                                              <  1:479041/149‑1 (MQ=255)
                                                                                                         gAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCTATGAAAACCTCCAGGTATAGTCAGAATTTTTCCTGCTACCAGCCACCCGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATca                                           <  2:115806/149‑1 (MQ=255)
                                                                                                         gAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAAGTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATca                                           <  1:794817/149‑1 (MQ=255)
                                                                                                          aaaaCAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTCTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATcaa                                          >  1:878264/1‑149 (MQ=255)
                                                                                                           aaaCAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATcaaa                                         <  1:547792/149‑1 (MQ=255)
                                                                                                           aaaCAGCGTCATCACCGAACGTTTGTTGGCCGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATcaaa                                         <  1:812335/149‑1 (MQ=255)
                                                                                                                 cGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGCTGTAGTCAGAATGTTTACTGCTACCAGCCACAAGGTGGCCGGTCAGAAGTTGTGTTACACCATAAGGAACGACTCTCTTTGTTCGAAAATcaaaacaaaa                                   >  2:334785/1‑149 (MQ=255)
                                                                                                                    cATTACCGAACGTGTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCCCTTTGTTCGAAAATCAAACAAAAAATg                                >  2:753357/1‑149 (MQ=255)
                                                                                                                     aTTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCACACAAAAAATGa                               >  2:389803/1‑149 (MQ=255)
                                                                                                                     aTTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGa                               >  1:1197832/1‑149 (MQ=255)
                                                                                                                     aTTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGa                               >  1:771086/1‑149 (MQ=255)
                                                                                                                         ccGAACGTGTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCaa                           >  2:699252/1‑149 (MQ=255)
                                                                                                                            aaCGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCACTAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCAATAc                        >  1:111960/1‑149 (MQ=255)
                                                                                                                            aaCGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGCGCaaaat                        >  2:923472/1‑146 (MQ=255)
                                                                                                                             aCGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGCACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCAATAcc                       >  2:800234/1‑149 (MQ=255)
                                                                                                                                tttGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAAGGAGCAAAACCCGa                    >  1:635309/1‑149 (MQ=255)
                                                                                                                                       gCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCACCAAGGCACGCCCCCCTTTGTCCGACAATCAAACAAAAACTGAGCACTACCCGACATTTgg             >  2:404956/1‑149 (MQ=255)
                                                                                                                                       gCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCTGACGAAAGATGGGCAATACCCGGCGTTTgg             >  1:1030273/1‑149 (MQ=255)
                                                                                                                                       gCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCAATACCCGACATTTgg             <  1:50242/149‑1 (MQ=255)
                                                                                                                                                  cATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCAATACCCGACATTTGGGCAGAAAATTg  >  2:1088109/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AGATCCTGAGGCGGATTGTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTGTGTTACCCAATAAGGAACGACTCTCTTTGTTCGAAAATCAAACAAAAAATGAGCAATACCCGACATTTGGGCAGAAAATTG  >  CP009273/3370631‑3370925

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