Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,366,855 C→A 7.5% 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 7.5% 242.7 / 6.5 96intergenic (‑82/+36)dcuA/aspAC4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base C (48/40);  new base A (3/5);  total (51/45)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAG  >  NC_000913/4366707‑4367001
                                                                                                                                                    |                                                                                                                                                  
cTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAACGGCACGTCATc                                                                                                                                                    <  1:902602/149‑1 (MQ=255)
       ccccAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTGTTTAACAAGTTGATATTAGATTGTTTTTTTTAAGTTACTGCTCACAAGAAAAAAGGCCCGTCATCTGACGTg                                                                                                                                             >  1:794605/1‑149 (MQ=255)
       ccccAATCTGGCGCCCAAGAAGATCACCAGCAAAACTATGATGAGTTCTACAACTAGCATATAAGCCTACCTTGTTTTTTAACAAGTTGATATTAGATTTTAATATTTAAGTTACTGCTAACAAGAAAAAAGGCACGTAAACTGACCTg                                                                                                                                             >  1:417469/1‑149 (MQ=255)
          cAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATTATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCt                                                                                                                                          <  1:683815/149‑1 (MQ=255)
           aaTCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCtt                                                                                                                                         <  1:610501/149‑1 (MQ=255)
             tCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCtttt                                                                                                                                       <  2:367134/149‑1 (MQ=255)
              cTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGACTTCTACAACTAGCATATTCGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCttttt                                                                                                                                      <  1:565276/149‑1 (MQ=255)
                ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                                    >  2:261092/1‑149 (MQ=255)
                ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTa                                                                                                                                    >  2:344773/1‑149 (MQ=255)
                ggCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCttttttc                                                                                                                                    <  2:247686/149‑2 (MQ=255)
                 gcgcCCAAGAAGATCGCCAGCAAATCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTATTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGGCATCTGACGTGCCTTTTTTAt                                                                                                                                   <  2:683508/149‑1 (MQ=255)
                 gcgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCttttttct                                                                                                                                   <  1:397593/149‑3 (MQ=255)
                 gcgcCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCAGATGACGTGCCttttttct                                                                                                                                   <  2:995874/149‑3 (MQ=255)
                    cccAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTg                                                                                                                                <  2:1018400/149‑1 (MQ=255)
                      cAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGta                                                                                                                              >  1:333510/1‑149 (MQ=255)
                        agaagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtact                                                                                                                            >  1:1024658/1‑149 (MQ=255)
                          aagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtactac                                                                                                                          <  1:192856/149‑1 (MQ=255)
                          aagaTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtactac                                                                                                                          >  2:1184646/1‑149 (MQ=255)
                           ataTCGCCAGCAAAACTATTATGAGTTATACAACTAGCATATTAGCCTTCCATGTTTTTTAACAAGTTGATAATAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTAcc                                                                                                                         <  1:398639/147‑1 (MQ=255)
                             aTCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTTCTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCt                                                                                                                       >  1:76063/1‑149 (MQ=255)
                              tCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCAAGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTg                                                                                                                      >  1:717911/1‑149 (MQ=255)
                                gCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGCCGTGCCTTTTTTATTTGTACTACCCTGTa                                                                                                                    >  2:791660/1‑149 (MQ=255)
                                 ccAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGTCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTAc                                                                                                                   >  1:537330/1‑149 (MQ=255)
                                 ccAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTAc                                                                                                                   >  1:322498/1‑149 (MQ=255)
                                    gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAt                                                                                                                >  1:640132/1‑149 (MQ=255)
                                    gCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATAGTTATTTTTAAGTTACTGCTAACAAGAAAAAAGGCACGTCATCAGACGTGCCTTTTATATATGTACTACCCTGTACGAt                                                                                                                >  2:374885/1‑149 (MQ=255)
                                        aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACATGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                            >  2:94480/1‑149 (MQ=255)
                                        aaCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACt                                                                                                            >  2:474125/1‑149 (MQ=255)
                                         aCTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTg                                                                                                           >  2:683034/1‑149 (MQ=255)
                                          cTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGt                                                                                                          >  1:985795/1‑149 (MQ=255)
                                           tATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCATGTACGATTACTGtt                                                                                                         >  1:1852/1‑149 (MQ=255)
                                               atgaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATATGACGTGCCTTTTTTATTTGTACTACCCTGTACGAGTACTGGTCGCt                                                                                                     >  2:818456/1‑149 (MQ=255)
                                               atgaGTTCTAAAACTAGCATATTAGCCTTCCTTGTTGTTATACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCt                                                                                                     <  2:698076/149‑1 (MQ=255)
                                                 gaGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                   <  1:1092733/149‑1 (MQ=255)
                                                  aGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtc                                                                                                  <  2:1024658/149‑1 (MQ=255)
                                                      cTACAACTAGCATATTAGCCTTGCTAGTTTTTTAACAAGTTGATTTTAGATTGTTCTTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                              <  2:987574/149‑1 (MQ=255)
                                                      cTACAACTAGAATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcatca                                                                                              <  1:5109/149‑1 (MQ=255)
                                                       tACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAg                                                                                             >  1:31106/1‑149 (MQ=255)
                                                         cAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGta                                                                                           <  2:322498/149‑1 (MQ=255)
                                                         cAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGta                                                                                           <  2:275144/149‑1 (MQ=255)
                                                           aCTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTCGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGtata                                                                                         >  1:584960/1‑149 (MQ=255)
                                                             tAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGc                                                                                       >  1:185435/1‑149 (MQ=255)
                                                             tAGCATATTAGCCTTCATTGTTTTTTAACAAGTTGATATTAGATAGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGAAGTGACTTTTTTTTTTGTACTAGCCTTTACTATTAATTTTCTATTTAATCAGtatatc                                                                                       >  1:713499/1‑147 (MQ=255)
                                                                cATATTATCCTTCCTTGTTTTTTAACTAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGtt                                                                                    <  2:270651/149‑1 (MQ=255)
                                                                 atatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGttt                                                                                   <  2:1093800/149‑1 (MQ=255)
                                                                  tatTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGtttt                                                                                  >  2:562314/1‑149 (MQ=255)
                                                                     tAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGAATGTTATATTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCt                                                                               <  1:269313/149‑1 (MQ=255)
                                                                      aGCCTTCCTTGTTTTTTATCAAGTTGATATTAGATTGTTATTTTTAAGTTACTTCTCACAAGAAAACAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTTCCCTGTTCGTTTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                              <  2:417469/149‑1 (MQ=255)
                                                                      aGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtt                                                                              <  1:844532/149‑1 (MQ=255)
                                                                       ggcttccttGTTTTTTTACAAGTTGATTTTTGATTTTTATTTTTTAGTCACTGCTCCAAAGTAAATATGCAGGTCTTCTGACGTGGCTTTTTTATTTGTACTACCCTGTACGATTACTGTTTGCTTTCATCAGTATAGCGTTTTGCttt                                                                             <  1:674479/147‑1 (MQ=255)
                                                                       gCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGGCAGCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCttt                                                                             <  2:210786/149‑1 (MQ=255)
                                                                        ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTTCTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTg                                                                            >  1:976466/1‑149 (MQ=255)
                                                                        ccttccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTAATCTGACGTGCCGTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCtttt                                                                            >  1:559613/1‑148 (MQ=255)
                                                                           tccttGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTaa                                                                         <  2:632160/149‑1 (MQ=255)
                                                                              ttGTTTTATAAAAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGcc                                                                      <  2:698295/149‑1 (MQ=255)
                                                                               tGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTCAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCg                                                                     <  1:1004476/149‑1 (MQ=255)
                                                                                tttttttAACAAGTGGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGAACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCgg                                                                    <  2:558573/148‑1 (MQ=255)
                                                                                 ttttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCggg                                                                   >  1:358317/1‑149 (MQ=255)
                                                                                   ttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                 <  1:1013597/149‑1 (MQ=255)
                                                                                   ttttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTg                                                                 >  2:433930/1‑149 (MQ=255)
                                                                                    tttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGc                                                                <  1:1015325/149‑1 (MQ=255)
                                                                                    tttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGc                                                                >  2:945472/1‑149 (MQ=255)
                                                                                     ttAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCa                                                               >  1:232080/1‑149 (MQ=255)
                                                                                      tAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCTTTTTGCTTTGTAAGCCGGGTGCAt                                                              >  2:784602/1‑149 (MQ=255)
                                                                                          aaGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGa                                                          >  1:160562/1‑149 (MQ=255)
                                                                                             ttGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTc                                                       <  1:1057871/149‑1 (MQ=255)
                                                                                             tagATATTAGATTGATATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCGTATTTATTTGTACTACCCTGTACTATTAATGTTCGATTTCATCAGTATAGCGTTATGCTTTGTAAGCCGGGTGCATCAGaata                                                       >  2:994087/3‑146 (MQ=255)
                                                                                              tGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTGATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCt                                                      >  2:110480/1‑149 (MQ=255)
                                                                                                ttatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                    <  1:148650/148‑1 (MQ=255)
                                                                                                atatTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGt                                                    <  2:640132/149‑1 (MQ=255)
                                                                                                  atTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACTTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGTATCAGATTCTGTAc                                                  <  1:395941/149‑1 (MQ=255)
                                                                                                       aTTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGAGGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATATCGTTTTGCTTTGTAAGCCGGGTGCAGCAGATTCTGTACGgaga                                             >  1:1058220/1‑149 (MQ=255)
                                                                                                           ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                         >  2:152864/1‑149 (MQ=255)
                                                                                                           ttATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTACGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAat                                         >  1:134781/1‑149 (MQ=255)
                                                                                                               ttttAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCg                                     <  1:883034/149‑1 (MQ=255)
                                                                                                                 ttAATTTTCTGCTCACAATAAAAAAGGCACGTCATCTGACGTTCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTc                                   <  1:1164961/149‑1 (MQ=255)
                                                                                                                   aaGTTACTGCTCACAACAAAAAATGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATCTCGTCaa                                 <  2:376578/149‑1 (MQ=255)
                                                                                                                      ttACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGtt                              >  1:316408/1‑149 (MQ=255)
                                                                                                                        aCTGCTCACAAGAAAAACGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTATCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTcc                            <  2:1125354/149‑1 (MQ=255)
                                                                                                                          tGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGc                          >  1:622017/1‑149 (MQ=255)
                                                                                                                            cTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCtt                        <  1:791660/149‑1 (MQ=255)
                                                                                                                              acaaataaaaaaGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTTCGGAGAAAATATGGTCAAGTTCCGCTTCa                      >  1:499734/7‑149 (MQ=39)
                                                                                                                               acaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                   <  1:419004/71‑1 (MQ=255)
                                                                                                                               acaAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCttt                                                                                                   >  2:419004/1‑71 (MQ=255)
                                                                                                                                caAGAAAAATGGGACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTTAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGt                    <  2:794605/149‑1 (MQ=255)
                                                                                                                                 aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                   >  2:528633/1‑149 (MQ=255)
                                                                                                                                 aaGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTc                   >  1:649665/1‑149 (MQ=255)
                                                                                                                                    aaaaaaGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCAGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAAc                >  2:689696/1‑149 (MQ=255)
                                                                                                                                         aGGCACGTCAGATGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTAGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGAATCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCATTCAACAGAcc           >  2:673687/1‑149 (MQ=255)
                                                                                                                                           gCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTACATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACcgc         >  1:456589/1‑149 (MQ=255)
                                                                                                                                             aCGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCATATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGt       >  1:812553/1‑149 (MQ=255)
                                                                                                                                               gTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc     <  2:31106/149‑1 (MQ=255)
                                                                                                                                               gTCATCTGACGTGCCTTTTTTATTTGTACTACACTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTc     >  1:557769/1‑149 (MQ=255)
                                                                                                                                                tCATCTGCCGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTcc    >  1:141801/1‑149 (MQ=255)
                                                                                                                                                 cATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCa   <  2:76063/149‑1 (MQ=255)
                                                                                                                                                  aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg  <  1:562314/149‑1 (MQ=255)
                                                                                                                                                  aTCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGATTTCATCAGTATAGCGTTTTGCTTTGTACGCCGGGTGCATCAGATTCTGTACGGAGTAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAg  >  1:554452/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCATCAGTATAGCGTTTTGCTTTGTAAGCCGGGTGCATCAGATTCTGTACGGAGAAAATATCGTCAAGTTCCGCTTCAGTCAACAGACCGCGTTCCAG  >  NC_000913/4366707‑4367001

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