Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 9.9% 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,8550CA 9.9% 352.0 / 19.8 127intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (58/56);  new base A (7/6);  total (65/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

Only 100 of 129 total aligned reads displayed.

CCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAC  >  NC_000913/4366706‑4367004
                                                                                                                                                       |                                                                                                                                                     
    ttCCCCAGAATCTGGCGCAGAAGAAGTTCTCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTACTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCGGAc                                                                                                                                                   <  1:101831/150‑1 (MQ=255)
       cccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTg                                                                                                                                                >  2:1136412/1‑150 (MQ=255)
        ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGc                                                                                                                                               >  1:1277334/1‑150 (MQ=255)
        ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGGTCTACACCTAGCAGATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGc                                                                                                                                               >  1:508612/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)
            aaTCTGGCGCCCAAGCAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTCACAAGTTGATATTAGATTGTTAGTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCCCGTCATCTGACGTGCCttt                                                                                                                                           <  1:849312/150‑1 (MQ=255)
            aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCAGATGACGTGCCttt                                                                                                                                           <  2:1173632/150‑1 (MQ=255)
            aaTCTGGCCCCCTAGAAGATCGCCAGCAAAACTATGATGAGTTCTAAAACTAGAATATTAGCCTTCCTTGTTTTTTAACAAGTTGTTATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCttt                                                                                                                                           <  2:245409/150‑1 (MQ=255)
              tCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCAGATGACGTGCCttttt                                                                                                                                         <  2:1240635/150‑1 (MQ=255)
               cTGGCGCCCAAGAAGATCTCCAGCAAAACTATGATGAGTTCTACAACTATCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCtttttt                                                                                                                                        <  2:1344267/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)
                                 gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTAc                                                                                                                      >  2:1068641/1‑150 (MQ=255)
                                 gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTAc                                                                                                                      >  1:763269/1‑150 (MQ=255)
                                    aGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                   >  1:1642219/1‑150 (MQ=255)
                                        aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                               >  1:391181/1‑150 (MQ=255)
                                        aaaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                               <  1:413503/150‑1 (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)
                                                   aGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                    <  2:466059/150‑1 (MQ=255)
                                                   aGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtca                                                                                                    >  1:1377115/1‑150 (MQ=255)
                                                     ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                  <  2:769377/150‑1 (MQ=255)
                                                     ttCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatc                                                                                                  <  2:762796/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‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtat                                                                                             <  2:360734/150‑1 (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)
                                                                   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                                                                         >  2:59911/1‑150 (MQ=255)
                                                                              cttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                         >  1:637197/1‑150 (MQ=255)
                                                                                tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                       >  1:1608318/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)
                                                                                  tattttAACAAGTTGATATAAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGAGCCTATTTTATTTGTACTACCCTGTACGATTACTGTTCTCTTTCATCAGTTTAGCGTTTTGCTTTGTAAGCCgggg                                                                     >  2:1067073/3‑149 (MQ=255)
                                                                                    ttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGc                                                                   <  1:924149/150‑1 (MQ=255)
                                                                                     tttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                                  >  2:442227/1‑150 (MQ=255)
                                                                                       tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATc                                                                >  2:18124/1‑150 (MQ=255)
                                                                                        aaCAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCa                                                               >  1:524618/1‑150 (MQ=255)
                                                                                           aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAt                                                            <  2:391181/150‑1 (MQ=255)
                                                                                               tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                        <  1:1130302/150‑1 (MQ=255)
                                                                                               tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTg                                                        <  2:1495300/150‑1 (MQ=255)
                                                                                                gATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                       <  2:528811/150‑1 (MQ=255)
                                                                                                 atatTAGATAGTTATTTTTAAGATACTGCTCACAAG‑‑AAAAAAGGCACGTCATATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTACGCTTTCATCAGAATATCGTTTTGCTTTGTAAGCACGGTGCAAAAGATTCTGTa                                                      >  2:76218/1‑150 (MQ=255)
                                                                                                  tatTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCTGGTGAATCAGATTCTGTaa                                                     >  1:1402712/1‑149 (MQ=255)
                                                                                                  tatTAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTAc                                                     <  1:1188537/150‑1 (MQ=255)
                                                                                                    ttAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACgg                                                   <  2:932867/150‑1 (MQ=255)
                                                                                                    ttAGATTGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACgg                                                   <  2:1700922/150‑1 (MQ=255)
                                                                                                          tGTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGaaaa                                             >  2:529555/1‑150 (MQ=255)
                                                                                                           gTTATTTTTAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                            >  1:177029/1‑150 (MQ=255)
                                                                                                               tttttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCGTTTTTTATTTGTACTACCCTGTACGAATACTGTTCGCTATCATCAGTATAGCGTTTTGCTGTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCg                                        >  1:1506938/1‑150 (MQ=255)
                                                                                                                ttttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGt                                       >  1:1678890/1‑150 (MQ=255)
                                                                                                                ttttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTAATACCCAGTACGATTACTGTTCGCTTTCATCAGTATAGCGATTTGCTTTGTAAGCCGGGTGCATCAGATTATGTACGGAGAAAATATCGt                                       >  1:356376/1‑150 (MQ=255)
                                                                                                                 tttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                      <  1:1136412/150‑1 (MQ=255)
                                                                                                                  ttAAGTTACTGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCa                                     <  1:781304/150‑1 (MQ=255)
                                                                                                                           tGCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCt                            <  2:906910/150‑1 (MQ=255)
                                                                                                                            gCTCACAAG‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                           >  1:837558/1‑150 (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)
                                                                                                                                         aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGa                >  1:794124/1‑150 (MQ=255)
                                                                                                                                         aaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATAACGTATCGTTCCGCTTCAGTCTACAGa                >  2:335408/1‑150 (MQ=255)
                                                                                                                                           aaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGAcc              >  1:941844/1‑150 (MQ=255)
                                                                                                                                             ggCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc            >  1:1655537/1‑150 (MQ=255)
                                                                                                                                              ggaCGTCATCTGACGTGCATTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg           <  1:1272828/148‑1 (MQ=255)
                                                                                                                                              gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg           <  2:752869/150‑1 (MQ=255)
                                                                                                                                              gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgcg           <  2:629703/150‑1 (MQ=255)
                                                                                                                                               cACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGCCCGCGt          >  2:862646/1‑150 (MQ=255)
                                                                                                                                                aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGtt         <  1:1389670/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‑‑AAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAGAAC  >  NC_000913/4366706‑4367004

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: