Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,427 C→T 5.7% intergenic (+31/+155) 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,4270CT 5.7% 342.7 / 4.6 102intergenic (+31/+155)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base C (49/47);  new base T (3/3);  total (52/50)
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.95e-01

Only 100 of 102 total aligned reads displayed.

GTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATC  >  NC_000913/3330282‑3330573
                                                                                                                                                  |                                                                                                                                                  
gTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCGGGGTTTc                                                                                                                                                 <  3:107728/149‑1 (MQ=255)
 tttGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCa                                                                                                                                                <  2:937265/149‑1 (MQ=255)
     tgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGatt                                                                                                                                            >  3:293313/1‑149 (MQ=255)
      gCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGatta                                                                                                                                           <  1:239705/149‑1 (MQ=255)
      gCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGatta                                                                                                                                           >  2:232697/1‑149 (MQ=255)
        aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                                         <  1:552881/149‑1 (MQ=255)
        aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                                         <  1:204767/149‑1 (MQ=255)
        aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                                         >  3:637406/1‑149 (MQ=255)
        aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                                                                                                                          >  1:218398/1‑148 (MQ=255)
        aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                                                                                                                          <  2:218398/148‑1 (MQ=255)
          tatCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGa                                                                                                                                       <  3:205880/149‑1 (MQ=255)
           atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGAc                                                                                                                                      >  3:475638/1‑149 (MQ=255)
                tCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGt                                                                                                                                 <  1:29762/149‑1 (MQ=255)
                  tGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgc                                                                                                                               >  3:247618/1‑149 (MQ=255)
                    gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                                                                             >  2:433657/1‑149 (MQ=255)
                    gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                                                                             >  4:737922/1‑149 (MQ=255)
                      tATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCtt                                                                                                                           >  2:157104/1‑149 (MQ=255)
                       aTGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCttt                                                                                                                          >  4:338191/1‑149 (MQ=255)
                        tGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTg                                                                                                                         <  1:369539/149‑1 (MQ=255)
                        tGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTg                                                                                                                         <  1:555530/149‑1 (MQ=255)
                         gCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                        <  1:243776/149‑1 (MQ=255)
                            gTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                                     >  2:332386/1‑149 (MQ=255)
                            gTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                                     >  1:110964/1‑149 (MQ=255)
                               gAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGc                                                                                                                  <  1:905084/149‑1 (MQ=255)
                                   cTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGa                                                                                                              <  1:332386/149‑1 (MQ=255)
                                     gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                            >  4:113158/1‑149 (MQ=255)
                                         aTCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        >  1:403728/1‑149 (MQ=255)
                                         aTCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  1:433657/149‑1 (MQ=255)
                                              cGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                                                   <  4:637406/149‑1 (MQ=255)
                                                  aTCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGcc                                                                                               <  1:448209/149‑1 (MQ=255)
                                                     gCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                                                                            >  4:57120/1‑149 (MQ=255)
                                                     gCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                                                                            >  4:793686/1‑149 (MQ=255)
                                                      ccGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt                                                                                           >  3:39434/1‑149 (MQ=255)
                                                        gTATTACGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAATCCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTGTcc                                                                                         <  4:233886/149‑1 (MQ=255)
                                                           ttCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                                      <  4:884647/149‑1 (MQ=255)
                                                              cGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                   <  4:261754/149‑1 (MQ=255)
                                                              cGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                                   >  3:791204/1‑149 (MQ=255)
                                                                        ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                         >  4:794187/1‑149 (MQ=255)
                                                                        ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                         >  1:308758/1‑149 (MQ=255)
                                                                        ttttGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTg                                                                         >  3:810779/1‑149 (MQ=255)
                                                                             aaaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAAc                                                                                                  <  4:655994/119‑1 (MQ=255)
                                                                             aaaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAAc                                                                                                  >  3:655994/1‑119 (MQ=255)
                                                                              aaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                                   >  2:918152/1‑149 (MQ=255)
                                                                                 ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                <  4:293313/149‑1 (MQ=255)
                                                                                 ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                <  2:249469/149‑1 (MQ=255)
                                                                                 ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                <  4:247618/149‑1 (MQ=255)
                                                                                 ccGTTTGTATAAAGATATTTATCAGAAAAATTAGTCAAAAAGAAACCCCTGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                <  3:794187/149‑1 (MQ=255)
                                                                                    tttGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACat                                                             >  1:762553/1‑149 (MQ=255)
                                                                                      tGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACatat                                                           <  1:354357/149‑1 (MQ=255)
                                                                                         ataAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGc                                                        >  4:129194/1‑149 (MQ=255)
                                                                                          taAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCa                                                       >  1:149662/1‑149 (MQ=255)
                                                                                           aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                <  3:378226/139‑1 (MQ=255)
                                                                                           aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                >  4:378226/1‑139 (MQ=255)
                                                                                           aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCaa                                                      <  1:23142/149‑1 (MQ=255)
                                                                                           aaaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCaa                                                      >  1:619219/1‑149 (MQ=255)
                                                                                            aaGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAt                                                     >  2:59728/1‑149 (MQ=255)
                                                                                             aGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACTAAATCGTGCAAATTCAACATATTGCAAtt                                                    <  3:793686/149‑1 (MQ=255)
                                                                                             aGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                                                    <  3:57120/149‑1 (MQ=255)
                                                                                             aGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAAtt                                                    >  3:851647/1‑149 (MQ=255)
                                                                                               atatTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTCTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATtct                                                  >  3:873213/1‑149 (MQ=255)
                                                                                                  tttcttaGAACAATTAGTCAAAAAGAACCCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                               <  4:798965/143‑1 (MQ=255)
                                                                                                  tttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                               <  4:851647/149‑1 (MQ=255)
                                                                                                   ttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                              <  3:123025/149‑1 (MQ=255)
                                                                                                   ttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                              <  3:57162/149‑1 (MQ=255)
                                                                                                     aTCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                            <  1:232697/149‑1 (MQ=255)
                                                                                                         gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                                        >  4:89857/1‑149 (MQ=255)
                                                                                                         gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                                        <  1:59728/149‑1 (MQ=255)
                                                                                                             aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                                     >  2:651452/1‑68 (MQ=255)
                                                                                                             aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                                     <  1:651452/68‑1 (MQ=255)
                                                                                                             aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                                                   <  1:163936/86‑1 (MQ=255)
                                                                                                             aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                                                   >  2:163936/1‑86 (MQ=255)
                                                                                                             aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGAt                                    >  2:92270/1‑149 (MQ=255)
                                                                                                                   tCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                            <  3:188113/103‑1 (MQ=255)
                                                                                                                   tCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATc                                                                            >  4:188113/1‑103 (MQ=255)
                                                                                                                   tCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                >  4:777777/1‑115 (MQ=255)
                                                                                                                   tCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                                <  3:777777/115‑1 (MQ=255)
                                                                                                                   ataatctgAAGCCCCGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGt                              >  4:419468/8‑149 (MQ=18)
                                                                                                                     aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                        >  3:840723/1‑57 (MQ=255)
                                                                                                                     aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                        <  4:840723/57‑1 (MQ=255)
                                                                                                                     aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                            >  4:424036/1‑149 (MQ=255)
                                                                                                                     aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                            >  1:204388/1‑149 (MQ=255)
                                                                                                                     aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                            >  3:579448/1‑149 (MQ=255)
                                                                                                                      atctgAAGCCCCGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc                           <  3:419468/145‑1 (MQ=18)
                                                                                                                      aaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgc                           >  2:412734/1‑149 (MQ=255)
                                                                                                                        ctgAAGCCCCGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAgcgc                         >  2:655378/3‑149 (MQ=18)
                                                                                                                          gAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                       >  3:885465/1‑149 (MQ=255)
                                                                                                                          gAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                       >  1:435716/1‑149 (MQ=255)
                                                                                                                           aagccccGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                               >  3:160949/4‑124 (MQ=17)
                                                                                                                           aagccccGGC‑CATGTGCCGGGGTTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCtt                                               <  4:160949/121‑1 (MQ=18)
                                                                                                                           aaaCCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATc                      <  1:92270/149‑1 (MQ=255)
                                                                                                                            aaCCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCa                     <  1:43486/149‑1 (MQ=255)
                                                                                                                               cccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGc                  >  3:833840/1‑149 (MQ=255)
                                                                                                                                ccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                 <  1:595222/149‑1 (MQ=255)
                                                                                                                                ccGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCt                 >  4:128/1‑149 (MQ=255)
                                                                                                                                     acaTG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTa                             >  3:890230/1‑132 (MQ=255)
                                                                                                                                     acaTG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTa                             <  4:890230/132‑1 (MQ=255)
                                                                                                                                     acaTG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGAttt            <  3:113158/149‑1 (MQ=255)
                                                                                                                                       aTG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTgg          >  2:438002/1‑149 (MQ=255)
                                                                                                                                        tG‑GCCGGGGCTTCAGATTATTCACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGc         <  4:873213/149‑1 (MQ=255)
                                                                                                                                                ggCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTatc  <  4:475638/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGCGTTTATC  >  NC_000913/3330282‑3330573

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