Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,281,717 A→C intergenic (+92/‑259) agaA → / → agaS pseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase/tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,281,7170AC100.0% 274.3 / NA 74intergenic (+92/‑259)agaA/agaSpseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase/tagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (34/40);  total (34/40)

ATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTCCATTGAACTTTCAG  >  NC_000913/3281570‑3281863
                                                                                                                                                   |                                                                                                                                                  
aTAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCAcg                                                                                                                                                 >  2:326426/1‑151 (MQ=255)
  aGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCAcgcg                                                                                                                                               >  2:283748/1‑151 (MQ=255)
  aGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCAcgc                                                                                                                                                <  2:800605/150‑1 (MQ=255)
     gggCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCt                                                                                                                                            <  2:510611/151‑1 (MQ=255)
      ggCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCtt                                                                                                                                           >  1:147033/1‑151 (MQ=255)
       gCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCtt                                                                                                                                           >  2:670416/1‑150 (MQ=255)
          aCATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATcc                                                                                                                                       <  1:294298/151‑1 (MQ=255)
             tgtgCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGCTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGAc                                                                                                                                    >  1:501108/1‑151 (MQ=255)
                  aacaaATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACa                                                                                                                               <  1:739129/151‑1 (MQ=255)
                    caaATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAg                                                                                                                              >  1:619963/1‑150 (MQ=255)
                     aaaTCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGtt                                                                                                                            >  1:799528/1‑151 (MQ=255)
                            gATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                                                     <  1:788754/151‑1 (MQ=255)
                                 aGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAg                                                                                                                <  2:152221/151‑1 (MQ=255)
                                     tCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCg                                                                                                            >  2:680605/1‑151 (MQ=255)
                                        aTTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGAt                                                                                                         >  1:456971/1‑151 (MQ=255)
                                           aGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATaa                                                                                                       >  1:866603/1‑150 (MQ=255)
                                                     tGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgcgc                                                                                            >  2:1118227/1‑151 (MQ=255)
                                                      gATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCACACGCGTCTTATCCGACCTACTGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgcgc                                                                                            <  2:652030/150‑1 (MQ=255)
                                                       aTAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGcgccg                                                                                          <  2:36775/151‑1 (MQ=255)
                                                        tAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGcgccgc                                                                                         >  1:1156299/1‑151 (MQ=255)
                                                            ttGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATcc                                                                                     >  2:166676/1‑151 (MQ=255)
                                                             tGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATcc                                                                                     >  2:978847/1‑150 (MQ=255)
                                                              gCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAg                                                                                   <  1:145859/151‑1 (MQ=255)
                                                              gCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAg                                                                                   <  1:936338/151‑1 (MQ=255)
                                                                 ccTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCaa                                                                                <  1:752564/151‑1 (MQ=255)
                                                                    ttATTGGGACTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATcc                                                                             >  1:610893/1‑151 (MQ=255)
                                                                      aTTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATccc                                                                            <  2:1139133/150‑1 (MQ=255)
                                                                         gggCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCtttt                                                                        <  2:81876/151‑1 (MQ=255)
                                                                         gggCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCttt                                                                         <  1:1063047/150‑1 (MQ=255)
                                                                            ccTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCt                                                                     <  1:1098261/151‑1 (MQ=255)
                                                                             cTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCtt                                                                    >  2:882066/1‑151 (MQ=255)
                                                                              ttCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTc                                                                   >  2:333940/1‑151 (MQ=255)
                                                                                cACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTcc                                                                  <  2:757009/150‑1 (MQ=255)
                                                                                 aCTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCtt                                                                <  1:793648/151‑1 (MQ=255)
                                                                                    tCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCttt                                                               <  1:750613/149‑1 (MQ=255)
                                                                                    tCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAt                                                             >  2:750668/1‑151 (MQ=255)
                                                                                      ccccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtct                                                           <  1:10032/151‑1 (MQ=255)
                                                                                       cccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCCTCCAGCAATCCCTTTTGCTTCCTTTAtctt                                                          >  1:664578/1‑151 (MQ=255)
                                                                                       cccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctt                                                          <  1:742125/151‑1 (MQ=255)
                                                                                       cccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctt                                                          <  1:525547/151‑1 (MQ=255)
                                                                                       cccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctt                                                          >  2:1072494/1‑151 (MQ=255)
                                                                                         cGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctttt                                                        <  1:577514/151‑1 (MQ=255)
                                                                                             aGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttcttt                                                    <  1:1081652/151‑1 (MQ=255)
                                                                                             aGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttct                                                      <  1:1110162/149‑1 (MQ=255)
                                                                                               gCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTtttctttc                                                   >  1:467074/1‑150 (MQ=255)
                                                                                               gCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCa                                                  <  2:657106/151‑1 (MQ=255)
                                                                                                ccTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCa                                                  >  1:630437/1‑150 (MQ=255)
                                                                                                   ttctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGGCGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACg                                               >  1:77120/1‑150 (MQ=255)
                                                                                                   ttctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACg                                               <  2:1141227/150‑1 (MQ=255)
                                                                                                       ttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGAt                                             <  2:280040/148‑1 (MQ=255)
                                                                                                       ttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATca                                           <  1:1079991/150‑1 (MQ=255)
                                                                                                        tttctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATcac                                          >  2:860224/1‑150 (MQ=255)
                                                                                                        tttctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATcac                                          >  2:864510/1‑150 (MQ=255)
                                                                                                          tctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACaa                                        <  2:1000657/150‑1 (MQ=255)
                                                                                                                gTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGt                                 >  1:131408/1‑151 (MQ=255)
                                                                                                                   ttGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGtttt                              >  1:997098/1‑151 (MQ=255)
                                                                                                                   ttGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGtttt                              >  2:1171379/1‑151 (MQ=255)
                                                                                                                    tGATCTGTGCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGttt                               <  1:272106/149‑1 (MQ=255)
                                                                                                                         tgtgCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCt                        >  2:183519/1‑151 (MQ=255)
                                                                                                                          gtgCGGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTAtt                           <  1:342730/147‑1 (MQ=255)
                                                                                                                          gtgCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTAtt                           <  1:981503/147‑1 (MQ=255)
                                                                                                                          gtgCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTc                         <  1:475726/149‑1 (MQ=255)
                                                                                                                            gCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATTCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAGCCCTTTTGCTTCCTTTACCTGTTCTTTCAGCGATCACAAATTTCGTTTTATTTCttt                      <  2:262761/150‑1 (MQ=255)
                                                                                                                            gCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCtttt                     <  2:896175/151‑1 (MQ=255)
                                                                                                                            gCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCtttt                     <  2:981639/151‑1 (MQ=255)
                                                                                                                            gCAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCtttt                     <  1:990289/151‑1 (MQ=255)
                                                                                                                             cAGCGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCttttt                    >  1:213605/1‑151 (MQ=255)
                                                                                                                                cGGTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTtct                 >  1:327002/1‑151 (MQ=255)
                                                                                                                                 ggTGTCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTtctc                >  1:53212/1‑151 (MQ=255)
                                                                                                                                   tgtCGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTcc               >  1:535832/1‑150 (MQ=255)
                                                                                                                                      cGGATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTCCAtt            <  1:1012723/150‑1 (MQ=255)
                                                                                                                                       ggATGCGACGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTCCATTGa          <  1:729965/151‑1 (MQ=255)
                                                                                                                                               cGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTCCATTGAACTTTc    >  1:257736/1‑149 (MQ=255)
                                                                                                                                               cGCTCACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTCCATTGAACTTTCAg  <  2:619933/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTCGTTTTATTTCTTTTTTCTCCATTGAACTTTCAG  >  NC_000913/3281570‑3281863

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