Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,413,968 A→G 17.8% 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,9680AG17.8% 326.4 / 58.9 148intergenic (‑41/+16)ralR/recTRac prophage; restriction alleviation protein/Rac prophage; recombination and repair protein
Reads supporting (aligned to +/- strand):  ref base A (59/62);  new base G (8/18);  total (67/81)
Fisher's exact test for biased strand distribution p-value = 1.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

Only 100 of 156 total aligned reads displayed.

TTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTC  >  NC_000913/1413869‑1414067
                                                                                                   |                                                                                                   
ttCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTa                                                                                                     <  1:2450046/100‑1 (MQ=255)
  cACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTATCTCAGAAAAAATATTTCAAGTTGGCGGTGCATTaca                                                                                                   >  1:2513584/1‑100 (MQ=255)
   aCCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTacac                                                                                                  <  2:1144838/100‑1 (MQ=255)
       ttccgatctgACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGCTGGCGGTGTAATGCACCGcc                                                                                              <  1:3348544/91‑1 (MQ=255)
        tACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCa                                                                                             >  1:680547/1‑100 (MQ=255)
        tACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCa                                                                                             >  2:2211758/1‑100 (MQ=255)
          cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaac                                                                                           <  1:1561257/100‑3 (MQ=255)
          cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaac                                                                                           <  2:3655258/100‑3 (MQ=255)
          cTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAgg                                                                                           >  2:1282024/1‑100 (MQ=255)
           tggtGGACAACCACAAAATGTACATGGTTTAACATTTTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                          <  1:2291026/97‑4 (MQ=255)
           tGATGGACAACCACAGAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGc                                                                                          <  2:3679645/100‑1 (MQ=255)
           tGATGGACAACCACAAAATGTACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                          <  2:1095154/100‑4 (MQ=255)
           tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                          <  2:217302/100‑4 (MQ=255)
           tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaact                                                                                          <  2:3128730/100‑4 (MQ=255)
           tGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGc                                                                                          <  1:2202035/100‑1 (MQ=255)
            gATGGACAATCACAAAATGGACATGGTTTAACATTGTTATATCTCATGATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaactt                                                                                         <  1:2536874/100‑5 (MQ=255)
              tGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaacttga                                                                                       <  2:680550/100‑7 (MQ=38)
                gACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAat                                                                                     >  1:1953067/1‑100 (MQ=255)
                 aCAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGTAATGCACCGCCaacttgaaat                                                                                    <  1:977635/100‑10 (MQ=37)
                   aaCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattat                                                                                  <  1:1552886/100‑1 (MQ=255)
                    aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                 >  1:3450426/1‑100 (MQ=255)
                    aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                 >  1:293240/1‑100 (MQ=255)
                    aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                 >  2:1985131/1‑100 (MQ=255)
                    aCCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAattatt                                                                                 >  2:73773/1‑100 (MQ=255)
                     ccACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTc                                                                                >  1:113666/1‑100 (MQ=255)
                      cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcc                                                                               >  2:2931062/1‑100 (MQ=255)
                      cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcc                                                                               <  2:464978/100‑1 (MQ=255)
                      cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcc                                                                               >  2:2932047/1‑100 (MQ=255)
                      cacaAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTcc                                                                               >  2:3259393/1‑100 (MQ=255)
                        caAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCATGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctc                                                                             <  1:2861444/100‑1 (MQ=255)
                        caAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCctc                                                                             >  1:3396322/1‑100 (MQ=255)
                          aaaTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTg                                                                           <  1:43114/100‑1 (MQ=255)
                            aTGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGaa                                                                         <  2:1489692/100‑1 (MQ=255)
                               gACATGGTTTAACATTGTCATCTCTCATACTTTTTCTAATAAAAAATATTTCATGTTATCGGTGCATTACACCGCCAGGCTGAATTATTCCTTTTAAtt                                                                       >  1:2000849/1‑99 (MQ=255)
                                 cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATc                                                                    <  2:359144/100‑1 (MQ=255)
                                 cATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATc                                                                    >  1:1666766/1‑100 (MQ=255)
                                  aTGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCg                                                                   <  1:2503265/100‑1 (MQ=255)
                                  aTGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTAGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATc                                                                    <  1:3620230/99‑1 (MQ=255)
                                   tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGa                                                                  >  2:1724004/1‑100 (MQ=255)
                                   tGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGa                                                                  >  2:2145969/1‑100 (MQ=255)
                                     gTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGAtt                                                                >  2:2598890/1‑100 (MQ=255)
                                      tttAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTa                                                               <  2:2756736/100‑1 (MQ=255)
                                       ttAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTac                                                              >  2:3281213/1‑100 (MQ=255)
                                         aaCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTTCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTacac                                                            <  1:2041760/100‑1 (MQ=255)
                                          aCATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCACCAGGCTGAATTATTCCTCTGAATTATCGATTACACt                                                           >  1:471856/1‑100 (MQ=255)
                                           cATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                          <  2:3034349/100‑1 (MQ=255)
                                           cATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTg                                                          >  1:1145015/1‑100 (MQ=255)
                                              tGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAt                                                       <  2:3337081/100‑1 (MQ=255)
                                              tGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAt                                                       <  2:3442027/100‑1 (MQ=255)
                                               gTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTAtt                                                      >  2:3571099/1‑100 (MQ=255)
                                                   tatCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcccc                                                  <  2:2544659/100‑1 (MQ=255)
                                                     tctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCgg                                                <  1:2337140/100‑1 (MQ=255)
                                                     tctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACAACGCCAGGCTGAATTATTCTTCTGAATTATCGATTACACTGTATTctc                                                   >  1:3728646/1‑95 (MQ=255)
                                                      ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTcccc                                                  >  1:328001/1‑97 (MQ=255)
                                                      ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                               <  2:609411/100‑1 (MQ=255)
                                                      ctcATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGt                                               >  1:13226/1‑100 (MQ=255)
                                                        cATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTa                                             <  2:113214/100‑1 (MQ=255)
                                                         aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCTGGTTaa                                            >  1:3085599/1‑100 (MQ=255)
                                                         aTAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCACCGGTTaa                                            >  1:3199165/1‑100 (MQ=255)
                                                          tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAt                                           <  2:2031705/100‑1 (MQ=255)
                                                          tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAt                                           >  2:3339167/1‑100 (MQ=255)
                                                          tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAt                                           <  2:1953077/100‑1 (MQ=255)
                                                          tAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAAt                                           <  2:91746/100‑1 (MQ=255)
                                                              ttttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACag                                       <  1:1894947/100‑1 (MQ=255)
                                                              ttttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACag                                       >  1:2467102/1‑100 (MQ=255)
                                                              ttttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACa                                        >  1:3007250/1‑99 (MQ=255)
                                                                ttCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACagag                                     <  1:1869289/100‑1 (MQ=255)
                                                                   tcATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGaa                                  >  1:2103391/1‑100 (MQ=255)
                                                                    cATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAAt                                 >  2:2389663/1‑100 (MQ=255)
                                                                     aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAAt                                 >  2:2835308/1‑99 (MQ=255)
                                                                     aTAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATc                                >  2:2353622/1‑100 (MQ=255)
                                                                      tAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCt                               >  2:1913273/1‑100 (MQ=255)
                                                                      tAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCt                               >  2:2630738/1‑100 (MQ=255)
                                                                      tAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATAGATTACACTGTATTCCCCGGTTAATACAGAGGAATCt                               >  1:2523730/1‑100 (MQ=255)
                                                                        aaaaaTATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGc                             <  2:3378980/100‑1 (MQ=255)
                                                                              atTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                       <  1:632754/100‑1 (MQ=255)
                                                                              atTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATc                       <  2:900595/100‑1 (MQ=255)
                                                                               tttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCg                      <  2:1811570/100‑1 (MQ=255)
                                                                                ttCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGa                     >  1:1679064/1‑100 (MQ=255)
                                                                                 tcagcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGAt                    >  2:610423/7‑100 (MQ=38)
                                                                                  cAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTACCCGGTTAATACAGAGGAATCTGCAGGATCGAtt                   >  1:3367596/1‑100 (MQ=255)
                                                                                    gcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGt                 >  1:197277/4‑100 (MQ=255)
                                                                                    gcctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGACGAATCTGCAGGATCGATTGt                 >  1:197266/4‑100 (MQ=255)
                                                                                     cctGGCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTc                >  2:957653/3‑100 (MQ=255)
                                                                                       tGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAg              <  2:3275021/100‑1 (MQ=255)
                                                                                        ggCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt             >  2:2041651/1‑100 (MQ=255)
                                                                                        ggCGGTGTAATGCACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt             >  2:1731430/1‑100 (MQ=255)
                                                                                        ggCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt             >  2:1860440/1‑100 (MQ=255)
                                                                                        ggCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGt             <  1:3664529/100‑1 (MQ=255)
                                                                                         gCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTg            <  2:148607/100‑1 (MQ=255)
                                                                                          cGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTgg           >  2:1788238/1‑100 (MQ=255)
                                                                                           ggTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGt          <  1:1282018/100‑1 (MQ=255)
                                                                                               cATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCt      >  2:1763754/1‑100 (MQ=255)
                                                                                               cATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCt      >  1:506483/1‑100 (MQ=255)
                                                                                                aTTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt     <  2:1666773/100‑1 (MQ=255)
                                                                                                 ttACACTTCTAGGCTTAATTATTCCTCTGAATTATCGATTAGACTGTATTCCCCGGTTAATACATATGAATCTTCAGGTTCGATTGTCAGTGGTTTCttt    <  1:1362811/100‑1 (MQ=255)
                                                                                                 ttACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCttt    <  1:1752366/100‑1 (MQ=255)
                                                                                                 ttACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtt     <  1:129892/99‑1 (MQ=255)
                                                                                                  tACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCtttt   <  1:1569404/100‑1 (MQ=255)
                                                                                                   acacCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTc  >  2:1825679/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TTCACCGTTACTGATGGACAACCACAAAATGGACATGGTTTAACATTGTCATATCTCATAATTTTTCTCATAAAAAATATTTCAAGTTGGCGGTGCATTACACCGCCAGGCTGAATTATTCCTCTGAATTATCGATTACACTGTATTCCCCGGTTAATACAGAGGAATCTGCAGGATCGATTGTCAGTGGTTCCTTTTC  >  NC_000913/1413869‑1414067

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