Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,854 T→G 10.1% intergenic (‑81/+37) 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,8540TG10.1% 346.0 / 24.4 127intergenic (‑81/+37)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base T (58/56);  new base G (6/7);  total (64/63)
Fisher's exact test for biased strand distribution p-value = 7.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 128 total aligned reads displayed.

CCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366706‑4367003
                                                                                                                                                      |                                                                                                                                                     
ccTATTCCCCCCATTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATc                                                                                                                                                      <  1:522315/150‑1 (MQ=255)
    ttCCCCAGAATCTGGCGCAGAAGAAGTTCTCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTACTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCGGAc                                                                                                                                                  <  1:101831/150‑1 (MQ=255)
     tCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCAGATGACg                                                                                                                                                 <  2:1006198/150‑1 (MQ=255)
       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTg                                                                                                                                               >  2:1136412/1‑150 (MQ=255)
        ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGc                                                                                                                                              >  1:1277334/1‑150 (MQ=255)
          ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCt                                                                                                                                            <  1:1011369/150‑1 (MQ=255)
          ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCt                                                                                                                                            <  2:1181277/150‑1 (MQ=255)
          ccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCt                                                                                                                                            >  2:1158107/1‑150 (MQ=255)
           cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCtt                                                                                                                                           >  2:1389670/1‑150 (MQ=255)
            aaTCTGGCGCCCAAGCAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTCACAAGTTGATATTAGATTGTTAGTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCCCGTCATCTGACGTGCCttt                                                                                                                                          <  1:849312/150‑1 (MQ=255)
              tCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCAGATGACGTGCCttttt                                                                                                                                        <  2:1240635/150‑1 (MQ=255)
                tGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                                      >  1:1014299/1‑150 (MQ=255)
                 ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCAGATGACGTGCCttttttct                                                                                                                                     <  1:1341064/150‑3 (MQ=255)
                   cgcCCAAGTAGATCTCCAGCAAAACTCTGAGGAGAGCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGGTGATATGAGCTTGTTATTTTTCAGTTACTGCTCACCAG‑‑AAAAAAGGCACGTCAGATGACGTGCCTTTTTTCttg                                                                                                                                   <  2:329733/150‑2 (MQ=37)
                   cgaccAAGAAGATCGGCAGCAAAACTATTTTCAGTTCTACTTCTTTCATATTAGCCTTCCTTGTTTTTTAACAAGTTTATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTAttt                                                                                                                                   <  1:188246/147‑1 (MQ=255)
                    gcCCAAGAAGATCGCCCGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTAtttt                                                                                                                                  >  1:1334959/1‑149 (MQ=255)
                         agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtacta                                                                                                                             >  1:932867/1‑150 (MQ=255)
                             gaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCt                                                                                                                         >  1:518343/1‑150 (MQ=255)
                              aTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTg                                                                                                                        <  1:600792/150‑1 (MQ=255)
                                 gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTAc                                                                                                                     >  2:208322/1‑150 (MQ=255)
                                    aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                  >  1:1642219/1‑150 (MQ=255)
                                        aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                              >  1:391181/1‑150 (MQ=255)
                                        aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                              >  1:1716414/1‑150 (MQ=255)
                                          aCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGt                                                                                                            >  2:1224438/1‑150 (MQ=255)
                                            tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTc                                                                                                          <  1:180334/150‑1 (MQ=255)
                                              tgatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGc                                                                                                        <  1:1662363/150‑1 (MQ=255)
                                               gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCATTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                       >  1:1635081/1‑150 (MQ=255)
                                               gatgaGTTCTACAACTAGCATATTAGCCTTCCTTGATTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTTACGTGCGTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                       >  1:694128/1‑150 (MQ=255)
                                                  gaGTTCTACAACTAGCATATTCGCCTTCCTTGTTTTTTAACAAGTTGATATTCGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtc                                                                                                    <  2:191705/150‑1 (MQ=255)
                                                   aGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                   <  2:466059/150‑1 (MQ=255)
                                                     ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                 <  2:769377/150‑1 (MQ=255)
                                                        tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                              <  1:412239/150‑1 (MQ=255)
                                                        tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                              >  1:492717/1‑150 (MQ=255)
                                                        tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                              <  1:1435813/150‑1 (MQ=255)
                                                        aaCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGt                                                                                              <  2:702889/149‑1 (MQ=255)
                                                         aCAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTTACGTGCCTTTTTTATTTGTACTACCCTGTAAGATTACTGTTCGCTTTAATCAGta                                                                                             >  1:1504646/1‑150 (MQ=255)
                                                          cAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTAAGATTACTGTTCGCTTTCATCAGtat                                                                                            >  1:1301237/1‑150 (MQ=255)
                                                            aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAg                                                                                          <  2:1536048/150‑1 (MQ=255)
                                                             cTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                         >  1:466040/1‑150 (MQ=255)
                                                                 cATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAGGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGttt                                                                                     >  2:221695/1‑150 (MQ=255)
                                                                   tatTAGCCTTCCTTTTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTg                                                                                   <  1:1375394/150‑1 (MQ=255)
                                                                   tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGttttt                                                                                   >  2:1132295/1‑149 (MQ=255)
                                                                   tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTg                                                                                   >  1:16229/1‑150 (MQ=255)
                                                                   tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTg                                                                                   >  2:746266/1‑150 (MQ=255)
                                                                    atTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGc                                                                                  <  1:519285/150‑1 (MQ=255)
                                                                         ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCCCAAGAAAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGGTTTGCgtt                                                                               <  2:271854/150‑1 (MQ=255)
                                                                         ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGATTTGt                                                                             >  2:666008/1‑150 (MQ=255)
                                                                              cttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                        >  1:637197/1‑150 (MQ=255)
                                                                                tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTAATTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATATCGTTTTGCTTTGTACGCCgg                                                                      >  2:81309/1‑150 (MQ=255)
                                                                                 gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGTCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                     <  1:1158107/150‑1 (MQ=255)
                                                                                 gTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                     <  2:1479127/150‑1 (MQ=255)
                                                                                  ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGt                                                                    <  1:67440/150‑1 (MQ=255)
                                                                                  tattttAACAAGTTGATATAAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGAGCCTATTTTATTTGTACTACCCTGTACGATTACTGTTCTCTTTCATCAGTTTAGCGTTTTGCTTTGTAAGCCgggg                                                                    >  2:1067073/3‑149 (MQ=255)
                                                                                    ttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGc                                                                  <  1:924149/150‑1 (MQ=255)
                                                                                     tttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                                 >  2:442227/1‑150 (MQ=255)
                                                                                      ttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCAt                                                                >  1:484832/1‑150 (MQ=255)
                                                                                           aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAt                                                           <  2:391181/150‑1 (MQ=255)
                                                                                               tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                       <  1:1130302/150‑1 (MQ=255)
                                                                                                gATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                      <  2:528811/150‑1 (MQ=255)
                                                                                                 atatTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTa                                                     >  2:602859/1‑150 (MQ=255)
                                                                                                 atatTAGATTGTAATTTTTATGTAACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTTACGTGCCTTTTTTATTTGTACTACCATGTACGATTAATGTTCGCTTTCATCATTATAGCGTATTGATTTGTAACCCGGTTGAATAAGATTCtctc                                                     >  2:246166/1‑147 (MQ=255)
                                                                                                 atatTAGATAGTTATTTTTAAGATACTGCTCACAAG‑‑AAAAAAGGCACGTCATATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTACGCTTTCATCAGAATATCGTTTTGCTTTGTAAGCACGGTGCAAAAGATTCTGTa                                                     >  2:76218/1‑150 (MQ=255)
                                                                                                  tatTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCTGGTGAATCAGATTCTGTaa                                                    >  1:1402712/1‑149 (MQ=255)
                                                                                                    ttAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACgg                                                  <  2:1700922/150‑1 (MQ=255)
                                                                                                    ttAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACgg                                                  <  2:932867/150‑1 (MQ=255)
                                                                                                      aGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGgag                                                >  2:367686/1‑150 (MQ=255)
                                                                                                        aTTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaa                                              >  1:828905/1‑150 (MQ=255)
                                                                                                          tGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaaaa                                            >  2:529555/1‑150 (MQ=255)
                                                                                                               tttttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCGTTTTTTATTTGTACTACCCTGTACGAATACTGTTCGCTATCATCAGTATAGCGTTTTGCTGTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCg                                       >  1:1506938/1‑150 (MQ=255)
                                                                                                                ttttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGt                                      >  1:1678890/1‑150 (MQ=255)
                                                                                                                 tttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                     <  1:1136412/150‑1 (MQ=255)
                                                                                                                  ttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCa                                    <  1:1169981/150‑1 (MQ=255)
                                                                                                                  ttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCa                                    <  1:781304/150‑1 (MQ=255)
                                                                                                                        tACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTcc                              <  1:240826/150‑1 (MQ=255)
                                                                                                                           tGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCt                           <  2:906910/150‑1 (MQ=255)
                                                                                                                                 caAG‑‑ACAAAATGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     <  2:798280/150‑1 (MQ=255)
                                                                                                                                 caAG‑‑AAAAAAGGCACGTCATCTGACGTGCATTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                     <  2:1635081/150‑1 (MQ=255)
                                                                                                                                  aaG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCa                    >  2:1391542/1‑150 (MQ=255)
                                                                                                                                    g‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                  <  1:1068641/150‑1 (MQ=255)
                                                                                                                                    g‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                  <  1:208322/150‑1 (MQ=255)
                                                                                                                                    g‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                  <  2:1014299/150‑1 (MQ=255)
                                                                                                                                    g‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                  >  1:979443/1‑150 (MQ=255)
                                                                                                                                     ataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCaa                   >  1:1495714/3‑150 (MQ=255)
                                                                                                                                        aaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg                <  2:1277334/150‑1 (MQ=255)
                                                                                                                                        aaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAg                >  1:1602518/1‑150 (MQ=255)
                                                                                                                                         aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa               >  1:794124/1‑150 (MQ=255)
                                                                                                                                         aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATAACGTATCGTTCCGCTTCAGTCTACAGa               >  2:335408/1‑150 (MQ=255)
                                                                                                                                           aaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc             >  2:153308/1‑150 (MQ=255)
                                                                                                                                             ggCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc           >  1:1655537/1‑150 (MQ=255)
                                                                                                                                              ggaCGTCATCTGACGTGCATTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg          <  1:1272828/148‑1 (MQ=255)
                                                                                                                                              gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg          <  2:629703/150‑1 (MQ=255)
                                                                                                                                              gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg          <  2:752869/150‑1 (MQ=255)
                                                                                                                                               cACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGCCCGCGt         >  2:862646/1‑150 (MQ=255)
                                                                                                                                                aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGtt        <  1:1389670/150‑1 (MQ=255)
                                                                                                                                                 cGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc       <  2:763269/150‑1 (MQ=255)
                                                                                                                                                 cGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc       >  1:837691/1‑150 (MQ=255)
                                                                                                                                                  gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTcc      <  1:442227/150‑1 (MQ=255)
                                                                                                                                                   tCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCGGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCa     >  2:557053/1‑150 (MQ=255)
                                                                                                                                                    cATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg    <  1:1335882/150‑1 (MQ=255)
                                                                                                                                                      tCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGaa  >  1:387080/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAA  >  NC_000913/4366706‑4367003

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: