Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,576,245 T→A pseudogene (2111/2721 nt) yjiV → putative uncharacterized protein YjiV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,576,2450TA97.0% 224.3 / ‑0.8 67pseudogene (2111/2721 nt)yjiVputative uncharacterized protein YjiV
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (28/37);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

AACGATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGATCCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCGTTCACTATACCGATCGCTAT  >  NC_000913/4576108‑4576378
                                                                                                                                         |                                                                                                                                     
aaCGATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAAcc                                                                                                                                     <  2:748297/140‑1 (MQ=255)
 aCGATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACca                                                                                                                                    >  2:523951/1‑140 (MQ=255)
   gATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGACCGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACcaca                                                                                                                                  <  2:762165/140‑1 (MQ=255)
   gATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACcaca                                                                                                                                  <  1:434647/140‑1 (MQ=255)
     ttCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACaaa                                                                                                                                <  2:47176/140‑1 (MQ=255)
        cGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAg                                                                                                                             <  2:175020/140‑1 (MQ=255)
        cGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAg                                                                                                                             >  2:625618/1‑140 (MQ=255)
        cGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAg                                                                                                                             >  2:663198/1‑140 (MQ=255)
               ggAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATg                                                                                                                      >  1:623328/1‑140 (MQ=255)
                     gCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCAc                                                                                                                >  2:544168/1‑140 (MQ=255)
                      cACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACt                                                                                                               >  2:328424/1‑140 (MQ=255)
                       aCCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTg                                                                                                              >  2:845591/1‑140 (MQ=255)
                        ccTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGt                                                                                                             >  2:660109/1‑140 (MQ=255)
                         cTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTc                                                                                                            >  2:844576/1‑140 (MQ=255)
                           gAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGc                                                                                                          >  1:204946/1‑140 (MQ=255)
                           gAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTAGc                                                                                                          >  1:672039/1‑140 (MQ=255)
                                      tggtggTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               <  1:523951/140‑1 (MQ=255)
                                      tggtggTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               >  1:429195/1‑140 (MQ=255)
                                      tggtggTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               <  2:693615/140‑1 (MQ=255)
                                      tggtggTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               <  1:267145/140‑1 (MQ=255)
                                      tggtggTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               <  1:663198/140‑1 (MQ=255)
                                      tggtggTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               >  2:267145/1‑140 (MQ=255)
                                       ggtggtGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAAc                                                                                              <  1:328424/140‑1 (MQ=255)
                                        gtggtgCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               >  1:477841/1‑138 (MQ=255)
                                        gtggtgCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               <  2:477841/138‑1 (MQ=255)
                                        gtggtgCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACg                                                                                             >  1:70677/1‑140 (MQ=255)
                                        gtggtgCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACg                                                                                             <  2:373346/140‑1 (MQ=255)
                                        gtggtgCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGACAACg                                                                                             <  1:163984/140‑1 (MQ=255)
                                            tgCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTc                                                                                         <  1:844576/140‑1 (MQ=255)
                                              cGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               <  2:188409/132‑1 (MQ=255)
                                              cGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCaa                                                                                               >  1:188409/1‑132 (MQ=255)
                                                  gTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGAc                                                                                   >  1:28906/1‑140 (MQ=255)
                                                  gTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGAc                                                                                   >  1:767661/1‑140 (MQ=255)
                                                    ctcGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGt                                                                                 <  2:204946/140‑1 (MQ=255)
                                                     tcGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTg                                                                                <  1:660109/140‑1 (MQ=255)
                                                      cGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGc                                                                               <  2:306337/140‑1 (MQ=255)
                                                        gTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCtt                                                                                          <  1:146084/127‑1 (MQ=255)
                                                        gTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCtt                                                                                          >  2:146084/1‑127 (MQ=255)
                                                             aaaaCGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTaa                                                                        <  2:429195/140‑1 (MQ=255)
                                                              aaaCGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAAc                                                                       <  2:740589/140‑1 (MQ=255)
                                                                  ggTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGAcc                                                                   <  1:625618/140‑1 (MQ=255)
                                                                      gAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATga                                                               >  2:369971/1‑140 (MQ=255)
                                                                        aCTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAGTGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATgaag                                                             <  2:789224/140‑1 (MQ=255)
                                                                        aCTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATgaag                                                             <  2:623328/140‑1 (MQ=255)
                                                                            aaCGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAgc                                                         <  2:728046/140‑1 (MQ=255)
                                                                                aGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTAACGACCATGAAGAAGCGCAGt                                                    <  2:28906/140‑1 (MQ=255)
                                                                                  tttATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTc                                                   >  2:799362/1‑140 (MQ=255)
                                                                                   ttATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTAGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCa                                                  <  2:672039/140‑1 (MQ=255)
                                                                                    tATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCAc                                                 <  1:369971/140‑1 (MQ=255)
                                                                                           ccGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTTAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGaa                                          <  2:717645/140‑1 (MQ=255)
                                                                                           ccGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGaa                                          >  2:763125/1‑140 (MQ=255)
                                                                                           ccGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGaa                                          >  2:761049/1‑140 (MQ=255)
                                                                                           ccGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGaa                                          >  1:327112/1‑140 (MQ=255)
                                                                                              gCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATaa                                       <  2:599421/140‑1 (MQ=255)
                                                                                               cTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAac                                      <  2:767661/140‑1 (MQ=255)
                                                                                                  aaGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACAcg                                   <  2:323862/140‑1 (MQ=255)
                                                                                                   aGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGATCCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACAcgc                                  >  2:441921/1‑140 (MQ=255)
                                                                                                        ccaACGCTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAAtt                             >  1:827716/1‑140 (MQ=255)
                                                                                                              cTGTCGAGCGGGTCAAAGACATTCAGATCCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGc                       <  1:441921/140‑1 (MQ=255)
                                                                                                              cTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGc                       <  1:47372/140‑1 (MQ=255)
                                                                                                              cTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGc                       >  1:724953/1‑140 (MQ=255)
                                                                                                              cTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGc                       <  1:761049/140‑1 (MQ=255)
                                                                                                              cTGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGc                       <  2:327112/140‑1 (MQ=255)
                                                                                                               tGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCg                      >  2:58255/1‑140 (MQ=255)
                                                                                                               tGTCGAGCGGGTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCg                      <  2:740618/140‑1 (MQ=255)
                                                                                                                        ggTCAAAGACATTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCGTTCACtata             <  1:845591/140‑1 (MQ=255)
                                                                                                                                  aTTCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCGTTCACTATACCGATCGCTa   <  1:58255/140‑1 (MQ=255)
                                                                                                                                   ttCAGAACCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCGTTCACTATACCGATCGCTAt  >  2:41696/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                     
AACGATTCCGGGGCCGGAGTGGCACCTGAATGACGAACTGGTGGTGCGTAGTCTCGGTTATAAAACGGTAGAACTGAACGAGTTTATTCTGCCGGCTAAAGCCACCAACGCTGTCGAGCGGGTCAAAGACATTCAGATCCACAAACAGTTGAATGGACCACTGTCGCAGTTTGGGCAACGCTTCTGGGACGTGCTGTTTAACGACCATGAAGAAGCGCAGTCACTGATGAATAACACGCGAATTACCGGCGTTCACTATACCGATCGCTAT  >  NC_000913/4576108‑4576378

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