Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,330,419:1 +T 7.7% 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 7.7% 432.8 / 31.6 104intergenic (+23/+163)dacB/obgED‑alanyl‑D‑alanine carboxypeptidase/GTPase involved in cell partioning and DNA repair
Reads supporting (aligned to +/- strand):  ref base . (50/46);  new base T (4/4);  total (54/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.86e-01

Only 100 of 105 total aligned reads displayed.

CAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC  >  NC_000913/3330271‑3330566
                                                                                                                                                     |                                                                                                                                                   
cAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATg                                                                                                                                                      >  2:204767/1‑149 (MQ=255)
     aaTGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCgg                                                                                                                                                >  3:677732/1‑149 (MQ=255)
       tGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGC‑CATGTGCCgggg                                                                                                                                              >  4:107728/1‑149 (MQ=255)
           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:204767/149‑1 (MQ=255)
                   aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattatt                                                                                                                                  <  1:552881/149‑1 (MQ=255)
                   aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                                                                                                                   <  2:218398/148‑1 (MQ=255)
                   aaTATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGattat                                                                                                                                   >  1:218398/1‑148 (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:224213/149‑1 (MQ=255)
                             tGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgc                                                                                                                        >  3:247618/1‑149 (MQ=255)
                               gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                                                                      >  4:737922/1‑149 (MQ=255)
                               gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTgcgc                                                                                                                      >  2:433657/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:793686/1‑149 (MQ=255)
                                                                gCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTa                                                                                     >  4:57120/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                                                                                           >  3:655994/1‑119 (MQ=255)
                                                                                        aaaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAAc                                                                                           <  4:655994/119‑1 (MQ=255)
                                                                                         aaGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAAtt                                                            >  2:918152/1‑149 (MQ=255)
                                                                                            ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                         <  2:249469/149‑1 (MQ=255)
                                                                                            ccGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCaa                                                         <  4:293313/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:57162/149‑1 (MQ=255)
                                                                                                              ttATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTg                                       <  3:123025/149‑1 (MQ=255)
                                                                                                                aTCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTa                                     <  1:232697/149‑1 (MQ=255)
                                                                                                                    gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                                 <  1:59728/149‑1 (MQ=255)
                                                                                                                    gAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGcc                                 >  4:89857/1‑149 (MQ=255)
                                                                                                                        aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTa                                                                                                              <  1:651452/68‑1 (MQ=255)
                                                                                                                        aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                                            >  2:163936/1‑86 (MQ=255)
                                                                                                                        aaTTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTaaa                                                                                            <  1:163936/86‑1 (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                                                                                                                 <  4:840723/57‑1 (MQ=255)
                                                                                                                                aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                                                                                                                 >  3:840723/1‑57 (MQ=255)
                                                                                                                                aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  3:579448/1‑149 (MQ=255)
                                                                                                                                aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  1:204388/1‑149 (MQ=255)
                                                                                                                                aaaaaGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAg                     >  4:424036/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                >  1:435716/1‑149 (MQ=255)
                                                                                                                                     gAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCAt                >  3:885465/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)
                                                                                                                                                     |                                                                                                                                                   
CAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATG‑GCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAACATATTGCAATTCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCTGATTTGGC  >  NC_000913/3330271‑3330566

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: ‑