Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,419:1 +T 5.2% intergenic (+23/+163) dacB → / ← obgE D‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,330,4191.T 5.2% 813.8 / 40.2 191intergenic (+23/+163)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base . (96/85);  new base T (4/6);  total (100/91)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-02

Only 100 of 195 total aligned reads displayed.

AACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCG  >  NC_000913/3330272‑3330567
                                                                                                                                                    |                                                                                                                                                    
aaCGAATGGCGTTTGTGCAATATCTTTCTGGGTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATgg                                                                                                                                                     >  4:1059912/1‑149 (MQ=255)
aaCGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATgg                                                                                                                                                     >  2:966960/1‑149 (MQ=255)
    aaTGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCgg                                                                                                                                                 <  3:991984/149‑1 (MQ=255)
       ggCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGc                                                                                                                                              <  4:250702/149‑1 (MQ=255)
               tgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGatt                                                                                                                                      <  2:1131873/149‑1 (MQ=255)
               tgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGatt                                                                                                                                      >  2:445632/1‑149 (MQ=255)
               tgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCCATGTGCCGGGGTTTcttttt                                                                                                                                      <  4:629648/149‑6 (MQ=255)
                gCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGatta                                                                                                                                     <  2:632489/149‑1 (MQ=255)
                 cAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattat                                                                                                                                    <  1:793459/149‑1 (MQ=255)
                 cAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattat                                                                                                                                    <  1:1660225/149‑1 (MQ=255)
                 cAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCCATGTGCCGGGGTTTctttttga                                                                                                                                    <  4:978958/149‑8 (MQ=255)
                  aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   >  1:1380822/1‑149 (MQ=255)
                  aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   >  1:1505737/1‑149 (MQ=255)
                  aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   >  2:185213/1‑149 (MQ=255)
                  aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   <  4:719919/149‑1 (MQ=255)
                  aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   >  2:1040152/1‑149 (MQ=255)
                    tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGa                                                                                                                                 >  2:1606067/1‑149 (MQ=255)
                     atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCCATGTGCCgggg                                                                                                                                               <  1:1595478/134‑1 (MQ=255)
                     atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCCATGTGCCgggg                                                                                                                                               >  2:1595478/1‑134 (MQ=255)
                        tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTGGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             <  3:1243612/149‑1 (MQ=255)
                        tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   >  2:106403/1‑143 (MQ=255)
                        tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGattatt                                                                                                                                   <  1:106403/143‑1 (MQ=255)
                        tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             >  1:533737/1‑149 (MQ=255)
                        tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             >  4:232907/1‑149 (MQ=255)
                        tttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGCGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGGCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             >  1:81676/1‑149 (MQ=255)
                        tttCTGGCTAGGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             <  4:94012/149‑1 (MQ=255)
                              gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                                                                       >  1:64942/1‑149 (MQ=255)
                                tATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCtt                                                                                                                     >  3:1718191/1‑149 (MQ=255)
                                tATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTGTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCtt                                                                                                                     >  1:431251/1‑149 (MQ=255)
                                 aTGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCttt                                                                                                                    <  1:1114270/149‑1 (MQ=255)
                                   gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                  >  1:203098/1‑149 (MQ=255)
                                   gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                  <  2:250387/149‑1 (MQ=255)
                                   gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                  >  4:112079/1‑149 (MQ=255)
                                   gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                  <  1:275254/149‑1 (MQ=255)
                                    cAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                                                                                 <  1:1206195/149‑1 (MQ=255)
                                    cAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                                                                                 <  3:927373/149‑1 (MQ=255)
                                     aGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCa                                                                                                                                          >  4:305228/1‑123 (MQ=255)
                                     aGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGGTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGttt                                                                                                                >  2:1702402/1‑149 (MQ=255)
                                      gTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                               <  4:1090637/149‑1 (MQ=255)
                                       tAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTAt                                                                                                              >  3:1116330/1‑149 (MQ=255)
                                       tAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTAt                                                                                                              <  2:334153/149‑1 (MQ=255)
                                             cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGttt                                                                                                                >  2:1290316/1‑141 (MQ=255)
                                             cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGa                                                                                                        >  1:704006/1‑149 (MQ=255)
                                              tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGAt                                                                                                       >  2:529844/1‑149 (MQ=255)
                                                     cAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                >  2:1505674/1‑149 (MQ=255)
                                                       gCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaa                                                                                              >  2:880344/1‑149 (MQ=255)
                                                         gTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAAc                                                                                            <  2:877536/149‑1 (MQ=255)
                                                         gTAATCGCCGTATTCCGGTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGTCGTAAAc                                                                                            <  4:1459067/149‑1 (MQ=255)
                                                                ccGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATCAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                     >  3:215458/1‑149 (MQ=255)
                                                                    aTTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                                                                                 <  2:1368311/117‑1 (MQ=255)
                                                                    aTTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGtt                                                                                                                 >  1:1368311/1‑117 (MQ=255)
                                                                         gTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                            >  1:750586/1‑149 (MQ=255)
                                                                         gTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                            >  3:1069782/1‑149 (MQ=255)
                                                                               gCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                      <  1:346905/149‑1 (MQ=255)
                                                                               gCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                      <  2:459376/149‑1 (MQ=255)
                                                                                 gTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                    <  3:627614/149‑1 (MQ=255)
                                                                                      gAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                               >  1:491642/1‑149 (MQ=255)
                                                                                        aaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                             >  2:1302254/1‑149 (MQ=255)
                                                                                         aGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTc                                                            <  2:390622/149‑1 (MQ=255)
                                                                                          gCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                           >  2:1089972/1‑149 (MQ=255)
                                                                                           ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             <  1:1324751/82‑1 (MQ=255)
                                                                                           ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACaaa                                                                                                                             >  2:1324751/1‑82 (MQ=255)
                                                                                              tttGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                    <  4:48545/136‑1 (MQ=255)
                                                                                              tttGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGt                                                                    >  3:48545/1‑136 (MQ=255)
                                                                                              tttGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACat                                                       >  4:1658163/1‑149 (MQ=255)
                                                                                                  tataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTg                                                   >  2:240340/1‑149 (MQ=255)
                                                                                                   ataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                  >  3:733597/1‑149 (MQ=255)
                                                                                                    taAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCa                                                 >  3:1218247/1‑149 (MQ=255)
                                                                                                     aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCt                                                                                                                      >  2:1521262/1‑79 (MQ=255)
                                                                                                     aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCaa                                                <  4:1218247/149‑1 (MQ=255)
                                                                                                     aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCaa                                                >  1:1390721/1‑149 (MQ=255)
                                                                                                      aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                                               >  3:1262585/1‑149 (MQ=255)
                                                                                                          tatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                           <  3:273296/149‑1 (MQ=255)
                                                                                                            tttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                           <  4:825841/147‑1 (MQ=255)
                                                                                                            tttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                           >  3:825841/1‑147 (MQ=255)
                                                                                                               aTCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                      >  3:957169/1‑149 (MQ=255)
                                                                                                                  aGAACAATTAGTCAAAAAGAAGCCCCGGCACATGGCCGGGTCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGc                                   <  1:102092/149‑1 (MQ=255)
                                                                                                                   gAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                         <  4:514508/94‑1 (MQ=255)
                                                                                                                   gAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                         >  3:514508/1‑94 (MQ=255)
                                                                                                                    aaCAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTCTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAGTTCAACATATTGCAATTCTCTTGTAGGCCt                                 >  1:1136063/1‑149 (MQ=255)
                                                                                                                    aaCAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCt                                 <  2:1175132/149‑1 (MQ=255)
                                                                                                                      cAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGa                               <  2:710898/149‑1 (MQ=255)
                                                                                                                       aaTTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                              >  3:363401/1‑149 (MQ=255)
                                                                                                                       aaTTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                              <  3:1124679/149‑1 (MQ=255)
                                                                                                                       aaTTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                              <  4:527941/149‑1 (MQ=255)
                                                                                                                        aTTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATa                             >  4:1079806/1‑149 (MQ=255)
                                                                                                                             tCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                           >  1:1500698/1‑130 (MQ=255)
                                                                                                                             tCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtctc                                           <  2:1500698/130‑1 (MQ=255)
                                                                                                                             tCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGt                        >  4:413003/1‑149 (MQ=255)
                                                                                                                             tCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGt                        >  3:1576825/1‑149 (MQ=255)
                                                                                                                               aaaaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                                  >  1:924276/1‑137 (MQ=255)
                                                                                                                               aaaaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATaa                            <  4:1346897/143‑1 (MQ=255)
                                                                                                                                aaaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc                     >  1:1535076/1‑149 (MQ=255)
                                                                                                                                 aaaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGTGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcg                    >  3:46113/1‑149 (MQ=255)
                                                                                                                                  aaGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcgc                   >  1:1288290/1‑149 (MQ=255)
                                                                                                                                   tgAAGCCCCGGCCATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCa                  >  3:1044350/2‑149 (MQ=18)
                                                                                                                                     aagccccGGCCATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                              >  1:1298261/4‑135 (MQ=18)
                                                                                                                                         cacGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc            <  3:1395363/147‑1 (MQ=255)
                                                                                                                                                caTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTg     <  4:1558762/149‑1 (MQ=255)
                                                                                                                                                 aTGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTgg    >  3:1532421/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCG  >  NC_000913/3330272‑3330567

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