Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,968 A→G 15.1% intergenic (‑41/+16) 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,9680AG15.1% 492.0 / 72.9 205intergenic (‑41/+16)ralR/recTRac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base A (74/100);  new base G (17/14);  total (91/114)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 209 total aligned reads displayed.

TTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTT  >  NC_000913/1413869‑1414066
                                                                                                   |                                                                                                  
ttCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTa                                                                                                    <  1:3711617/100‑1 (MQ=255)
   aCCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTacac                                                                                                 <  2:3460860/100‑1 (MQ=255)
    ccGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAcc                                                                                                >  2:345/1‑100 (MQ=255)
    ccGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAcc                                                                                                <  2:2766123/100‑1 (MQ=255)
     cGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCg                                                                                               >  1:2783787/1‑100 (MQ=255)
      gTTACTGATGTACAACCACAAAATATACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAcc                                                                                                <  1:858831/98‑1 (MQ=255)
         aCTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAg                                                                                           <  2:3142668/100‑1 (MQ=255)
         aCTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAg                                                                                           <  2:613183/100‑1 (MQ=255)
          cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAgg                                                                                          >  2:4292065/1‑100 (MQ=255)
           tgatctaCAACCACAAAATGTACATTGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                         <  1:4259424/94‑4 (MQ=255)
           tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                         <  2:3337219/100‑4 (MQ=255)
           tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGc                                                                                         <  2:654613/100‑1 (MQ=255)
             aTGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaacttg                                                                                       <  2:974429/100‑6 (MQ=38)
              tGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTAATAAAAAATATTTAAAGTTGGCGGTGCATTACACCGCCAGGATGa                                                                                      >  1:3740369/1‑100 (MQ=255)
              tGGACAACCACAAAATGGACATGGTTTAACATTGTCATACCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGCCGGTGCACTACACCGCCAGCCTg                                                                                       >  1:619419/1‑99 (MQ=255)
               ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGaa                                                                                     >  2:2242111/1‑100 (MQ=255)
               ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGaa                                                                                     <  2:45959/100‑1 (MQ=255)
               ggACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGCTGGCGGTGTAATGCACCGCCaacttgaa                                                                                     <  2:2923742/100‑8 (MQ=37)
                gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                    <  1:1832786/100‑1 (MQ=255)
                gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                    >  1:1671372/1‑100 (MQ=255)
                gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                    <  2:2192502/100‑1 (MQ=255)
                gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTATATCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaacttgaaa                                                                                    <  1:815511/100‑9 (MQ=37)
                 aCAACCAAAAAATGAACATGGTATAACATTCTCATAGATTATAATTTTTCAAATAAAAAATATTTCAAGTTGGCGGTGTAATGCTCCGCCaac                                                                                          <  1:500084/93‑3 (MQ=255)
                  cAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAatta                                                                                  >  1:398518/1‑100 (MQ=255)
                    aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                >  2:4214090/1‑100 (MQ=255)
                     ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                >  1:1020718/1‑99 (MQ=255)
                     ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTc                                                                               >  2:2113411/1‑100 (MQ=255)
                     ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTc                                                                               >  1:1205269/1‑100 (MQ=255)
                     ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTc                                                                               <  1:4134507/100‑1 (MQ=255)
                      cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTAAAAACAAAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcc                                                                              <  1:3209298/100‑1 (MQ=255)
                       acaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCct                                                                             >  1:662739/1‑100 (MQ=255)
                         aaaaTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctct                                                                           <  2:75042/100‑1 (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)
                                  aTGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                  <  2:3356571/100‑1 (MQ=255)
                                  aTGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                  <  2:828277/100‑1 (MQ=255)
                                   tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGa                                                                 >  2:3885813/1‑100 (MQ=255)
                                   tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGa                                                                 >  2:959701/1‑100 (MQ=255)
                                       ttAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTac                                                             >  1:3422963/1‑100 (MQ=255)
                                         aagaTTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTacac                                                           <  1:3325469/97‑1 (MQ=255)
                                         aaCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTacac                                                           <  2:3777533/100‑1 (MQ=255)
                                          aCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACt                                                          >  1:1255060/1‑100 (MQ=255)
                                           cATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                         <  2:2340989/100‑1 (MQ=255)
                                           cATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                         <  2:398519/100‑1 (MQ=255)
                                               gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAtt                                                     >  2:2146275/1‑100 (MQ=255)
                                               gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAtt                                                     <  1:1392519/100‑1 (MQ=255)
                                                 cATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcc                                                   >  2:2473513/1‑100 (MQ=255)
                                                   tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcccc                                                 <  1:3836011/100‑1 (MQ=255)
                                                      ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                              >  2:119503/1‑100 (MQ=255)
                                                      ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                              <  2:4118051/100‑1 (MQ=255)
                                                       tcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                             <  2:662741/100‑1 (MQ=255)
                                                       taaTAATTTTTCTCATCCAAAATATTTCCAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGtt                                             <  1:363273/98‑1 (MQ=255)
                                                          tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAt                                          <  1:3160713/100‑1 (MQ=255)
                                                            aTTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATAc                                        <  1:2241700/100‑1 (MQ=255)
                                                               tttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACaga                                     >  1:1109521/1‑100 (MQ=255)
                                                               tttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACaga                                     >  1:1171290/1‑100 (MQ=255)
                                                                ttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                    <  2:2515724/100‑1 (MQ=255)
                                                                 tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                   <  2:869473/100‑1 (MQ=255)
                                                                 tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                   <  1:157688/100‑1 (MQ=255)
                                                                 tctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAgg                                   <  2:3069001/100‑1 (MQ=255)
                                                                  ctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTTAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                    <  1:2703302/98‑1 (MQ=255)
                                                                  ctcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGa                                  >  2:3661478/1‑100 (MQ=255)
                                                                     aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATc                               >  1:2979054/1‑100 (MQ=255)
                                                                     aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATc                               >  1:849178/1‑100 (MQ=255)
                                                                       aaaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTg                             >  2:3173003/1‑100 (MQ=255)
                                                                         aaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCa                           <  2:351752/100‑1 (MQ=255)
                                                                           aaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAgg                         <  2:690864/100‑1 (MQ=255)
                                                                             tatTTCAAGTTGGCGGTGCATTACTCCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAg                          <  2:1661948/97‑1 (MQ=255)
                                                                              atTTCAAGTTGGCGGTGCATTACACCGCGAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                      <  2:2073874/100‑1 (MQ=255)
                                                                              atTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                      <  2:111963/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                 >  2:1320308/5‑100 (MQ=255)
                                                                                   agcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                 >  1:1711013/5‑100 (MQ=255)
                                                                                    gcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                 >  1:2066299/4‑99 (MQ=255)
                                                                                     gTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc               <  2:440248/100‑1 (MQ=255)
                                                                                     gTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc               <  2:1361427/100‑1 (MQ=255)
                                                                                     cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc               >  2:2041072/3‑100 (MQ=255)
                                                                                     cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc               >  2:3944292/3‑100 (MQ=255)
                                                                                      ttGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCa              <  1:2844289/100‑1 (MQ=255)
                                                                                      ctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCATGATCGATTGTCa              >  1:2619924/2‑100 (MQ=255)
                                                                                      ctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTg                 >  1:3261484/2‑97 (MQ=255)
                                                                                        ggCGGTGTAATGCACCGCCAGGCTGAATTATTCGTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATAGATTGTCAGt            >  1:2886781/1‑100 (MQ=255)
                                                                                        ggCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt            >  2:1252324/1‑100 (MQ=255)
                                                                                        ggCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt            >  2:538823/1‑100 (MQ=255)
                                                                                         gCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTg           <  2:3227831/100‑1 (MQ=255)
                                                                                          cGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTgg          <  1:1451582/100‑1 (MQ=255)
                                                                                          cGGTGCATTACACCGCCAGGCTGAATTATTCCCCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTgg          <  1:4301633/100‑1 (MQ=255)
                                                                                           ggTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGt         <  2:1334280/100‑1 (MQ=255)
                                                                                           ggTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGt         >  2:2537580/1‑100 (MQ=255)
                                                                                             tGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTc       <  2:2866278/100‑1 (MQ=255)
                                                                                             tGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATGACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTc       <  1:667828/100‑1 (MQ=255)
                                                                                               cATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCt     <  2:1413379/100‑1 (MQ=255)
                                                                                                aTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt    >  1:769859/1‑100 (MQ=255)
                                                                                                aTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt    >  2:1520463/1‑100 (MQ=255)
                                                                                                 ttACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt   <  2:4380768/100‑1 (MQ=255)
                                                                                                 ttACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt   <  1:1339034/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTT  >  NC_000913/1413869‑1414066

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