Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,826 +A intergenic (+1/+44) yjgL → / ← argI hypothetical protein/ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,8231.A99.5% 734.6 / ‑4.9 192*617K (TAA→AAA) yjgLhypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/0);  major base A (121/69);  minor base C (1/0);  total (123/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910
                                                                                             |                                                                                       
ttAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCat                                                                                          <  2:1138737/93‑1 (MQ=255)
  aaTGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataa                                                                                       >  2:5724938/1‑93 (MQ=255)
       aCAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCataaaataa                                                                                  >  1:4966393/1‑89 (MQ=255)
         aGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtg                                                                                 >  2:1649397/1‑93 (MQ=255)
           cGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAAtgtg                                                                               >  2:1327593/1‑93 (MQ=255)
            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  2:5260245/1‑93 (MQ=255)
            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  2:517581/1‑93 (MQ=255)
            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  2:2272855/1‑93 (MQ=255)
            gTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGa                                                                              >  1:2711521/1‑93 (MQ=255)
             tttAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGAt                                                                             >  2:3626169/1‑93 (MQ=255)
               taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGc                                                                           <  1:1133059/93‑1 (MQ=255)
               taaataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGc                                                                           <  1:5541059/93‑1 (MQ=255)
                  ataaaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCgg                                                                        >  2:2202743/1‑93 (MQ=255)
                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  1:3217460/1‑93 (MQ=255)
                     aaCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGAt                                                                     >  1:6133612/1‑93 (MQ=255)
                      aCTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATg                                                                    >  1:2404922/1‑93 (MQ=255)
                        tAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGt                                                                  >  2:1416606/1‑93 (MQ=255)
                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 <  1:3886221/93‑1 (MQ=255)
                         aaaTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGtt                                                                 <  1:6568318/93‑1 (MQ=255)
                           aTCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTg                                                               <  1:5349522‑M1/93‑2 (MQ=255)
                              cACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtat                                                            <  2:2836155‑M1/93‑5 (MQ=255)
                               aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatt                                                           >  1:1942719‑M1/1‑88 (MQ=255)
                               aCCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTGgtatt                                                           >  2:2458105‑M1/1‑88 (MQ=255)
                                 cGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                         >  1:1165131‑M1/1‑86 (MQ=255)
                                 cGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttc                                                         >  2:3154548‑M1/1‑86 (MQ=255)
                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        >  2:5099378‑M1/1‑85 (MQ=255)
                                  ggAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcc                                                        <  1:68753‑M1/93‑9 (MQ=255)
                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccg                                                       >  1:4677507‑M1/1‑84 (MQ=255)
                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       <  1:158812‑M1/93‑10 (MQ=255)
                                   gAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttccc                                                       >  2:4562305‑M1/1‑84 (MQ=255)
                                    aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      <  1:3641537‑M1/93‑11 (MQ=255)
                                    aaaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccg                                                      <  1:4551905‑M1/93‑11 (MQ=255)
                                     aaaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATATGATGCCGGGATGGTTTgtatttcccgg                                                     >  2:592722‑M1/1‑82 (MQ=255)
                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    >  1:4955981‑M1/1‑81 (MQ=255)
                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    >  1:712652‑M1/1‑81 (MQ=255)
                                      aaTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggc                                                    >  1:5947005‑M1/1‑81 (MQ=255)
                                       aTTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggca                                                   >  2:3842210‑M1/1‑80 (MQ=255)
                                        tttACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcat                                                  >  2:1734849‑M1/1‑79 (MQ=255)
                                          tACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatct                                                >  2:5548076‑M1/1‑77 (MQ=255)
                                           aCGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctt                                               >  2:2696955‑M1/1‑76 (MQ=255)
                                            cGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttt                                              >  1:5698435‑M1/1‑75 (MQ=255)
                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  2:4465661‑M1/1‑74 (MQ=255)
                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  1:3106203‑M1/93‑20 (MQ=255)
                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             <  1:4126139‑M1/93‑20 (MQ=255)
                                             gTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatcttta                                             >  1:5758545‑M1/1‑74 (MQ=255)
                                              tataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttat                                            <  1:2329251‑M1/93‑21 (MQ=255)
                                               ataGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttata                                           <  1:6430987‑M1/93‑22 (MQ=255)
                                                 aGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatatc                                         <  2:3286819‑M1/93‑24 (MQ=255)
                                                  gCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcg                                        <  1:2130077‑M1/93‑25 (MQ=255)
                                                    aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      >  1:2369783‑M1/1‑67 (MQ=255)
                                                    aaTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgat                                      <  2:487109‑M1/93‑27 (MQ=255)
                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     <  1:1103351‑M1/93‑28 (MQ=255)
                                                     aTAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgata                                     >  2:2062158‑M1/1‑66 (MQ=255)
                                                      tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  2:6065255‑M1/1‑65 (MQ=255)
                                                      tAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatag                                    >  1:1245829‑M1/1‑65 (MQ=255)
                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatcgcgatagc                                   <  1:6468093‑M1/93‑30 (MQ=255)
                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                   <  1:3801134‑M1/93‑30 (MQ=255)
                                                       aGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagc                                   <  1:2010479‑M1/93‑30 (MQ=255)
                                                        gAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagca                                  >  2:700465‑M1/1‑63 (MQ=255)
                                                         aaCACTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaa                                 >  2:1559795‑M1/1‑62 (MQ=255)
                                                          acacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaat                                >  1:2887758‑M1/1‑61 (MQ=255)
                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               <  1:4292587‑M1/93‑34 (MQ=255)
                                                           cacTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatt                               <  1:535913‑M1/93‑34 (MQ=255)
                                                            acTTTGGGTGGAAGAATAGACCTATCACTGCATCAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              >  2:6106261‑M1/1‑59 (MQ=255)
                                                            acTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaatta                              <  2:4662162‑M1/93‑35 (MQ=255)
                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:1735553‑M1/1‑58 (MQ=255)
                                                             cTTTGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattat                             >  1:6195660‑M1/1‑58 (MQ=255)
                                                                tGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                          <  2:351398‑M1/93‑39 (MQ=255)
                                                                tGGGTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattattta                          >  2:3997096‑M1/1‑55 (MQ=255)
                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:1848419‑M1/93‑40 (MQ=255)
                                                                 gggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttac                         <  2:1082099‑M1/93‑40 (MQ=255)
                                                                  ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        >  1:4472437‑M1/1‑53 (MQ=255)
                                                                  ggTGGAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttact                        >  2:3584891‑M1/1‑53 (MQ=255)
                                                                     ggAAGAATAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgag                     >  1:3670581‑M1/1‑50 (MQ=255)
                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                    >  1:1690875‑M1/1‑49 (MQ=255)
                                                                      gaagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagc                    <  2:3759696‑M1/93‑45 (MQ=255)
                                                                       aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg                   >  2:6301648‑M1/1‑48 (MQ=255)
                                                                       aagaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcg                   >  1:3487432‑M1/1‑48 (MQ=255)
                                                                        agaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgt                  >  1:4467252‑M1/1‑47 (MQ=255)
                                                                         gaaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtc                 <  2:4098179‑M1/93‑48 (MQ=255)
                                                                          aaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcg                <  2:71031‑M1/93‑49 (MQ=255)
                                                                          aaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcg                >  1:3530430‑M1/1‑45 (MQ=255)
                                                                          aaTAGACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcg                >  1:2671255‑M1/1‑45 (MQ=255)
                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            <  2:2094703‑M1/93‑53 (MQ=255)
                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            <  2:1495756‑M1/93‑53 (MQ=255)
                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            >  1:5889026‑M1/1‑41 (MQ=255)
                                                                              gACCTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgac            >  1:5098136‑M1/1‑41 (MQ=255)
                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          <  2:4218645‑M1/93‑55 (MQ=255)
                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          >  1:944002‑M1/1‑39 (MQ=255)
                                                                                ccTATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgacca          <  1:5433183‑M1/93‑55 (MQ=255)
                                                                                  tATCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccata        >  1:2080441‑M1/1‑37 (MQ=255)
                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac      >  2:4360695‑M1/1‑35 (MQ=255)
                                                                                    tCACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcac      >  2:5511674‑M1/1‑35 (MQ=255)
                                                                                     cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacc     <  2:4254008‑M1/93‑60 (MQ=255)
                                                                                     cACTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcacc     <  1:5185376‑M1/93‑60 (MQ=255)
                                                                                       cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   >  1:4521874‑M1/1‑32 (MQ=255)
                                                                                       cTGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgc   >  1:1403829‑M1/1‑32 (MQ=255)
                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  <  2:5402646‑M1/93‑63 (MQ=255)
                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  <  2:6521698‑M1/93‑63 (MQ=255)
                                                                                        tGCATAAAATAATGTGATGCCGGGATGGTTTgtatttcccggcatctttatagcgatagcaattatttactgagcgtcgcgaccatcaccgct  >  1:3369667‑M1/1‑31 (MQ=255)
                                                                                             |                                                                                       
TTAATGAACAGCGTTTAAATAAACTAAATCCACCGGAAAATTTACGTATAGCAATAGAACACTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTAT  >  NC_000913/4474731‑4474910

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