Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,211,308 G→T intergenic (+129/‑105) mcrA → / → icdC putative 5‑methylcytosine/5‑hydroxymethylcytosine‑specific restriction nuclease; 5‑methylcytosine DNA binding protein/pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,211,3080GT100.0% 334.4 / NA 86intergenic (+129/‑105)mcrA/icdCputative 5‑methylcytosine/5‑hydroxymethylcytosine‑specific restriction nuclease; 5‑methylcytosine DNA binding protein/pseudogene, isocitrate dehydrogenase C‑terminal gene fragment; idcC' is a 54 codon 3' gene fragment created during e14 prophage insertion
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (48/38);  total (48/38)

AAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTG  >  NC_000913/1211161‑1211455
                                                                                                                                                   |                                                                                                                                                   
aaaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATtttt                                                                                                                                                  >  2:257623/1‑151 (MQ=255)
 aaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTta                                                                                                                                                 <  1:1142321/151‑1 (MQ=255)
  aCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTtat                                                                                                                                                >  2:583156/1‑151 (MQ=255)
    cGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTtata                                                                                                                                               >  2:547713/1‑150 (MQ=255)
     gATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAAt                                                                                                                                             <  1:1050182/151‑1 (MQ=255)
            gAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATg                                                                                                                                        >  1:1193617/1‑149 (MQ=255)
             aaaTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGt                                                                                                                                     >  2:718407/1‑151 (MQ=255)
              aaTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGtt                                                                                                                                    >  1:70320/1‑151 (MQ=255)
                tAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTtat                                                                                                                                  >  2:215946/1‑151 (MQ=255)
                 aaaaTTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTtat                                                                                                                                  >  2:538146/1‑150 (MQ=255)
                  aaaTTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTtat                                                                                                                                  >  2:947814/1‑149 (MQ=255)
                     ttatttatTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGaa                                                                                                                              >  1:684337/1‑150 (MQ=255)
                       atttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAAt                                                                                                                            <  2:969272/150‑1 (MQ=255)
                       atttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATg                                                                                                                           >  1:532521/1‑151 (MQ=255)
                        tttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGc                                                                                                                          >  1:828177/1‑151 (MQ=255)
                         ttattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGc                                                                                                                          >  2:427311/1‑150 (MQ=255)
                         ttattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGc                                                                                                                          <  2:18351/150‑1 (MQ=255)
                          tattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCt                                                                                                                         <  1:813587/150‑1 (MQ=255)
                          tattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCt                                                                                                                         <  2:97665/150‑1 (MQ=255)
                           attaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTaa                                                                                                                       <  2:752275/151‑1 (MQ=255)
                              aaaGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAgg                                                                                                                     <  2:78040/150‑1 (MQ=255)
                                   cacaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAAt                                                                                                               <  2:392144/151‑1 (MQ=255)
                                    acaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTATCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATg                                                                                                              >  1:632038/1‑151 (MQ=255)
                                       tttAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGt                                                                                                           >  2:707633/1‑151 (MQ=255)
                                          aaGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtat                                                                                                        >  2:734844/1‑151 (MQ=255)
                                          aaGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtat                                                                                                        <  1:914618/151‑1 (MQ=255)
                                           aGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtata                                                                                                       <  1:696/151‑1 (MQ=255)
                                           aGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtata                                                                                                       <  1:335074/151‑1 (MQ=255)
                                           aGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtat                                                                                                        >  2:1137964/1‑150 (MQ=255)
                                             aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTatat                                                                                                      <  1:7968/150‑1 (MQ=255)
                                             aCGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATg                                                                                                     >  1:518030/1‑151 (MQ=255)
                                                 aaTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCaa                                                                                                  <  1:378890/150‑1 (MQ=255)
                                                  aTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAAt                                                                                                <  2:389037/151‑1 (MQ=255)
                                                  aTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAAt                                                                                                >  2:882451/1‑151 (MQ=255)
                                                    aCCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTa                                                                                              <  1:881716/151‑1 (MQ=255)
                                                      ccTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGt                                                                                            <  1:554281/151‑1 (MQ=255)
                                                       cTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGt                                                                                            >  1:705336/1‑150 (MQ=255)
                                                         aCAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaaa                                                                                         <  2:650030/151‑1 (MQ=255)
                                                         aCAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaaa                                                                                         <  1:304261/151‑1 (MQ=255)
                                                         aCAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTa                                                                                           <  2:634485/149‑1 (MQ=255)
                                                                 aaaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGtt                                                                                 <  2:103587/151‑1 (MQ=255)
                                                                  aaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                                <  1:638428/151‑1 (MQ=255)
                                                                   aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGc                                                                               <  2:42393/151‑1 (MQ=255)
                                                                       ttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTAt                                                                            <  1:414590/150‑1 (MQ=255)
                                                                           ctGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTc                                                                         >  2:63902/1‑149 (MQ=255)
                                                                             gATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTc                                                                         >  2:195161/1‑147 (MQ=255)
                                                                                  tAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGAt                                                                     <  2:965712/146‑1 (MQ=255)
                                                                                    aCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACg                                                              >  2:180688/1‑151 (MQ=255)
                                                                                     cTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACg                                                              >  1:735302/1‑150 (MQ=255)
                                                                                      ttCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAt                                                            >  1:762035/1‑151 (MQ=255)
                                                                                       tCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGa                                                             <  2:1033358/149‑1 (MQ=255)
                                                                                       tCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAAGGCTAAGGGTTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAt                                                            <  1:886108/150‑1 (MQ=255)
                                                                                         tGCAAATGTTAACTGCTATTTTTATGCTAATAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAtt                                                           >  2:451825/1‑149 (MQ=255)
                                                                                         tGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAttt                                                          >  1:1125733/1‑150 (MQ=255)
                                                                                          gCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAttt                                                          >  2:527019/1‑149 (MQ=255)
                                                                                          gCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAttt                                                          >  1:1070098/1‑149 (MQ=255)
                                                                                            aaaTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTc                                                      >  2:409226/1‑151 (MQ=255)
                                                                                             aaTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCa                                                     <  1:253323/151‑1 (MQ=255)
                                                                                                   aaCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaa                                                <  2:703932/150‑1 (MQ=255)
                                                                                                    aCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaact                                              >  2:783492/1‑151 (MQ=255)
                                                                                                     cTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaact                                              >  1:797599/1‑150 (MQ=255)
                                                                                                        cTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtg                                          >  2:832337/1‑151 (MQ=255)
                                                                                                          aTTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgct                                        <  1:696123/151‑1 (MQ=255)
                                                                                                          aTTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgct                                        >  1:672026/1‑151 (MQ=255)
                                                                                                            ttttATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgctgc                                      >  2:1125844/1‑151 (MQ=255)
                                                                                                                aTGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGcc                                   >  1:846316/1‑150 (MQ=255)
                                                                                                                  gCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCata                                <  2:951755/151‑1 (MQ=255)
                                                                                                                   cTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGGTGCGCCatat                               >  1:745074/1‑151 (MQ=255)
                                                                                                                    tAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATg                              >  2:783890/1‑151 (MQ=255)
                                                                                                                     aaaaaTGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATgg                             <  2:220333/151‑1 (MQ=255)
                                                                                                                       aaaTGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGt                           <  1:1106720/151‑1 (MQ=255)
                                                                                                                       aaaTGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGt                           >  1:519962/1‑151 (MQ=255)
                                                                                                                          tGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTgg                        >  2:830815/1‑151 (MQ=255)
                                                                                                                             ttATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACt                     >  2:975065/1‑151 (MQ=255)
                                                                                                                              tATCTAAACTCAAAAACACATTTTTATAATCAATGAGTTTTAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGATGGACt                     >  1:531485/1‑150 (MQ=255)
                                                                                                                               aTCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTg                    >  2:857589/1‑150 (MQ=255)
                                                                                                                               aTCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGa                   <  1:1104951/151‑1 (MQ=255)
                                                                                                                                tCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGa                   >  1:397717/1‑150 (MQ=255)
                                                                                                                                  aaaaCTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGc                <  1:122877/151‑1 (MQ=255)
                                                                                                                                      cTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCt            >  2:598261/1‑151 (MQ=255)
                                                                                                                                      cTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCt            >  2:236992/1‑151 (MQ=255)
                                                                                                                                       tCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTg           >  2:567595/1‑151 (MQ=255)
                                                                                                                                       tCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTg           >  2:482379/1‑151 (MQ=255)
                                                                                                                                        cAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAATTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCt            <  2:854202/149‑1 (MQ=255)
                                                                                                                                         aaaaaCACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGAc         >  2:16311/1‑151 (MQ=255)
                                                                                                                                           aaaCACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCt       <  1:231034/151‑1 (MQ=255)
                                                                                                                                                caTTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTg  <  2:407168/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTG  >  NC_000913/1211161‑1211455

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