Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,306,862 C→A E97* (GAA→TAA)  yciU ← DUF440 domain‑containing protein YciU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,306,8620CA98.4% 217.5 / ‑2.8 62E97* (GAA→TAA) yciUDUF440 domain‑containing protein YciU
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (37/24);  total (38/24)
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.58e-01

TTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGAT  >  NC_000913/1306726‑1306996
                                                                                                                                        |                                                                                                                                          
tttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTagcg                                                                                                                                         <  1:1174746/140‑1 (MQ=255)
 ttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCgg                                                                                                                                        <  1:585306/140‑1 (MQ=255)
 ttCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCgg                                                                                                                                        <  1:1428260/140‑1 (MQ=255)
  tCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGc                                                                                                                                       <  2:409017/140‑1 (MQ=255)
          tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGa                                                                                                                               >  2:1388049/1‑140 (MQ=255)
          tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGa                                                                                                                               <  2:1270541/140‑1 (MQ=255)
          tCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGa                                                                                                                               >  1:1270541/1‑140 (MQ=255)
            cGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAAt                                                                                                                             >  2:1342534/1‑140 (MQ=255)
             gTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATg                                                                                                                            <  1:1388049/140‑1 (MQ=255)
                    gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCaa                                                                                                                     >  2:53005/1‑140 (MQ=255)
                    gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCaa                                                                                                                     >  2:1075689/1‑140 (MQ=255)
                    gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCaa                                                                                                                     >  2:852892/1‑140 (MQ=255)
                    gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCaa                                                                                                                     >  2:870234/1‑140 (MQ=255)
                    gCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGAAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCaa                                                                                                                     >  2:1222091/1‑140 (MQ=255)
                             tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGt                                                                                                            <  1:870234/140‑1 (MQ=255)
                             tCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGt                                                                                                            <  1:1222091/140‑1 (MQ=255)
                                    gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTc                                                                                                     <  1:1075689/140‑1 (MQ=255)
                                    gCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTc                                                                                                     <  1:366987/140‑1 (MQ=255)
                                      tttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCAc                                                                                                   <  2:255180/140‑1 (MQ=255)
                                      tttAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGACTTCAc                                                                                                   <  1:1356834/140‑1 (MQ=255)
                                            aaCAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCtt                                                                                             <  1:53005/140‑1 (MQ=255)
                                              cAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCa                                                                                           <  1:852892/140‑1 (MQ=255)
                                              cAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCa                                                                                           <  1:1342534/140‑1 (MQ=255)
                                                aTAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACt                                                                                         >  1:1259713/1‑140 (MQ=255)
                                                 tAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTg                                                                                        >  2:1068922/1‑140 (MQ=255)
                                                        tGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCa                                                                             <  1:1068922/140‑1 (MQ=255)
                                                         gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCa                                                                             <  1:1083732/139‑1 (MQ=255)
                                                         gACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCa                                                                             >  2:1083732/1‑139 (MQ=255)
                                                          aCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGa                                                                           >  1:1134632/1‑140 (MQ=255)
                                                          aCAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGa                                                                           <  1:608884/140‑1 (MQ=255)
                                                           cAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGAc                                                                          >  2:608513/1‑140 (MQ=255)
                                                           cAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGAc                                                                          <  2:142435/140‑1 (MQ=255)
                                                               tGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAAt                                                                      >  2:113566/1‑140 (MQ=255)
                                                                gTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATc                                                                     >  1:155473/1‑140 (MQ=255)
                                                                   aGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATcacc                                                                  >  1:283577/1‑140 (MQ=255)
                                                                    gCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATcacca                                                                 >  2:280326/1‑140 (MQ=255)
                                                                    gCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATcacca                                                                 >  1:305062/1‑140 (MQ=255)
                                                                          tttCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCgg                                                           >  2:465032/1‑140 (MQ=255)
                                                                                           ttatttatTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTc                                          <  2:283577/140‑1 (MQ=255)
                                                                                            tatttattCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCa                                         >  2:1391617/1‑140 (MQ=255)
                                                                                                     ccGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAAt                                <  1:280326/140‑1 (MQ=255)
                                                                                                     ccGCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAAt                                >  2:1366915/1‑140 (MQ=255)
                                                                                                       gCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGTGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCa                              >  1:253802/1‑140 (MQ=255)
                                                                                                       gCCAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCa                              >  1:818703/1‑140 (MQ=255)
                                                                                                        ccAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCaa                             >  1:299623/1‑140 (MQ=255)
                                                                                                        ccAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCaa                             >  1:298379/1‑140 (MQ=255)
                                                                                                         cAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAAc                            <  1:706745/140‑1 (MQ=255)
                                                                                                         cAGATGATATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAAc                            <  1:1366915/140‑1 (MQ=255)
                                                                                                            atgatATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATg                         >  2:1395660/1‑140 (MQ=255)
                                                                                                            atgatATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATg                         >  1:609969/1‑140 (MQ=255)
                                                                                                            atgatATGGCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATg                         >  2:1031815/1‑140 (MQ=255)
                                                                                                                   ggCAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGc                  >  2:1692301/1‑140 (MQ=255)
                                                                                                                     cAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCt                    >  1:582152/1‑140 (MQ=255)
                                                                                                                     cAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGTCCGCCCGCCAGACCAATCACCACTTCGGAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTg                   >  1:933363/1‑140 (MQ=255)
                                                                                                                     cAAAGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCa                >  1:52418/1‑140 (MQ=255)
                                                                                                                      aaaGTTTGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAg               >  1:1240234/1‑140 (MQ=255)
                                                                                                                          tttGTGATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCtt           >  1:6262/1‑140 (MQ=255)
                                                                                                                              tgATCTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTccgc       >  2:957856/1‑140 (MQ=255)
                                                                                                                                  cTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGa   <  2:6262/140‑1 (MQ=255)
                                                                                                                                  cTTTTTAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGa   <  1:465032/140‑1 (MQ=255)
                                                                                                                                   tttttCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGAt  >  1:1696373/1‑140 (MQ=255)
                                                                                                                                   tttttAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGAt  >  1:484907/1‑140 (MQ=255)
                                                                                                                                   tttttAGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGAt  >  2:977103/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                          
TTTCAGACACTCCGTTTCTTGCCTGAAAGTCGATCCGCTTTAAAAACAATAAGGGCTGACAGTTGTCAGCCCTTTTTCACGCTAAAAGCGATTATTTATTCCCGCCAGATGATATGGCAAAGTTTGTGATCTTTTTCGCGGCATAACAGAATGCGGGCAAAAACATCGTTGATTTCACCATCTTCACTGT‑‑‑‑CCGCCAGACCAATCACCACTTCGGCAAAAAAGTCCGGGTTCAAATCGAAATCAACATGTTCCTGCCAGTCTTCCGCCGGAT  >  NC_000913/1306726‑1306996

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: