Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,966 T→A 15.3% intergenic (‑39/+18) ralR ← / ← recT Rac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,413,9660TA15.3% 533.0 / 75.5 203intergenic (‑39/+18)ralR/recTRac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base T (71/100);  new base A (17/14);  total (89/114)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01

Only 100 of 217 total aligned reads displayed.

CTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTT  >  NC_000913/1413867‑1414065
                                                                                                   |                                                                                                   
 tttCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCAtt                                                                                                    <  2:3593709/100‑1 (MQ=255)
 tttCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCAtt                                                                                                    <  2:3987374/100‑1 (MQ=255)
   tCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTac                                                                                                  <  1:117765/100‑1 (MQ=255)
    cACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCAGTGCATTaca                                                                                                 >  1:827384/1‑100 (MQ=255)
     aCCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTacac                                                                                                <  2:3460860/100‑1 (MQ=255)
      ccGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAcc                                                                                               >  2:345/1‑100 (MQ=255)
       cGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCg                                                                                              <  2:3892816/100‑1 (MQ=255)
       cGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCg                                                                                              >  1:2783787/1‑100 (MQ=255)
         ttACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGcc                                                                                            <  2:832812/100‑1 (MQ=255)
            cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaac                                                                                         <  2:1191241/100‑3 (MQ=255)
            cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAgg                                                                                         <  2:2907103/100‑1 (MQ=255)
            cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAgg                                                                                         >  2:4292065/1‑100 (MQ=255)
             tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                        <  2:3337219/100‑4 (MQ=255)
             tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGc                                                                                        <  2:654613/100‑1 (MQ=255)
              gATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaactt                                                                                       <  2:540128/100‑5 (MQ=255)
                tGGACAACCACAAAATGGACATGGTTTAACATTGTCATACCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGCCGGTGCACTACACCGCCAGCCTg                                                                                      >  1:619419/1‑99 (MQ=255)
                 ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAACAATATTTCAAGTTGGCGGTGCATTACACCGCCaacttgca                                                                                    >  2:3123637/1‑93 (MQ=255)
                  gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                   <  1:1832786/100‑1 (MQ=255)
                  gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                   >  1:1671372/1‑100 (MQ=255)
                   aCAACCAAAAAATGAACATGGTATAACATTCTCATAGATTATAATTTTTCAAATAAAAAATATTTCAAGTTGGCGGTGTAATGCTCCGCCaac                                                                                         <  1:500084/93‑3 (MQ=255)
                    cAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAatta                                                                                 >  2:2953562/1‑100 (MQ=255)
                    cAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAatta                                                                                 >  1:398518/1‑100 (MQ=255)
                      aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                               >  2:4214090/1‑100 (MQ=255)
                      aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                               >  2:451055/1‑100 (MQ=255)
                       ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                               >  1:1020718/1‑99 (MQ=255)
                       ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTc                                                                              <  1:4134507/100‑1 (MQ=255)
                        cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcc                                                                             <  1:4235019/100‑1 (MQ=255)
                         acaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCct                                                                            >  1:662739/1‑100 (MQ=255)
                          caAAATGTACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctc                                                                           <  1:3376994/100‑1 (MQ=255)
                            aaaTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTg                                                                         <  2:1070858/100‑1 (MQ=255)
                            aaaTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTg                                                                         <  2:2428904/100‑1 (MQ=255)
                            aaaTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACAGCCAGGCTGAATTATTCCCCTg                                                                         >  1:1036359/1‑100 (MQ=255)
                               tGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAt                                                                      <  2:1905359/100‑1 (MQ=255)
                                ggACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAAtt                                                                     >  1:446785/1‑100 (MQ=255)
                                ggACATGGTTTAACATTGTCATATCTCATAATTTATAAAAAAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGaa                                                                       <  1:2089042/98‑1 (MQ=255)
                                   cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATc                                                                  <  2:1292183/100‑1 (MQ=255)
                                   cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATc                                                                  >  2:2838953/1‑100 (MQ=255)
                                    aTGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                 <  2:3356571/100‑1 (MQ=255)
                                           aaCACTGTCATATCTCATAATTTTTCGTATAAAAAATATTTCATGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTAAATTATCGATTacac                                                          <  1:3760434/100‑1 (MQ=255)
                                             cATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                        <  2:2165896/100‑1 (MQ=255)
                                             cATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                        >  2:2487999/1‑100 (MQ=255)
                                              aTTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGt                                                       >  2:3948895/1‑100 (MQ=255)
                                              aTTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGt                                                       >  2:924368/1‑100 (MQ=255)
                                                 gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAtt                                                    >  2:157689/1‑100 (MQ=255)
                                                 gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAtt                                                    <  1:1392519/100‑1 (MQ=255)
                                                 gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAtt                                                    >  2:639855/1‑100 (MQ=255)
                                                   cATAGCTCATAACTTTTCTCATCAAACATATTTCAAGTTGGCGGTGCATTACACCGCCGGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcc                                                  <  1:3293233/100‑1 (MQ=255)
                                                     tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcccc                                                <  1:3206566/100‑1 (MQ=255)
                                                       tctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCgg                                              <  1:1270623/100‑1 (MQ=255)
                                                        ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                             <  2:2321551/100‑1 (MQ=255)
                                                        ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                             <  2:4118051/100‑1 (MQ=255)
                                                         tcATAATTTTTCTCATAAAAAATATTTCAATTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                            <  1:851036/100‑1 (MQ=255)
                                                         tcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                            <  2:662741/100‑1 (MQ=255)
                                                         tcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                            <  2:971398/100‑1 (MQ=255)
                                                           aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTaa                                          <  2:146945/100‑1 (MQ=255)
                                                           aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTaa                                          >  1:3158835/1‑100 (MQ=255)
                                                              aTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAc                                       <  1:2241700/100‑1 (MQ=255)
                                                                  ttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATGCagag                                   >  2:1692151/1‑100 (MQ=255)
                                                                   tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                  <  1:157688/100‑1 (MQ=255)
                                                                   tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                  <  1:3211199/100‑1 (MQ=255)
                                                                   tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                  <  2:2221759/100‑1 (MQ=255)
                                                                    ctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTTAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                   <  1:2703302/98‑1 (MQ=255)
                                                                       aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATc                              >  1:2979054/1‑100 (MQ=255)
                                                                        tAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCt                             >  2:4386142/1‑100 (MQ=255)
                                                                         aaaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTg                            <  2:374349/100‑1 (MQ=255)
                                                                         aaaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTg                            >  2:3173003/1‑100 (MQ=255)
                                                                          aaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                           >  1:3521450/1‑100 (MQ=255)
                                                                          aaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                           <  1:4292034/100‑1 (MQ=255)
                                                                           aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCa                          <  2:351752/100‑1 (MQ=255)
                                                                           aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCa                          >  2:723098/1‑100 (MQ=255)
                                                                           aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCa                          <  1:3651654/100‑1 (MQ=255)
                                                                            aaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAg                         >  2:2124781/1‑100 (MQ=255)
                                                                             aaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAgg                        <  1:2953546/100‑1 (MQ=255)
                                                                                atTTCAAGTTGGCGGTGCATTACACCGCGAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                     <  2:2073874/100‑1 (MQ=255)
                                                                                 attcagcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGCATTATCGATTACACTGTATTCCACGGTTACTACAGAGGAATCTGCAGGATc                     >  1:4046469/9‑99 (MQ=25)
                                                                                  ttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGa                   <  1:4129447/100‑1 (MQ=255)
                                                                                    cAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGAtt                 <  1:98507/100‑1 (MQ=255)
                                                                                     agcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                >  1:1711013/5‑100 (MQ=255)
                                                                                     agcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                >  2:1320308/5‑100 (MQ=255)
                                                                                      gcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                >  1:2066299/4‑99 (MQ=255)
                                                                                       gTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc              <  2:1361427/100‑1 (MQ=255)
                                                                                       cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTGTTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc              >  1:1220462/3‑100 (MQ=255)
                                                                                       cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc              >  1:947007/3‑100 (MQ=255)
                                                                                       cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc              >  2:2041072/3‑100 (MQ=255)
                                                                                        ctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCATGATCGATTGTCa             >  1:2619924/2‑100 (MQ=255)
                                                                                        ctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCa             >  2:827727/2‑100 (MQ=255)
                                                                                         tGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAg            >  1:1130223/1‑100 (MQ=255)
                                                                                          ggCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt           <  2:2728934/100‑1 (MQ=255)
                                                                                          ggCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt           >  2:1252324/1‑100 (MQ=255)
                                                                                           gCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTg          <  2:3227831/100‑1 (MQ=255)
                                                                                           gCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTg          <  2:386355/100‑1 (MQ=255)
                                                                                            cGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTgg         >  1:2866262/1‑100 (MQ=255)
                                                                                            cGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTgg         <  1:1451582/100‑1 (MQ=255)
                                                                                             ggTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGt        <  2:1334280/100‑1 (MQ=255)
                                                                                             ggTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGt        >  2:2537580/1‑100 (MQ=255)
                                                                                               tGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATGACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTc      <  1:667828/100‑1 (MQ=255)
                                                                                                  aTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt   >  1:769859/1‑100 (MQ=255)
                                                                                                  aTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt   <  1:3724053/100‑1 (MQ=255)
                                                                                                   ttACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt  >  2:516431/1‑100 (MQ=255)
                                                                                                   ttACACCGCCAGGCTCAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt  <  2:4335577/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CTTTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTT  >  NC_000913/1413867‑1414065

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