Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 708,658 T→C W368R (TGG→CGG)  ybfM → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913708,6580TC87.5% 721.0 / 82.4 263W368R (TGG→CGG) ybfMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/17);  new base C (112/118);  total (128/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 266 total aligned reads displayed.

CGCCTGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGCGTAC  >  NC_000913/708568‑708747
                                                                                          |                                                                                         
  ccTGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGaa                                                                                       >  2:3583050/1‑93 (MQ=255)
   cTGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAAt                                                                                      <  2:2433431/93‑1 (MQ=255)
   cTGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAAt                                                                                      <  2:626049/93‑1 (MQ=255)
   cTGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAAt                                                                                      >  1:3171536/1‑93 (MQ=255)
    tGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATc                                                                                     <  1:6123150/93‑1 (MQ=255)
     ggATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCt                                                                                    <  2:6334581/93‑1 (MQ=255)
     ggATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCt                                                                                    >  2:5753474/1‑93 (MQ=255)
       atatCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTc                                                                                  <  1:5527015/93‑1 (MQ=255)
         atCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCa                                                                                <  1:4352253/93‑1 (MQ=255)
          tCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAg                                                                               <  2:4758364/93‑1 (MQ=255)
           cTGGTGGGATAACCGTTCTGACTTCAACGCCAACGTCGAAAAAGCTGTCTGCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAgg                                                                              <  2:1109561/93‑1 (MQ=255)
           cTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAgg                                                                              <  2:845737/93‑1 (MQ=255)
                gggATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCg                                                                         >  2:5215251/1‑93 (MQ=255)
                  gATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGcc                                                                       <  1:3558157/93‑1 (MQ=255)
                    tAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCAt                                                                     >  1:978433/1‑93 (MQ=255)
                    tAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCAt                                                                     >  2:4878008/1‑93 (MQ=255)
                    tAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCAt                                                                     >  2:297975/1‑93 (MQ=255)
                     aaCCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATc                                                                    >  1:2311637/1‑93 (MQ=255)
                      aCCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCg                                                                   >  2:4307383/1‑93 (MQ=255)
                       ccGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCgg                                                                  <  1:1001234/93‑1 (MQ=255)
                       ccGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCgg                                                                  >  1:3478951/1‑93 (MQ=255)
                        cGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGgc                                                                 >  2:5628583/1‑93 (MQ=255)
                        cGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGgc                                                                 >  1:754242/1‑93 (MQ=255)
                         gTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGgcg                                                                <  1:3811298/93‑1 (MQ=255)
                         gTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGgcg                                                                >  2:2271043/1‑93 (MQ=255)
                         gTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGgcg                                                                >  1:551838/1‑93 (MQ=255)
                          ttCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGgcgc                                                               >  2:5928678/1‑93 (MQ=255)
                          ttCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGgcgc                                                               >  2:2132011/1‑93 (MQ=255)
                           tCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCt                                                              <  1:1041853/93‑1 (MQ=255)
                            cTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCtt                                                             >  2:3447777/1‑93 (MQ=255)
                             tGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTc                                                            <  2:176997/93‑1 (MQ=255)
                             cgACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTc                                                            <  1:6253949/92‑1 (MQ=255)
                              gACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTcc                                                           >  1:6444183/1‑93 (MQ=255)
                              gACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTcc                                                           <  1:5587313/93‑1 (MQ=255)
                               aCTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCt                                                          <  2:1730091/93‑1 (MQ=255)
                                cTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTa                                                         >  2:2061499/1‑93 (MQ=255)
                                   cAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGt                                                      <  2:3369092/93‑1 (MQ=255)
                                     aCGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAATCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGttt                                                    >  2:6088618/1‑93 (MQ=255)
                                     aCGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGttt                                                    >  1:5879265/1‑93 (MQ=255)
                                       gCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTAc                                                  >  1:5657045/1‑93 (MQ=255)
                                          aaCGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCa                                               >  2:3529830/1‑93 (MQ=255)
                                          aaCGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCa                                               <  1:1572346/93‑1 (MQ=255)
                                           aCGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCAt                                              <  2:135320/93‑1 (MQ=255)
                                           aCGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCAt                                              >  1:4554196/1‑93 (MQ=255)
                                             ggCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATgg                                            >  1:1150582/1‑93 (MQ=255)
                                             ggCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATgg                                            <  2:1922513/93‑1 (MQ=255)
                                              gCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATggg                                           <  2:892412/93‑1 (MQ=255)
                                               cGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGa                                          <  2:584574/93‑1 (MQ=255)
                                                gAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGAt                                         >  1:1771986/1‑93 (MQ=255)
                                                gAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGAt                                         <  1:5094731/93‑1 (MQ=255)
                                                 aaaaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATg                                        <  1:1878338/93‑1 (MQ=255)
                                                 aaaaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATg                                        >  2:1988832/1‑93 (MQ=255)
                                                 aaaaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATg                                        >  2:3846272/1‑93 (MQ=255)
                                                  aaaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCAATCGGCGCTTCCTACGTTTACGCATGGGATGc                                       <  2:5516843/93‑1 (MQ=255)
                                                   aaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCt                                      <  1:1447028/93‑1 (MQ=255)
                                                   aaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCt                                      >  1:1697537/1‑93 (MQ=255)
                                                   aaaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCt                                      >  2:3058237/1‑93 (MQ=255)
                                                    aaGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTa                                     <  2:4494770/93‑1 (MQ=255)
                                                     aGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGTCGCTTCCTACGTTTACGCATGGGCTGCTaa                                    >  2:1380380/1‑93 (MQ=255)
                                                      gCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTaaa                                   >  2:5285019/1‑93 (MQ=255)
                                                        ggTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAAcc                                 <  1:3886231/93‑1 (MQ=255)
                                                        ggTCTTCTTCGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAAcc                                 <  1:5010589/93‑1 (MQ=255)
                                                            ttcttcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCg                             <  2:4022864/93‑1 (MQ=255)
                                                             tcttcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGa                            <  1:245433/93‑1 (MQ=255)
                                                              cttcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGAc                           <  1:3716409/93‑1 (MQ=255)
                                                               ttcGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCGTGGGATGCTAAACCTGCGAcc                          <  1:2985949/93‑1 (MQ=255)
                                                               ttcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGAcc                          <  2:4590342/93‑1 (MQ=255)
                                                               ttcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGAcc                          <  1:6587376/93‑1 (MQ=255)
                                                               ttcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGAcc                          <  1:622584/93‑1 (MQ=255)
                                                                tcGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCt                         >  2:4180504/1‑93 (MQ=255)
                                                                tcGGTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGAGGCTAAACCTGCGACCt                         >  1:3351941/1‑93 (MQ=255)
                                                                  ggTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTgg                       >  1:6294889/1‑93 (MQ=255)
                                                                  ggTGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTgg                       <  2:741719/93‑1 (MQ=255)
                                                                    tGCGATGTATGACCTGAAAAGCCGGAATCTTCCGGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCa                     <  1:6022652/93‑1 (MQ=255)
                                                                    tGCGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCa                     <  1:6635168/93‑1 (MQ=255)
                                                                      ggATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCaga                   <  1:4742691/92‑1 (MQ=255)
                                                                      cGATGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCaga                   >  1:2660152/1‑93 (MQ=255)
                                                                       gATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCagag                  >  2:1604220/1‑93 (MQ=255)
                                                                        aTGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGc                 <  2:2918042/93‑1 (MQ=255)
                                                                         tGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCa                <  1:5798280/93‑1 (MQ=255)
                                                                         tGTATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCa                >  2:2952210/1‑93 (MQ=255)
                                                                           tATGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAAt              <  1:200557/93‑1 (MQ=255)
                                                                            aTGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATc             <  2:6530726/93‑1 (MQ=255)
                                                                            aTGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATc             <  2:2714925/93‑1 (MQ=255)
                                                                            aTGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATc             <  2:5085132/93‑1 (MQ=255)
                                                                             tGACCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATcc            <  2:2812862/93‑1 (MQ=255)
                                                                               aCCTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCgg          <  2:4646275/93‑1 (MQ=255)
                                                                                ccTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGa         <  2:5387173/93‑1 (MQ=255)
                                                                                ccTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGa         <  1:5412821/93‑1 (MQ=255)
                                                                                ccTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGa         <  1:4619048/93‑1 (MQ=255)
                                                                                 cTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGAt        <  1:739390/93‑1 (MQ=255)
                                                                                 cTGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGAt        >  2:638922/1‑93 (MQ=255)
                                                                                  tGAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATACTAAACCTGCGACCTGGCAGAGCAATCCGGATg       >  1:431500/1‑93 (MQ=255)
                                                                                   gAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGc      <  1:4137456/93‑1 (MQ=255)
                                                                                   gAAAAACCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGc      >  2:4086974/1‑93 (MQ=255)
                                                                                    aaaaaCCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGCg     <  2:515882/93‑1 (MQ=255)
                                                                                    aaaaaCCGGAATCTTCCAGGCTTAGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGCg     >  1:3583707/1‑93 (MQ=255)
                                                                                     aaaaCCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGCGt    <  1:5167837/93‑1 (MQ=255)
                                                                                     aaaaCCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTCGCAGAGCAATCCGGATGCGt    >  1:1317352/1‑93 (MQ=255)
                                                                                       aaCCGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGCGtac  <  2:3617691/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCCTGGATATCTGGTGGGATAACCGTTCTGACTTCAACGCCAACGGCGAAAAAGCGGTCTTCTTCGGTGCGATGTATGACCTGAAAAACTGGAATCTTCCAGGCTTCGCCATCGGCGCTTCCTACGTTTACGCATGGGATGCTAAACCTGCGACCTGGCAGAGCAATCCGGATGCGTAC  >  NC_000913/708568‑708747

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