Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,599,353 G→A A30V (GCC→GTC)  opgB ← OPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,599,3530GA97.1% 580.6 / 1.5 173A30V (GCC→GTC) opgBOPG periplasmic biosynthetic phosphoglycerol transferases I (membrane‑bound) and II (soluble)
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (85/83);  total (88/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 173 total aligned reads displayed.

ATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGA  >  NC_000913/4599256‑4599448
                                                                                                 |                                                                                               
 taGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTg                                                                                              >  2:212185/1‑100 (MQ=255)
   gagaACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCa                                                                                            >  1:3820360/1‑100 (MQ=255)
   gagaACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCa                                                                                            >  1:1234134/1‑100 (MQ=255)
   gagaACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCa                                                                                            >  1:1014745/1‑100 (MQ=255)
   gagaACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCa                                                                                            >  1:4074096/1‑100 (MQ=255)
   gagaACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCa                                                                                            >  1:5167074/1‑100 (MQ=255)
        ccGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAacca                                                                                       >  2:3893886/1‑100 (MQ=255)
        ccGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAacca                                                                                       <  2:2098738/100‑1 (MQ=255)
             tcgtTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAgg                                                                                  <  1:3177307/100‑1 (MQ=255)
               gtTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGgtg                                                                                >  2:1596774/1‑100 (MQ=255)
               gtTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGgtg                                                                                <  1:646042/100‑1 (MQ=255)
               gtTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGgtg                                                                                >  2:3324446/1‑100 (MQ=255)
               gtTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGgtg                                                                                >  1:4757782/1‑100 (MQ=255)
               gtTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAACCCCAGCACCGTTAACGTGACTGCAAACAACCAGgtg                                                                                >  2:4268178/1‑100 (MQ=255)
                 tAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGtt                                                                              <  2:455377/100‑1 (MQ=255)
                  aaTACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTa                                                                             <  2:1647039/100‑1 (MQ=255)
                  aaTACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTa                                                                             <  1:2567651/100‑1 (MQ=255)
                   aTACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTAc                                                                            <  2:907072/100‑1 (MQ=255)
                    tACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACg                                                                           >  2:3553774/1‑100 (MQ=255)
                     aCCAGCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCTCCGTTAACGTGACTGCAAACCACCAGGTGTTACGt                                                                          >  2:4640376/1‑100 (MQ=255)
                      ccATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTc                                                                         >  1:4253862/1‑100 (MQ=255)
                       cATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTcc                                                                        <  1:4537915/100‑1 (MQ=255)
                       cATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTcc                                                                        <  1:1555212/100‑1 (MQ=255)
                        aTCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTccc                                                                       >  1:4746227/1‑100 (MQ=255)
                         tCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCg                                                                      <  1:1102344/100‑1 (MQ=255)
                           gCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCt                                                                    >  1:5087231/1‑100 (MQ=255)
                            ccAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCtt                                                                   >  1:2183165/1‑100 (MQ=255)
                             cAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCttt                                                                  <  1:3413760/100‑1 (MQ=255)
                              aGTAAAATAGTCGCTGGCAAACAGGATGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTAGGTCCCGCTTTc                                                                 >  2:4379060/1‑100 (MQ=255)
                               tctaaaTAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTcc                                                                <  1:1388540/97‑1 (MQ=255)
                               gTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTcc                                                                <  2:4666668/100‑1 (MQ=255)
                               gTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTcn                                                                <  1:4422927/100‑2 (MQ=255)
                               gTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTcc                                                                >  2:5150042/1‑100 (MQ=255)
                               gTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCAAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTcc                                                                >  1:3393697/1‑100 (MQ=255)
                               ctAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGCTACGTCCCGCTTTcc                                                                <  1:4110091/99‑1 (MQ=255)
                                tAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCa                                                               <  1:1231820/100‑1 (MQ=255)
                                tAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCa                                                               <  1:3329044/100‑1 (MQ=255)
                                 aaaaTAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCAt                                                              <  1:2942030/100‑1 (MQ=255)
                                 aaaaTAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGCTACGTCCCGCTTTCCag                                                              >  2:4072890/1‑99 (MQ=255)
                                  aaaTAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATg                                                             <  1:2535287/100‑1 (MQ=255)
                                  aaaTAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATg                                                             <  1:5477264/100‑1 (MQ=255)
                                   aaTAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGc                                                            >  2:1721344/1‑100 (MQ=255)
                                       gTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTaa                                                        >  2:2047663/1‑100 (MQ=255)
                                          gCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATc                                                     <  2:4216029/100‑1 (MQ=255)
                                            tGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAg                                                   >  1:1136240/1‑100 (MQ=255)
                                             ggCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGc                                                  >  1:3256802/1‑100 (MQ=255)
                                             ggCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGc                                                  <  1:3157036/100‑1 (MQ=255)
                                             ggCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGc                                                  <  2:661137/100‑1 (MQ=255)
                                               cAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGcac                                                >  2:4319293/1‑100 (MQ=255)
                                                 aaCAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACag                                              >  1:2926622/1‑100 (MQ=255)
                                                  aCAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACTTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACaga                                             <  1:2547610/100‑1 (MQ=255)
                                                   cAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACagag                                            <  2:3728474/100‑1 (MQ=255)
                                                    aGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAgg                                           <  2:1151709/100‑1 (MQ=255)
                                                     ggggGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGc                                          <  1:3928712/100‑1 (MQ=255)
                                                     gggTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGc                                          >  2:2457935/1‑100 (MQ=255)
                                                        tGATATTTAAAACGACAAATAGCGCCAGCACCGTTAATGTGACTGAAAACCACTGGGTGTTACGTCCCGCTTTCCTTGCGTAAATCAGCACAGAGGCgag                                       <  2:2516696/100‑1 (MQ=255)
                                                        tGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCgag                                       >  1:3224887/1‑100 (MQ=255)
                                                        tGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCgag                                       <  1:4261005/100‑1 (MQ=255)
                                                        tGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCgag                                       <  1:1846203/100‑1 (MQ=255)
                                                        tGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGATACGTGCCGCTTGCCATGCGTAAATCAGCACAGAGGCgag                                       >  2:5249575/1‑100 (MQ=255)
                                                         gATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCgaga                                      >  1:4537526/1‑100 (MQ=255)
                                                            atTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGaaaa                                   >  2:3787022/1‑100 (MQ=255)
                                                            atTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGACCGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGaaaa                                   >  2:4229115/1‑100 (MQ=255)
                                                               tAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAgg                                >  1:3590889/1‑100 (MQ=255)
                                                               tAAAACGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAgg                                >  2:3136447/1‑100 (MQ=255)
                                                                aaaaCGACAAATAGCCCCAGCACCGTTAACGTTACTGCAAACCACCAGGTGTTACGTCCCGCTTTGCATGCGTAACTCAGCACAGAGGCGCGACAACggg                               <  2:5502159/100‑1 (MQ=255)
                                                                aaaaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAggg                               <  2:493067/100‑1 (MQ=255)
                                                                 aaaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGc                              <  1:69150/100‑1 (MQ=255)
                                                                 aaaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGc                              <  1:5517959/100‑1 (MQ=255)
                                                                 aaaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGc                              >  2:1792288/1‑100 (MQ=255)
                                                                  aaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCg                             >  2:2891307/1‑100 (MQ=255)
                                                                  aaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCg                             <  1:1850214/100‑1 (MQ=255)
                                                                  aaCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCg                             <  2:418520/100‑1 (MQ=255)
                                                                  aaCGACAAATAGCCCCAGCACCGTTAACGTCACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCg                             <  2:2409800/100‑1 (MQ=255)
                                                                   aCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGa                            >  2:3241247/1‑100 (MQ=255)
                                                                   aCGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGa                            >  1:2660/1‑100 (MQ=255)
                                                                    cGACAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGaa                           >  2:4958029/1‑100 (MQ=255)
                                                                      ataaaTAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTTCCGCTTTCCATGCGTAAATCGGCACAGAGGCGAGAAAAAGGGCGAAag                         <  2:5476419/98‑1 (MQ=255)
                                                                       cAAATAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAaga                        >  2:3243920/1‑100 (MQ=255)
                                                                         aaTAGCCCCAGCACCGTTAACGTGACTGTAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAagaga                      >  2:4778447/1‑100 (MQ=255)
                                                                         aaTAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAagaga                      >  1:2772277/1‑100 (MQ=255)
                                                                         aaTAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAagaga                      >  1:2236950/1‑100 (MQ=255)
                                                                         aaTAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAagaga                      >  2:2325698/1‑100 (MQ=255)
                                                                          aTAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAagagag                     <  1:4851318/100‑1 (MQ=255)
                                                                          aTAGCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAagagag                     <  1:677998/100‑1 (MQ=255)
                                                                             gCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGagtag                  >  2:3712244/1‑100 (MQ=255)
                                                                             gCCCCAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAAACAGCACAGAGGCGAGAAAAAGGGCGAAAGAGagtag                  >  1:5248837/1‑100 (MQ=255)
                                                                               cccAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGtt                <  2:3995064/100‑1 (MQ=255)
                                                                               cccAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGtt                >  1:1959333/1‑100 (MQ=255)
                                                                                ccAGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTc               <  2:4828505/100‑1 (MQ=255)
                                                                                  aGCACCGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTg             <  1:362821/100‑1 (MQ=255)
                                                                                      ccGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACaa         >  2:4014827/1‑100 (MQ=255)
                                                                                       cGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAAc        <  1:4609169/100‑1 (MQ=255)
                                                                                       cGTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAAc        >  1:5409029/1‑100 (MQ=255)
                                                                                        gTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACg       <  2:1348902/100‑1 (MQ=255)
                                                                                        gTTAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACg       >  1:4356356/1‑100 (MQ=255)
                                                                                          tAACGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACggg     <  1:4970445/100‑1 (MQ=255)
                                                                                           aaCGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGc    <  2:1462086/100‑1 (MQ=255)
                                                                                           aaCGTGACTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGc    <  1:2131274/100‑1 (MQ=255)
                                                                                             cGTGACTGCAAACCACCAGGGGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGGAGTTCTGACAACGGGCGa  >  2:1622798/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGA  >  NC_000913/4599256‑4599448

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: