Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,249,316 G→A 53.7% K116K (AAG→AAA)  yqjC → DUF1090 family putative periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,249,3160GA53.7% 10.8 / 221.3 168K116K (AAG→AAA) yqjCDUF1090 family putative periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (30/47);  new base A (38/52);  total (68/100)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.64e-01

Only 100 of 168 total aligned reads displayed.

CGCCAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTGCTGAGT  >  NC_000913/3249218‑3249414
                                                                                                  |                                                                                                  
cGCCAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGc                                                                                                   <  1:2026401/100‑1 (MQ=255)
 gCCAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAACGc                                                                                                   >  2:102430/1‑99 (MQ=255)
 gCCAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCt                                                                                                  <  1:3416260/100‑1 (MQ=255)
  ccAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTg                                                                                                 <  1:1885448/100‑1 (MQ=255)
  ccAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTg                                                                                                 <  1:1705399/100‑1 (MQ=255)
  ccAGCAAGATTTAGCCGAAGCGAAGCAAAAAAGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAAAAGCGCAGGAAGAGCTGAAAAAGCt                                                                                                  >  2:1039780/1‑99 (MQ=255)
   cAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTgg                                                                                                >  2:1853046/1‑100 (MQ=255)
   cAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTgg                                                                                                >  1:2876686/1‑100 (MQ=255)
      cAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAg                                                                                             >  2:3023808/1‑100 (MQ=255)
       aaGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAgc                                                                                            >  1:1133591/1‑100 (MQ=255)
         gATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAgcgc                                                                                          >  1:3037138/1‑100 (MQ=255)
            ttAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGcgcgcg                                                                                       <  1:1838631/100‑1 (MQ=255)
            ttAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGcgcgcg                                                                                       <  1:322667/100‑1 (MQ=255)
               gCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGact                                                                                    >  2:848871/1‑100 (MQ=255)
               gCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGact                                                                                    >  2:1437658/1‑100 (MQ=255)
                ccGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGacta                                                                                   >  1:3737755/1‑100 (MQ=255)
                 cGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGactac                                                                                  >  2:1086001/1‑100 (MQ=255)
                 cGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGactac                                                                                  <  2:68639/100‑1 (MQ=255)
                 cGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGactac                                                                                  <  2:594876/100‑1 (MQ=255)
                  gtggCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGactact                                                                                 >  1:2652674/3‑100 (MQ=255)
                    ggCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTaa                                                                               >  1:2231367/1‑100 (MQ=255)
                      cGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACt                                                                             <  2:2652650/100‑1 (MQ=255)
                       gAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTc                                                                            >  2:2926955/1‑100 (MQ=255)
                          gCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTcaca                                                                         <  2:493610/100‑1 (MQ=255)
                           cAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACaa                                                                        <  2:2785219/100‑1 (MQ=255)
                           cAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACaa                                                                        <  2:2785234/100‑1 (MQ=255)
                            aaaaaGGCGATGCCGATAAGATTGCCAACCGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAAt                                                                       >  2:1416955/1‑100 (MQ=255)
                            aaaaaGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAAt                                                                       >  1:491020/1‑100 (MQ=255)
                               aaGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAg                                                                     >  2:2825673/1‑99 (MQ=255)
                                aGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTc                                                                   <  2:622135/100‑1 (MQ=255)
                                aGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTc                                                                   <  2:2895561/100‑1 (MQ=255)
                                aGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTc                                                                   <  1:3357036/100‑1 (MQ=255)
                                  gCGATGCCGATAAGGTTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCac                                                                 <  2:3257747/100‑1 (MQ=255)
                                  gCGATGCAGATAAGATTGCCAAACGCGAACGGAAACTTGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCa                                                                  >  2:2087450/1‑99 (MQ=255)
                                   cGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCact                                                                <  1:2111903/100‑1 (MQ=255)
                                   cGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGAAAGCGAGCGACTACTAACTCACAATAGTCac                                                                 >  1:750319/1‑99 (MQ=255)
                                     aTGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCactac                                                              <  1:2834067/100‑1 (MQ=255)
                                     aTGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCacta                                                               <  1:3170839/99‑1 (MQ=255)
                                      tGCCGATCAGATTGCCAAACGCGAACGAAAACTTGCAGAAGCGCAGGAAGAGCTGAAAAAACCGGAAGCGCGCGACTACTAACTCACAATAGTCact                                                                >  2:3143362/1‑97 (MQ=255)
                                      tGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCactact                                                             >  1:2533076/1‑100 (MQ=255)
                                       gCCGATAAGATTGCCAAACGCGAACGGAAACTTGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACtactt                                                            <  1:2027956/100‑1 (MQ=255)
                                         cGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACtacttac                                                          >  1:1925111/1‑100 (MQ=255)
                                         cGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACtacttac                                                          >  2:685764/1‑100 (MQ=255)
                                         cGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACtacttac                                                          >  1:869327/1‑100 (MQ=255)
                                              aGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCAcc                                                     <  2:3466064/100‑1 (MQ=255)
                                              aGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCAcc                                                     <  1:2222284/100‑1 (MQ=255)
                                               gATTGCCCAACGCAAACAGAAAATGGCAGAAGAGCAGGAAGAGTTGACAAAGATGGAAGCGCGCGACTACTAAGAAAAAATAATAACTACTTACTCAcat                                                    >  2:3322462/1‑98 (MQ=255)
                                               gATTGCCAAACGCGAACGGAAACTGGCAGAAGGGCAGGATGAACTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCt                                                    <  2:2746693/100‑1 (MQ=255)
                                               gATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCt                                                    >  2:1142000/1‑100 (MQ=255)
                                                 ttGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTgg                                                  <  2:1100469/100‑1 (MQ=255)
                                                 ttGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTgg                                                  <  2:2520193/100‑1 (MQ=255)
                                                  tGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGga                                                 <  2:2475443/100‑1 (MQ=255)
                                                   gCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGgag                                                >  1:926151/1‑100 (MQ=255)
                                                   gCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGgag                                                >  2:189647/1‑100 (MQ=255)
                                                   gCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGgag                                                <  2:491013/100‑1 (MQ=255)
                                                    ccAAACGCGAACGGAAACTGGCAGAGGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGgaga                                               >  1:1035700/1‑100 (MQ=255)
                                                    ccAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGgaga                                               >  1:361930/1‑100 (MQ=255)
                                                     cAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGaa                                              >  1:582225/1‑100 (MQ=255)
                                                       aaCGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCTCAATAGTCACTACTTACTCACCTGGAGaan                                             <  2:3453156/99‑2 (MQ=255)
                                                        aCGCGAACGGAAACTGGCAGAAGCGCAGGAAGTGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAAc                                           <  1:2454180/100‑1 (MQ=255)
                                                         cgcgAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAACCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAAc                                           >  2:1224018/1‑99 (MQ=255)
                                                          gcgAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTa                                         <  2:1046314/100‑1 (MQ=255)
                                                           cgAACGGAATCTGGCAGAAGCGTAGTAAGAGCTGAAAAAGCTTGAAGCGCGCGACTACTAACTCACAATAGTCACTACTGACTCACCTGGAGAAAACTa                                         <  1:3139183/99‑1 (MQ=255)
                                                           cgAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTAt                                        <  1:139543/100‑1 (MQ=255)
                                                           cgAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTAACAATAGTCACTACTTACTCACATGGAGAAAACTAt                                        >  1:1570361/1‑100 (MQ=255)
                                                              aCGGAAACTGGCAGAATCGCAGGAAAAGCTTAAACCACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTATTCACCTGGAGAAAACTATGTc                                     <  2:3148370/100‑1 (MQ=255)
                                                               cGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCg                                    <  2:1385968/100‑1 (MQ=255)
                                                                 gAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGaa                                  <  2:3393481/100‑1 (MQ=255)
                                                                  aaaCTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGaaa                                 <  1:1162643/100‑1 (MQ=255)
                                                                   aaCTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAg                                <  1:1789149/100‑1 (MQ=255)
                                                                   aaCTGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAg                                <  1:535951/100‑1 (MQ=255)
                                                                    aCTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTGCTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGa                               <  2:1384576/100‑1 (MQ=255)
                                                                     caggCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGaa                              >  1:2639744/3‑100 (MQ=255)
                                                                      tGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTTGAAGCGCTCGACTACTAACTCATATTAGTCACTACTTACTCACCTGTAGTAAACTATGTCGAAAGAac                             <  2:1116135/100‑1 (MQ=255)
                                                                      tGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAac                             <  1:2502718/100‑1 (MQ=255)
                                                                      tGGCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAac                             <  1:1415836/100‑1 (MQ=255)
                                                                        gCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGTAGAAAACTATGTCGAAAGAacac                           <  2:561388/100‑1 (MQ=255)
                                                                        gCAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAacac                           <  1:1212722/100‑1 (MQ=255)
                                                                         cAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAacaca                          >  1:2627131/1‑100 (MQ=255)
                                                                         cAGAAGCGCAGGAAGAGCTGAAAAAACTGGAAGGGCGCGACTACTAACTCACAATAGTCACTACTTACTCCCCTGGAGAAAACTATGTCGAAAGAacaca                          >  2:2347571/1‑100 (MQ=255)
                                                                          aGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAacacac                         <  1:1664087/100‑1 (MQ=255)
                                                                            aaGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTa                       <  1:3651091/100‑1 (MQ=255)
                                                                             aGCGCAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTAc                      <  2:3755421/100‑1 (MQ=255)
                                                                              gcgcAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTa                       >  2:1510033/1‑98 (MQ=255)
                                                                                gcAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGa                   <  1:2392964/100‑1 (MQ=255)
                                                                                gcAGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGa                   >  1:2826748/1‑100 (MQ=255)
                                                                                  aGGAAGAGCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCATCTGGAGAAAACTATGTCGAAAGAACACACTACGGAAc                 <  2:2086000/100‑1 (MQ=255)
                                                                                   ggAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACa                >  1:130450/1‑100 (MQ=255)
                                                                                   ggAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACa                <  2:1276431/100‑1 (MQ=255)
                                                                                   ggAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACa                <  1:3058399/100‑1 (MQ=255)
                                                                                       gagCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTg            <  2:130446/100‑1 (MQ=255)
                                                                                        agCTTAAAAAACTGGAAGCGCGCGACTACTAACTCACAATACTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTg            <  2:917487/99‑1 (MQ=255)
                                                                                         gCTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCg          <  2:175029/100‑1 (MQ=255)
                                                                                          cTGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGt         >  1:1080137/1‑100 (MQ=255)
                                                                                           tGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTg        <  2:648966/100‑1 (MQ=255)
                                                                                           tGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTg        >  2:2022344/1‑100 (MQ=255)
                                                                                           tGAAAAAACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAAAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTg        >  1:2435474/1‑100 (MQ=255)
                                                                                              aaaaGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTGCTg     >  1:3171157/1‑100 (MQ=255)
                                                                                              aaaaACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTGCTg     >  2:3411977/1‑100 (MQ=255)
                                                                                              aaaaACTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTGCTg     >  1:3502214/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
CGCCAGCAAGATTTAGCCGAGGCGAAGCAAAAAGGCGATGCCGATAAGATTGCCAAACGCGAACGGAAACTGGCAGAAGCGCAGGAAGAGCTGAAAAAGCTGGAAGCGCGCGACTACTAACTCACAATAGTCACTACTTACTCACCTGGAGAAAACTATGTCGAAAGAACACACTACGGAACATCTGCGTGCTGAGT  >  NC_000913/3249218‑3249414

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