Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,007,473 G→A 5.7% intergenic (+21/‑37) ygeW → / → ygeX putative carbamoyltransferase/2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,007,4730GA 5.7% 787.7 / 23.8 243intergenic (+21/‑37)ygeW/ygeXputative carbamoyltransferase/2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (119/110);  new base A (6/8);  total (125/118)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 248 total aligned reads displayed.

TCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGA  >  NC_000913/3007325‑3007621
                                                                                                                                                    |                                                                                                                                                    
tCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTgg                                                                                                                                                      <  3:855426/149‑1 (MQ=255)
     gtaattgaagaagCTAGAGGGAAGCCGTATATTATCGCTGCGTTGCTATTGTACCGTAAATACGCCAACCCAGGTGCACTCCTCGAGCAACTGATGAAGGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGc                                                                                                                                                 <  1:1320947/142‑1 (MQ=255)
     acctacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGc                                                                                                                                                 >  1:1453158/1‑149 (MQ=255)
      ccgacAAAGAAGCCAGCTGGAAGCCTTACATCAGCGCCGCGATGATCCTGTCCCGTAAATACGCCAAGCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAAGAACCCCTCGGGCCGATCAGCTGGTCGGcc                                                                                                                                                <  2:233069/146‑1 (MQ=255)
        tacAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTTAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCct                                                                                                                                              >  4:1162214/1‑149 (MQ=255)
            aagaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAt                                                                                                                                          >  1:980117/1‑149 (MQ=255)
            aagaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTAt                                                                                                                                          <  1:1717487/149‑1 (MQ=255)
            aagaagCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcccgaggg                                                                                                                                          <  3:179902/149‑6 (MQ=255)
                 gCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCa                                                                                                                                     <  3:394388/149‑1 (MQ=255)
                    aGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCagag                                                                                                                                  >  1:1420064/1‑149 (MQ=255)
                        ggAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCa                                                                                                                              >  1:268235/1‑149 (MQ=255)
                         gAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaa                                                                                                                             >  2:460420/1‑149 (MQ=255)
                         gAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaa                                                                                                                             >  1:1152577/1‑149 (MQ=255)
                           aGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCaaaa                                                                                                                           <  2:324132/149‑1 (MQ=255)
                             ccTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagg                                                                                                                         <  4:249847/149‑1 (MQ=255)
                              cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCTTCTCTTTCCAGAGGcctaccgga                                                                                                                        <  2:717653/149‑8 (MQ=255)
                              cTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAagga                                                                                                                        <  3:61715/149‑1 (MQ=255)
                               ttATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGACCAAAaggat                                                                                                                       >  4:1103126/1‑149 (MQ=255)
                                tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAGCTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                                      <  1:1665692/149‑1 (MQ=255)
                                tataTCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                                      <  3:332757/149‑1 (MQ=255)
                                   atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatagga                                                                                                                   <  2:268235/149‑1 (MQ=255)
                                   atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggatagga                                                                                                                   <  3:1025249/149‑1 (MQ=255)
                                   atcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTACTTCCAGAGGCCAAAaggatagga                                                                                                                   <  1:717757/149‑1 (MQ=255)
                                    tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAatgataggat                                                                                                                  >  2:1438992/1‑149 (MQ=255)
                                    tcatcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggataggat                                                                                                                  <  4:172370/149‑1 (MQ=255)
                                       tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTcc                                                                                                                                      <  3:75676/126‑1 (MQ=255)
                                       tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTcc                                                                                                                                      >  4:75676/1‑126 (MQ=255)
                                       tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                                      >  4:164877/1‑142 (MQ=255)
                                       tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATg                                                                                                               >  3:1596253/1‑149 (MQ=255)
                                       tcGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGACCAGCTGATCGGcccgag                                                                                                                                            <  3:802834/120‑4 (MQ=255)
                                         gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTTGGGCCGACCAGCTGATCGGcccgagg                                                                                                                                           <  2:823802/119‑5 (MQ=255)
                                         gccgcGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTTGGGCCGACCAGCTGATCGGcccgagg                                                                                                                                           >  1:823802/1‑115 (MQ=255)
                                              gatgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTc                                                                                                                                                    >  3:1468168/1‑105 (MQ=255)
                                               atgatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAATGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGGTAGGATATGTCCGtttt                                                                                                       <  2:201191/149‑1 (MQ=255)
                                                 gatCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTtct                                                                                                     <  4:5108/149‑1 (MQ=255)
                                                   tCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCa                                                                                                   >  1:226437/1‑149 (MQ=255)
                                                     cTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACTCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCGATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCAtt                                                                                                 <  3:1252438/149‑1 (MQ=255)
                                                     cTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCAtt                                                                                                 <  2:1562557/149‑1 (MQ=255)
                                                        tCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGaa                                                                                              >  1:1501246/1‑149 (MQ=255)
                                                               aaTACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGt                                                                                                                                                     >  4:433196/1‑87 (MQ=255)
                                                               aaTACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGt                                                                                                                                                     <  3:433196/87‑1 (MQ=255)
                                                                 tACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTGCGTTTTCTCATTGAAGATTGatat                                                                                     >  1:1592844/1‑149 (MQ=255)
                                                                  aCGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATc                                                                                    >  3:54713/1‑149 (MQ=255)
                                                                    gCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGGTATGTCCGTTTTCTCATTGAAGATTGATATCGc                                                                                  >  4:123005/1‑149 (MQ=255)
                                                                     ccAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGCCCGTTTTCTCATTGAAGATTGATATCgtc                                                                                 <  3:1103126/149‑3 (MQ=255)
                                                                     ccAAACCAGGTGCACTGCTCGAGCAACTACTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGcc                                                                                 <  1:558077/149‑1 (MQ=255)
                                                                       aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                               >  3:1480335/1‑149 (MQ=255)
                                                                       aaaCCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                               <  4:208315/149‑1 (MQ=255)
                                                                            aGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACa                                                                          >  4:888212/1‑149 (MQ=255)
                                                                             ggTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACaa                                                                         <  2:1152577/149‑1 (MQ=255)
                                                                              gTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACaaa                                                                        >  2:160387/1‑149 (MQ=255)
                                                                                gCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAAtt                                                                      >  3:1430364/1‑149 (MQ=255)
                                                                                 cACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAAttt                                                                     <  3:34180/149‑1 (MQ=255)
                                                                                  aCTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAAtttt                                                                    <  2:1361277/149‑1 (MQ=255)
                                                                                        cGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACg                                                              <  3:1010451/149‑1 (MQ=255)
                                                                                         gAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCAtt                                                                                                 >  3:486651/1‑113 (MQ=255)
                                                                                         gAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGa                                                                               <  4:470313/131‑1 (MQ=255)
                                                                                           gCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCg                                                           <  4:1345433/149‑1 (MQ=255)
                                                                                               ctgctgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAAc                                                       >  3:172453/1‑149 (MQ=255)
                                                                                                 gctgAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTGTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACAt                                                     <  4:1462712/149‑1 (MQ=255)
                                                                                                          agCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAACGATAGGATATTTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTGCAACGGCGAAACATCACCGctct                                            <  4:648291/149‑1 (MQ=255)
                                                                                                           gcgcAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCtt                                           >  3:9621/1‑149 (MQ=255)
                                                                                                             gcAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCtttt                                         >  1:1185884/1‑149 (MQ=255)
                                                                                                              cAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTc                                        <  2:914459/149‑1 (MQ=255)
                                                                                                               aaGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGatat                                                                                                                >  1:446887/1‑76 (MQ=255)
                                                                                                                aGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCAGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGc                                      >  3:1668979/1‑149 (MQ=255)
                                                                                                                  aaCGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCaa                                    <  4:1547293/149‑1 (MQ=255)
                                                                                                                   aCGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCaaa                                   <  3:1344234/149‑1 (MQ=255)
                                                                                                                    cgcgTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAg                                  <  2:1592844/149‑1 (MQ=255)
                                                                                                                       gTGAAATAACCCCTCGTGCCGATCAGCTGGTCGGCCCTCTATTTCCAGATGCCAAAAGGATAGGCTATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTTTTTTCGCAAAGCCa                               >  4:82835/1‑149 (MQ=255)
                                                                                                                       gTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGAtt                                                                                          <  2:1337308/90‑1 (MQ=255)
                                                                                                                        tGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAg                              <  2:440254/149‑1 (MQ=255)
                                                                                                                        tGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAg                              >  2:377888/1‑149 (MQ=255)
                                                                                                                        tGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAg                              <  2:602419/149‑1 (MQ=255)
                                                                                                                         gAAATAACCCCTCGGGCCGATCAGCTGGTCAGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAgg                             <  4:1668979/149‑1 (MQ=255)
                                                                                                                          aaaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                                      >  3:1312540/1‑59 (MQ=255)
                                                                                                                          aaaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAaggata                                                                                                                      <  4:1312540/59‑1 (MQ=255)
                                                                                                                          aaaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGc                            >  1:1109055/1‑149 (MQ=255)
                                                                                                                          aaaTAACCCCTCGGGCCGATCAGCTAGTCGGCCCTCTATTTCCAGAGGCCAGAAGGATAGGATATGTCCGTCTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGc                            >  4:1074465/1‑149 (MQ=255)
                                                                                                                           aaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTc                                                                  >  4:108670/1‑110 (MQ=255)
                                                                                                                           aaTAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTc                                                                  <  3:108670/110‑1 (MQ=255)
                                                                                                                             tAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCaa                         >  2:6790/1‑149 (MQ=255)
                                                                                                                              aaCCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCaaa                        >  2:1261088/1‑149 (MQ=255)
                                                                                                                               aCCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAAc                       <  2:1485179/149‑1 (MQ=255)
                                                                                                                                 cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGTCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg                     >  4:1426839/1‑149 (MQ=255)
                                                                                                                                 cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg                     >  4:698828/1‑149 (MQ=255)
                                                                                                                                 cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg                     >  2:514002/1‑149 (MQ=255)
                                                                                                                                 cccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTg                     >  3:1182757/1‑149 (MQ=255)
                                                                                                                                  ccTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTgg                    <  3:666561/149‑1 (MQ=255)
                                                                                                                                   cTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgc                   >  4:1375348/1‑149 (MQ=255)
                                                                                                                                   cTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgc                   >  2:28993/1‑149 (MQ=255)
                                                                                                                                     cGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgc                 <  2:1501246/149‑1 (MQ=255)
                                                                                                                                     cGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgc                 <  3:1162214/149‑1 (MQ=255)
                                                                                                                                      gggCCGACCAGCTGATCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGgcgcg                >  1:1268537/1‑149 (MQ=255)
                                                                                                                                        gCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCAtt                                                                                                 >  4:767285/1‑66 (MQ=255)
                                                                                                                                             tCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTc         >  4:291787/1‑149 (MQ=255)
                                                                                                                                               aGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTcca       <  3:1375348/149‑1 (MQ=255)
                                                                                                                                                gCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTccac      <  1:648243/149‑1 (MQ=255)
                                                                                                                                                    gTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGa  >  3:894964/1‑149 (MQ=255)
                                                                                                                                                    gTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGa  >  1:904485/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGCTACCTACAAAGAAGCCAGCTGGAAGCCTTATATCATCGCCGCGATGATCCTGTCCCGTAAATACGCCAAACCAGGTGCACTGCTCGAGCAACTGCTGAAAGAAGCGCAAGAACGCGTGAAATAACCCCTCGGGCCGATCAGCTGGTCGGCCCTCTATTTCCAGAGGCCAAAAGGATAGGATATGTCCGTTTTCTCATTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGA  >  NC_000913/3007325‑3007621

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