Predicted mutation
evidence seq id position mutation annotation gene description
MC JC NC_000913 4,474,849 20 bp→GTATTTCCCGGCATCTTTATAGCGATAGCA intergenic (+24/+2) yjgL → / ← argI hypothetical protein/ornithine carbamoyltransferase 1

Missing coverage evidence...
   seq id start end size ←reads reads→ gene description
* * ÷ NC_000913 4474849 4474868 20 213 [0] [0] 245 yjgL/argI hypothetical protein/ornithine carbamoyltransferase 1

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 44748480 (0.000)153 (0.970)
+GTATTTCCCGGCATCTTTATAGCGATAGCA
76/124 0.5 100% intergenic (+23/+22) yjgL/argI hypothetical protein/ornithine carbamoyltransferase 1
?NC_000913 4474869 = 0 (0.000)intergenic (+44/+1) yjgL/argI hypothetical protein/ornithine carbamoyltransferase 1

Only 100 of 213 total aligned reads displayed.

AATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAA  >  NC_000913/4474757‑4474940
                                                                                            |                                                                                            
aaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGttt                                                                                              <  2:5853518/93‑1 (MQ=255)
   ccACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgta                                                                                           >  2:4352741‑M1/1‑90 (MQ=255)
    cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                                                          >  2:1787278‑M1/1‑89 (MQ=255)
    cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                                                          <  2:267727‑M1/93‑5 (MQ=255)
    cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                                                          <  2:3758268‑M1/93‑5 (MQ=255)
     aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                                                         <  2:3832858‑M1/93‑6 (MQ=255)
        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      <  1:5999544‑M1/93‑9 (MQ=255)
        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      <  2:620783‑M1/93‑9 (MQ=255)
        ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                                                      >  2:4801732‑M1/1‑85 (MQ=255)
         gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                                                     >  2:5327009‑M1/1‑84 (MQ=255)
               ttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatc                                                                               <  2:3711325‑M1/93‑16 (MQ=255)
                 aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                                                             <  1:356423‑M1/93‑18 (MQ=255)
                  cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                                                            >  2:914641‑M1/1‑75 (MQ=255)
                  cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                                                            >  2:3426649‑M1/1‑75 (MQ=255)
                   gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                                                           >  1:667289‑M1/1‑74 (MQ=255)
                     ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                                                         >  1:3266559‑M1/1‑72 (MQ=255)
                     ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                                                         >  1:5422106‑M1/1‑72 (MQ=255)
                     ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                                                         >  2:1186171‑M1/1‑72 (MQ=255)
                      taGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatag                                                                        >  2:3944705‑M1/1‑71 (MQ=255)
                       aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagc                                                                       <  1:390898‑M1/93‑24 (MQ=255)
                        gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                                                      >  1:5176017‑M1/1‑69 (MQ=255)
                          aaTAGAACACTTTGGTTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                                                    <  1:550735‑M1/93‑27 (MQ=255)
                          aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                                                    >  1:2342827‑M1/1‑67 (MQ=255)
                          aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                                                    >  2:3648206‑M1/1‑67 (MQ=255)
                          aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                                                    >  2:5531806‑M1/1‑67 (MQ=255)
                           aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                                                   >  2:5083436‑M1/1‑66 (MQ=255)
                           aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                                                   >  2:4169251‑M1/1‑66 (MQ=255)
                            tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                                                  >  1:3624455‑M1/1‑65 (MQ=255)
                            tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                                                  >  1:2996833‑M1/1‑65 (MQ=255)
                            tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                                                  >  2:4898934‑M1/1‑65 (MQ=255)
                            tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                                                  >  2:4252231‑M1/1‑65 (MQ=255)
                             aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                                                 <  1:742042‑M1/93‑30 (MQ=255)
                                acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                                              >  2:858132‑M1/1‑61 (MQ=255)
                                 cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccgcaatctttatagcgatagcaatt                                                             >  2:1703161‑M1/1‑60 (MQ=255)
                                  acTTTGGGTGGCAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                                                            >  1:3010662‑M1/1‑59 (MQ=255)
                                   cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                                                           <  1:4638897‑M1/93‑36 (MQ=255)
                                   cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                                                           >  1:1019782‑M1/1‑58 (MQ=255)
                                   cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                                                           >  2:2092363‑M1/1‑58 (MQ=255)
                                   cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                                                           >  1:924245‑M1/1‑58 (MQ=255)
                                   cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                                                           >  2:4079256‑M1/1‑58 (MQ=255)
                                   cTTTGGGAGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                                                           >  1:3787086‑M1/1‑58 (MQ=255)
                                     ttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt                                                         <  1:1776358‑M1/93‑38 (MQ=255)
                                     ttGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattt                                                         <  2:3095212‑M1/93‑38 (MQ=255)
                                       gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                                                       >  1:4597413‑M1/1‑54 (MQ=255)
                                       gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                                                       <  2:1768927‑M1/93‑40 (MQ=255)
                                          tGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactga                                                    <  2:3662633‑M1/93‑43 (MQ=255)
                                           ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                                                   >  1:1451664‑M1/1‑50 (MQ=255)
                                           ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                                                   >  1:4435566‑M1/1‑50 (MQ=255)
                                           ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttagagcgatagcaattatttactgag                                                   >  1:4301212‑M1/1‑50 (MQ=255)
                                            gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                                                  <  2:762748‑M1/93‑45 (MQ=255)
                                              agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgt                                                >  1:1622743‑M1/1‑47 (MQ=255)
                                                 aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc                                             >  2:2850927‑M1/1‑44 (MQ=255)
                                                 aTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgc                                             >  2:4216875‑M1/1‑44 (MQ=255)
                                                   aGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcga                                           >  2:432283‑M1/1‑42 (MQ=255)
                                                     aCCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacc                                         >  1:1890799‑M1/1‑40 (MQ=255)
                                                       cTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccat                                       <  2:3562327‑M1/93‑56 (MQ=255)
                                                        tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatc                                      >  2:4588762‑M1/1‑37 (MQ=255)
                                                         aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatca                                     >  1:3880359‑M1/1‑36 (MQ=255)
                                                         aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatca                                     <  2:1996891‑M1/93‑58 (MQ=255)
                                                         aTCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatca                                     <  2:219833‑M1/93‑58 (MQ=255)
                                                          tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac                                    >  1:5931729‑M1/1‑35 (MQ=255)
                                                             cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc                                 >  1:2966180‑M1/1‑32 (MQ=255)
                                                             cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc                                 <  2:5683864‑M1/93‑62 (MQ=255)
                                                              tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagtgtcgcgaccatcaccgct                                >  1:5051003‑M1/1‑31 (MQ=255)
                                                              tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacagct                                >  2:2280588‑M1/1‑31 (MQ=255)
                                                                cataaaataatgtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgcttt                              <  1:4097697‑M1/89‑65 (MQ=255)
                                                                 ataaaataatgtGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttg                             >  2:2771231‑M1/4‑28 (MQ=255)
                                                                 ataaaataatgtGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttg                             >  1:1353430‑M1/4‑28 (MQ=255)
                                                                  taaaataatgTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttga                            <  1:3899055‑M1/91‑67 (MQ=255)
                                                                  taaaataatgTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttga                            >  1:5697493‑M1/3‑27 (MQ=255)
                                                                  taaaataatgTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttga                            <  2:5395160‑M1/91‑67 (MQ=255)
                                                                    aaataaTGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgcta                          >  2:1350843‑M1/1‑25 (MQ=255)
                                                                     aataatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatag                         >  2:5522406‑M1/1‑24 (MQ=255)
                                                                       taatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagta                       <  2:3156491‑M1/93‑72 (MQ=255)
                                                                       taatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagta                       <  2:2321578‑M1/93‑72 (MQ=255)
                                                                       taatGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagta                       <  1:3125850‑M1/93‑72 (MQ=255)
                                                                         atGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatg                     <  2:3066953‑M1/93‑74 (MQ=255)
                                                                          tgtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgc                    >  2:5943174‑M1/1‑19 (MQ=255)
                                                                          tgtgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgc                    <  1:565517‑M1/93‑75 (MQ=255)
                                                                            tgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcat                  >  1:4502323‑M1/1‑17 (MQ=255)
                                                                            tgATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcat                  <  1:5243515‑M1/93‑77 (MQ=255)
                                                                             gATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcata                 <  2:1280305‑M1/93‑78 (MQ=255)
                                                                             gATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcata                 <  1:2848816‑M1/93‑78 (MQ=255)
                                                                             gATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcata                 <  2:1057858‑M1/93‑78 (MQ=255)
                                                                             gATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcata                 >  1:489869‑M1/1‑16 (MQ=255)
                                                                               tGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacg               >  2:3225763‑M1/1‑14 (MQ=255)
                                                                               tGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacg               <  1:2222394‑M1/93‑80 (MQ=255)
                                                                                 ccGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggt             <  1:794280‑M1/93‑82 (MQ=255)
                                                                                  cGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgcccatcaccgctttgatagtatgcatacggtt            >  1:3356923‑M1/1‑11 (MQ=255)
                                                                                   gggATGGTTTgtctttcccggcatctttatagcgatagcaattttttgctgagcgtcgcgaccatcaccgctttgatagtatgcatacggttt           <  1:5492425‑M1/93‑84 (MQ=255)
                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttt           >  2:5046117‑M1/1‑10 (MQ=255)
                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttt           >  2:5773742‑M1/1‑10 (MQ=255)
                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttt           >  2:5995645‑M1/1‑10 (MQ=255)
                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttt           <  1:4373156‑M1/93‑84 (MQ=255)
                                                                                   gggATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttt           >  1:1037774‑M1/1‑10 (MQ=255)
                                                                                     gATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttttc         <  1:2620783‑M1/93‑86 (MQ=255)
                                                                                     gATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttttc         >  1:5626865‑M1/1‑8 (MQ=255)
                                                                                         gTTTgtatttccgggcatctttatagcgatagcaattatttactgagcgacgcgaccatcaccgctttgatagtatgcatacggttttccgcc     <  2:4189988‑M1/93‑90 (MQ=255)
                                                                                           ttgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttttccgcctg   <  2:2980180‑M1/93‑92 (MQ=255)
                                                                                            tgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgctttgatagtatgcatacggttttccgcctga  <  2:2921527‑M1/93‑93 (MQ=255)
                                                                                            |                                                                                            
AATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAA  >  NC_000913/4474757‑4474940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: