Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 10.7% intergenic (‑82/+36) dcuA ← / ← aspA C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,366,8550CA10.7% 528.0 / 39.4 157intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (73/67);  new base A (7/10);  total (80/77)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 165 total aligned reads displayed.

CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366707‑4367003
                                                                                                                                                     |                                                                                                                                                    
cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                      <  4:532819/149‑1 (MQ=255)
cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                      >  3:295335/1‑149 (MQ=255)
   ttCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGa                                                                                                                                                   >  3:1494019/1‑149 (MQ=255)
      cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                <  1:1163178/149‑1 (MQ=255)
      cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGt                                                                                                                                                >  3:63217/1‑149 (MQ=255)
       ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  3:724420/1‑149 (MQ=255)
        cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              >  2:1380662/1‑149 (MQ=255)
        cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              >  4:1194824/1‑149 (MQ=255)
        cccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGAATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              <  2:1551324/149‑1 (MQ=255)
         ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGcc                                                                                                                                             <  2:588464/149‑1 (MQ=255)
           aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                           <  2:723509/149‑1 (MQ=255)
            aTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCttt                                                                                                                                          <  2:1351940/149‑1 (MQ=255)
               tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCt                                                                                                                                                     >  2:797439/1‑135 (MQ=255)
               tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCt                                                                                                                                                     <  1:797439/135‑1 (MQ=255)
                ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAGAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                                      <  3:1372942/149‑1 (MQ=255)
                  cgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtt                                                                                                                                    <  4:63217/149‑1 (MQ=255)
                   gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAATAAAAAAGGCACGTCAGATGACGTGCCTTTTTTCtt                                                                                                                                    <  4:1517993/149‑1 (MQ=38)
                   gcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   <  1:916919/149‑1 (MQ=255)
                     cgctgtagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGt                                                                                                                                 <  4:352632/143‑1 (MQ=255)
                     ccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAGAGGCACGTCATCTGACGTGCCTTTTTTATTTGt                                                                                                                                 <  3:1066277/149‑1 (MQ=255)
                      cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                >  3:540323/1‑149 (MQ=255)
                      cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                                <  1:1281141/149‑1 (MQ=255)
                       aagaagATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                               >  1:1024808/1‑149 (MQ=255)
                         gaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtacta                                                                                                                             >  4:885625/1‑149 (MQ=255)
                           agaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                           >  1:1423571/1‑149 (MQ=255)
                                    gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                     >  1:907244/1‑130 (MQ=255)
                                    gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAt                                                                                                                                     <  2:907244/130‑1 (MQ=255)
                                    gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                  <  3:407746/149‑1 (MQ=255)
                                           tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtttttt                                                                                                                                       <  2:601016/121‑1 (MQ=255)
                                           tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtttttt                                                                                                                                       >  1:601016/1‑121 (MQ=255)
                                               atgaGTTCTACAACTAGCATCTTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                          <  4:1174250/146‑1 (MQ=255)
                                               atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTGAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                          >  3:1174250/1‑146 (MQ=255)
                                                 gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAACCTGTACGATTACTGTTCGCttt                                                                                                     >  1:1524327/1‑149 (MQ=255)
                                                   gTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   <  3:566135/149‑1 (MQ=255)
                                                    ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   <  3:722733/116‑1 (MQ=255)
                                                    ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   >  4:722733/1‑116 (MQ=255)
                                                    ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat                                                                                                  >  1:1311811/1‑149 (MQ=255)
                                                      cTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                <  4:724420/149‑1 (MQ=255)
                                                       tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                               >  2:121535/1‑149 (MQ=255)
                                                       tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                               <  1:531507/149‑1 (MQ=255)
                                                          aaCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtat                                                                                            >  2:1190574/1‑149 (MQ=255)
                                                           aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtttt                                                                                                                                         <  1:677027/103‑1 (MQ=255)
                                                             tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  1:827281/1‑149 (MQ=255)
                                                             tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  4:350959/1‑149 (MQ=255)
                                                               gCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGt                                                                                       >  2:698498/1‑149 (MQ=255)
                                                                 atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGttt                                                                                     <  1:222367/149‑1 (MQ=255)
                                                                  tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGtttt                                                                                    >  1:1687044/1‑149 (MQ=255)
                                                                    ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                                  >  1:791389/1‑149 (MQ=255)
                                                                    ttAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                                  >  3:730429/1‑149 (MQ=255)
                                                                     tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                                 <  2:722375/149‑1 (MQ=255)
                                                                       gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                               <  1:1351651/149‑1 (MQ=255)
                                                                           tccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTaa                                                                           <  3:1688776/149‑1 (MQ=255)
                                                                              ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                               >  2:493574/1‑94 (MQ=255)
                                                                              ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtac                                                                                                                               <  1:493574/94‑1 (MQ=255)
                                                                              ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCTTCAGTATAGCGTTTTGCTTTGTAAGcc                                                                        >  1:711355/1‑149 (MQ=255)
                                                                              ttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                        <  2:637051/149‑1 (MQ=255)
                                                                               tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCg                                                                       >  2:886750/1‑149 (MQ=255)
                                                                                      tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                      <  2:15388/63‑1 (MQ=255)
                                                                                      tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATc                                                                                                                                                      >  1:15388/1‑63 (MQ=255)
                                                                                        aCAAGTTGAGATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTGATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTAGTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCa                                                              >  4:489206/1‑149 (MQ=255)
                                                                                          aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGa                                                            <  4:523695/149‑1 (MQ=255)
                                                                                            gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                <  3:845876/143‑1 (MQ=255)
                                                                                            gTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACATCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                >  4:845876/1‑143 (MQ=255)
                                                                                             ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGTTGCATCAGATTc                                                         <  4:852335/149‑1 (MQ=255)
                                                                                              tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCt                                                        <  2:1639699/149‑1 (MQ=255)
                                                                                                atatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                      <  4:1310478/149‑1 (MQ=255)
                                                                                                 tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                 >  1:937566/1‑105 (MQ=255)
                                                                                                 tatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                     <  2:827281/149‑1 (MQ=255)
                                                                                                  atTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                                    <  3:885625/149‑1 (MQ=255)
                                                                                                   ttAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                <  2:286673/104‑1 (MQ=255)
                                                                                                   ttAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                                >  1:286673/1‑104 (MQ=255)
                                                                                                   ttAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACg                                                   >  3:1505975/1‑149 (MQ=255)
                                                                                                     aGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                      >  1:443539/1‑128 (MQ=255)
                                                                                                     aGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGga                                                 <  3:1608252/149‑1 (MQ=255)
                                                                                                      gATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACg                                                   <  4:1512152/146‑1 (MQ=255)
                                                                                                      gATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgag                                                <  2:1307417/149‑1 (MQ=255)
                                                                                                       aTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgaga                                               >  1:751181/1‑149 (MQ=255)
                                                                                                       aTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgaga                                               >  4:1283064/1‑149 (MQ=255)
                                                                                                           ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                           >  2:664502/1‑149 (MQ=255)
                                                                                                              tttttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATc                                        <  2:791389/149‑1 (MQ=255)
                                                                                                                 ttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                     <  3:478969/149‑1 (MQ=255)
                                                                                                                 ttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                     <  3:1283064/149‑1 (MQ=255)
                                                                                                                      ttACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                                >  1:588468/1‑149 (MQ=255)
                                                                                                                        aCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTcc                              <  1:1190574/149‑1 (MQ=255)
                                                                                                                             tCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTc                         >  3:1030656/1‑149 (MQ=255)
                                                                                                                               acaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCaa                                   <  4:931616/137‑1 (MQ=255)
                                                                                                                               acaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCaa                                   >  3:931616/1‑137 (MQ=255)
                                                                                                                               acaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCa                        >  4:1630601/1‑148 (MQ=255)
                                                                                                                                caAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                      >  1:107163/1‑149 (MQ=255)
                                                                                                                                  aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     >  2:1450371/1‑149 (MQ=255)
                                                                                                                                        aaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa               <  3:1542129/149‑1 (MQ=255)
                                                                                                                                        aaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCg                             >  3:1067175/1‑135 (MQ=255)
                                                                                                                                         aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAc              <  2:184951/149‑1 (MQ=255)
                                                                                                                                          aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc             >  4:1118590/1‑149 (MQ=255)
                                                                                                                                          aGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc             >  3:849747/1‑149 (MQ=255)
                                                                                                                                          aGGCACGTCAGATGACGTGCCTTTTTTATCTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     >  3:943828/1‑141 (MQ=255)
                                                                                                                                            gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGGGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  3:1244516/1‑149 (MQ=255)
                                                                                                                                                   aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg    <  1:494448/149‑1 (MQ=255)
                                                                                                                                                     cTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa  <  4:144845/149‑1 (MQ=255)
                                                                                                                                                     cTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa  <  2:926702/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366707‑4367003

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