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,3080GT95.6% 343.4 / 7.4 91intergenic (+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 (2/2);  new base T (52/35);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.13e-01

AAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATT  >  NC_000913/1211161‑1211454
                                                                                                                                                   |                                                                                                                                                  
aaaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACAtttt                                                                                                                                                  <  2:532416/150‑1 (MQ=255)
aaaCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATtttt                                                                                                                                                 >  3:1368052/1‑151 (MQ=255)
  aCCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTtat                                                                                                                                               <  2:1277985/151‑1 (MQ=255)
   ccGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTtat                                                                                                                                               >  3:1149700/1‑150 (MQ=255)
    cGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTtata                                                                                                                                              >  3:130385/1‑150 (MQ=255)
     gATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAAt                                                                                                                                            >  2:25610/1‑151 (MQ=255)
        tACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCa                                                                                                                                          >  2:1152673/1‑150 (MQ=255)
        tACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCa                                                                                                                                          >  2:1612555/1‑150 (MQ=255)
         aCAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAAt                                                                                                                                        <  3:602734/151‑1 (MQ=255)
         aCAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAAt                                                                                                                                        >  2:340375/1‑151 (MQ=255)
         aCAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAAt                                                                                                                                        >  3:427664/1‑151 (MQ=255)
          cAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATg                                                                                                                                       <  2:781925/151‑1 (MQ=255)
            gAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATg                                                                                                                                       >  3:1522996/1‑149 (MQ=255)
             aaaTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGt                                                                                                                                    <  3:219346/151‑1 (MQ=255)
              aaTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGtt                                                                                                                                   >  2:347391/1‑151 (MQ=255)
               aTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGtt                                                                                                                                   >  2:984208/1‑150 (MQ=255)
               aTAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTta                                                                                                                                  <  3:1490611/151‑1 (MQ=255)
                 aaaaTTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTtat                                                                                                                                 >  3:997244/1‑150 (MQ=255)
                    attattTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGaa                                                                                                                             <  2:1062873/151‑1 (MQ=255)
                     ttatttatTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGaa                                                                                                                             >  2:1367613/1‑150 (MQ=255)
                      tatttattAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAAt                                                                                                                           >  2:898972/1‑151 (MQ=255)
                         ttattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCt                                                                                                                        >  2:420748/1‑151 (MQ=255)
                          tattaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCt                                                                                                                        <  3:288393/150‑1 (MQ=255)
                           attaAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTa                                                                                                                       >  2:1380413/1‑150 (MQ=255)
                             taAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAgg                                                                                                                    >  2:602005/1‑151 (MQ=255)
                                    acaTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATg                                                                                                             <  3:487223/151‑1 (MQ=255)
                                     caTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGa                                                                                                            <  2:1260991/151‑1 (MQ=255)
                                      aTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATg                                                                                                             >  3:1289742/1‑149 (MQ=255)
                                      aTTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGa                                                                                                            <  2:595549/150‑1 (MQ=255)
                                        ttAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGtt                                                                                                         >  2:995069/1‑151 (MQ=255)
                                          aaGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtat                                                                                                       >  3:359607/1‑151 (MQ=255)
                                           aGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtat                                                                                                       >  2:380161/1‑150 (MQ=255)
                                           aGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTtat                                                                                                       >  3:637720/1‑150 (MQ=255)
                                              cGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGc                                                                                                   <  2:436044/151‑1 (MQ=255)
                                               gTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCa                                                                                                  <  3:160662/151‑1 (MQ=255)
                                                 aaTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCaa                                                                                                 >  2:1269333/1‑150 (MQ=255)
                                                  aTACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAAt                                                                                               <  2:334364/151‑1 (MQ=255)
                                                    aCCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAAtt                                                                                              >  2:990503/1‑150 (MQ=255)
                                                     cccTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAg                                                                                            <  3:25608/151‑1 (MQ=255)
                                                         aCAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaaa                                                                                        <  3:655742/151‑1 (MQ=255)
                                                         aCAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTaa                                                                                         >  3:1241559/1‑150 (MQ=255)
                                                             ggTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAAtt                                                                                     >  2:253949/1‑150 (MQ=255)
                                                               tAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATg                                                                                  >  3:89782/1‑151 (MQ=255)
                                                                aaaaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGt                                                                                 >  2:191828/1‑151 (MQ=255)
                                                                  aaaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTg                                                                               <  3:185152/151‑1 (MQ=255)
                                                                   aaTTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGc                                                                              >  2:1066949/1‑151 (MQ=255)
                                                                     ttttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCt                                                                             >  3:1057995/1‑150 (MQ=255)
                                                                      tttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTa                                                                            <  3:506388/150‑1 (MQ=255)
                                                                       ttCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATg                                                                          >  3:1504175/1‑151 (MQ=255)
                                                                           ctGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTc                                                                        >  2:1498219/1‑149 (MQ=255)
                                                                           ctGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTc                                                                        <  2:1140596/149‑1 (MQ=255)
                                                                                cTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGt                                                                 >  3:266857/1‑151 (MQ=255)
                                                                                  tAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGtt                                                                >  2:1497915/1‑150 (MQ=255)
                                                                                    aCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACg                                                             <  3:1107292/151‑1 (MQ=255)
                                                                                    aCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACg                                                             <  2:973542/151‑1 (MQ=255)
                                                                                    aCTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACg                                                             >  3:867587/1‑151 (MQ=255)
                                                                                      ttCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAt                                                           >  3:727007/1‑151 (MQ=255)
                                                                                       tCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGa                                                            <  3:1604296/149‑1 (MQ=255)
                                                                                          gCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGAttt                                                         >  2:405641/1‑149 (MQ=255)
                                                                                          gCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTa                                                        <  3:166298/150‑1 (MQ=255)
                                                                                            aaaTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTc                                                     >  2:550411/1‑151 (MQ=255)
                                                                                             aaTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTc                                                     >  3:325770/1‑150 (MQ=255)
                                                                                             aaTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCa                                                    <  3:1000539/151‑1 (MQ=255)
                                                                                              aTGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCAt                                                   <  3:1220818/151‑1 (MQ=255)
                                                                                                    aCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaact                                             <  2:73414/151‑1 (MQ=255)
                                                                                                    aCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaact                                             <  3:542302/151‑1 (MQ=255)
                                                                                                    aCTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATctaact                                             <  2:508279/151‑1 (MQ=255)
                                                                                                       gCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAt                                          >  2:676179/1‑151 (MQ=255)
                                                                                                           tttttATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgctg                                      >  2:996115/1‑151 (MQ=255)
                                                                                                           tttttATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgctg                                      >  2:327325/1‑151 (MQ=255)
                                                                                                            ttttATGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAAtgctgc                                     >  3:1290250/1‑151 (MQ=255)
                                                                                                                aTGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCa                                 <  2:1149794/151‑1 (MQ=255)
                                                                                                                 tGCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGTTCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATACTCCGcc                                  <  3:841771/149‑1 (MQ=255)
                                                                                                                  gCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCat                                >  3:668103/1‑150 (MQ=255)
                                                                                                                  gCTAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCat                                >  2:1213121/1‑150 (MQ=255)
                                                                                                                    tAAAAATGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATg                             >  3:1082097/1‑151 (MQ=255)
                                                                                                                     aaaaaTGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATgg                            >  2:492817/1‑151 (MQ=255)
                                                                                                                       aaaTGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGt                          >  2:923128/1‑151 (MQ=255)
                                                                                                                        aaTGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGtt                         <  2:1628341/151‑1 (MQ=255)
                                                                                                                         aTGGTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTg                        >  2:1234133/1‑151 (MQ=255)
                                                                                                                           ggTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGa                      <  2:594027/151‑1 (MQ=255)
                                                                                                                           ggTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGa                      <  3:551959/151‑1 (MQ=255)
                                                                                                                            gTTATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGAc                     >  3:1574969/1‑151 (MQ=255)
                                                                                                                              tATCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACt                    <  3:984138/150‑1 (MQ=255)
                                                                                                                               aTCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTg                   >  3:1182917/1‑150 (MQ=255)
                                                                                                                                tCAAAACTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGa                  >  3:1385984/1‑150 (MQ=255)
                                                                                                                                     aCTCAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGc            >  2:1490235/1‑151 (MQ=255)
                                                                                                                                        cAAAAACACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTg          <  3:312078/150‑1 (MQ=255)
                                                                                                                                        cAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCt           <  2:649250/149‑1 (MQ=255)
                                                                                                                                            aaCACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg     >  3:272303/1‑151 (MQ=255)
                                                                                                                                            aaCACATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTg     >  2:11154/1‑151 (MQ=255)
                                                                                                                                             acacATTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGa    <  2:87670/151‑1 (MQ=255)
                                                                                                                                               acaTTTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGAtt  <  2:1628032/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AAACCGATTACAGAAATAAAATTATTTATTAAAGTCACATTTAAGACGTAATACCCTACAGGGTAAAAATTTTCTCTGATCTTAACTTCTGCAAATGTTAACTGCTATTTTTATGCTAAAAATGGTTATCAAAACTCAAAAACACATGTTTATAATCAATGAGTTATAGAAATGCTAAGGGCTAATGAGTTATATGCAAATTAGTAAAATTATGTTGCTATGTCAGATAGTTACGATTTAGTCATCTAACTAATGCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATT  >  NC_000913/1211161‑1211454

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