Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,972,674 G→T 7.2% E51D (GAG→GAT rffG → dTDP‑glucose 4,6‑dehydratase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6740GT 7.2% 301.6 / 6.6 90E51D (GAG→GATrffGdTDP‑glucose 4,6‑dehydratase 2
Reads supporting (aligned to +/- strand):  ref base G (38/45);  new base T (5/2);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAC  >  NC_000913/3972525‑3972820
                                                                                                                                                     |                                                                                                                                                  
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgag                                                                                                                                                    >  2:851514/1‑150 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                   >  1:930425/1‑150 (MQ=255)
  aaataTTCTGCTACCAGGTGGGGCCGGGGTTCTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGCGCGACGCGGTGGTAGTGGTCGATAAGCGGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagcg                                                                                                                                                  <  1:1457786/150‑1 (MQ=255)
  aaaaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagcg                                                                                                                                                  <  1:827487/150‑1 (MQ=255)
    aaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCt                                                                                                                                                <  1:802522/150‑1 (MQ=255)
    aaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCt                                                                                                                                                >  1:532647/1‑150 (MQ=255)
       ttCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTg                                                                                                                                             <  2:1082762/150‑1 (MQ=255)
       ttCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTg                                                                                                                                             <  1:1289849/150‑1 (MQ=255)
        tCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGc                                                                                                                                            <  1:1370004/150‑1 (MQ=255)
         cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                           >  1:643860/1‑150 (MQ=255)
          tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCTCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                          <  1:1084616/150‑1 (MQ=255)
          tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                          >  1:1064390/1‑150 (MQ=255)
          tGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGAGGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                          >  1:957446/1‑150 (MQ=255)
           gATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCtt                                                                                                                                         >  1:674096/1‑150 (MQ=255)
            aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                                                                                                        >  1:1160715/1‑150 (MQ=255)
             tAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       >  2:1251818/1‑150 (MQ=255)
              aaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                      <  2:177158/150‑1 (MQ=255)
                          cGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGatat                                                                                                                          >  2:347072/1‑150 (MQ=255)
                             gTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGCCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTg                                                                                                                       >  2:829174/1‑150 (MQ=255)
                              tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGCCGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGc                                                                                                                      >  2:1168142/1‑150 (MQ=255)
                                 aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       <  2:1381956/130‑1 (MQ=255)
                                 aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTg                                                                                                                                       >  1:1381956/1‑130 (MQ=255)
                                 aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGAt                                                                                                                   >  1:466904/1‑150 (MQ=255)
                                   tGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCg                                                                                                                 >  2:732513/1‑150 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCgg                                                                                                                <  1:909951/150‑1 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCgg                                                                                                                >  1:1086787/1‑150 (MQ=255)
                                     gCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCggg                                                                                                               <  2:592240/150‑1 (MQ=255)
                                           cgcTGGTGCGTTATATCATCAACGTAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAAc                                                                                                         <  1:564132/150‑1 (MQ=255)
                                            gcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACt                                                                                                        <  1:851514/150‑1 (MQ=255)
                                              tGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTgg                                                                                                      >  1:1559604/1‑150 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCGCGCTTTGCCTTTGAGAAAGTTGATATGTGCGCTCGGGCAGAACTGGc                                                                                                     >  1:392817/1‑150 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTgag                                                                                                                                     <  1:503933/118‑1 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTgag                                                                                                                                     >  2:503933/1‑118 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                                                                                                     <  2:990998/150‑1 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                                                                                                     >  1:642129/1‑150 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                                                                                                     <  1:1634319/150‑1 (MQ=255)
                                                gTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCa                                                                                                    <  2:930425/150‑1 (MQ=255)
                                                    gTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                <  1:742448/150‑1 (MQ=255)
                                                    gTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                <  2:1086787/150‑1 (MQ=255)
                                                    gTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                <  1:1003528/150‑1 (MQ=255)
                                                       atatCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTAt                                                                                             <  2:1090885/150‑1 (MQ=255)
                                                          tcatcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCa                                                                                          <  2:1160715/150‑1 (MQ=255)
                                                           catcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCAc                                                                                         >  2:6480/1‑150 (MQ=255)
                                                               aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAg                                                                                     <  2:368478/150‑1 (MQ=255)
                                                               aacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACAGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAg                                                                                     <  2:1531500/150‑1 (MQ=255)
                                                                acgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGc                                                                                    >  2:453886/1‑150 (MQ=255)
                                                                 cgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCa                                                                                   <  1:732513/150‑1 (MQ=255)
                                                                   aaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATc                                                                                 >  1:673522/1‑150 (MQ=255)
                                                                   aaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATc                                                                                 >  2:1031615/1‑150 (MQ=255)
                                                                     acgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAg                                                                               <  2:1174945/150‑1 (MQ=255)
                                                                       gaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGcc                                                                             <  1:1031615/150‑1 (MQ=255)
                                                                       gaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGcc                                                                             <  2:868061/150‑1 (MQ=255)
                                                                                ggtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTgtgt                                                                    <  2:643860/150‑1 (MQ=255)
                                                                                ggtggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTgtgt                                                                    <  1:1376533/150‑1 (MQ=255)
                                                                                  tggtAGTGGTCGATAAGCTGACCTACGCCGGAATCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGATCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCa                                                                  <  1:1040594/150‑1 (MQ=255)
                                                                                  tggtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCa                                                                  >  2:966013/1‑150 (MQ=255)
                                                                                   ggtAGTGGTCGATAAGCTGACCTAAGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCAt                                                                 >  2:61814/1‑150 (MQ=255)
                                                                                    gtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGTGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATg                                                                <  1:18135/150‑1 (MQ=255)
                                                                                    gtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATg                                                                <  1:627047/150‑1 (MQ=255)
                                                                                         ggTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCt                                                           >  2:947623/1‑150 (MQ=255)
                                                                                            cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGc                                                        >  1:883037/1‑150 (MQ=255)
                                                                                               tAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGc                                                     <  1:1251818/150‑1 (MQ=255)
                                                                                                aaGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCATAACTGGCACGCGTATTCACTGATCATCAGCCAGACTGTGTCATGCATCTGGCAGcc                                                    <  2:674096/150‑1 (MQ=255)
                                                                                                 aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                   >  1:812523/1‑150 (MQ=255)
                                                                                                 aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                   >  1:1152318/1‑150 (MQ=255)
                                                                                                  gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                   <  2:883249/149‑1 (MQ=255)
                                                                                                  gCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                   >  1:883249/1‑149 (MQ=255)
                                                                                                   cTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGaa                                                 <  2:532647/150‑1 (MQ=255)
                                                                                                     gACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAg                                               <  2:466904/150‑1 (MQ=255)
                                                                                                      aCCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGc                                              <  2:399498/150‑1 (MQ=255)
                                                                                                        cTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCa                                            >  1:1105793/1‑150 (MQ=255)
                                                                                                         tACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCAt                                           <  2:642129/150‑1 (MQ=255)
                                                                                                           cGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGt                                         <  1:1572849/150‑1 (MQ=255)
                                                                                                           cGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGt                                         <  1:1399340/150‑1 (MQ=255)
                                                                                                            gCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGtt                                        <  2:812523/150‑1 (MQ=255)
                                                                                                                  aaCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCAGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGt                                  >  1:16331/1‑150 (MQ=255)
                                                                                                                    ccTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                >  2:1033446/1‑150 (MQ=255)
                                                                                                                    ccTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                >  1:802099/1‑150 (MQ=255)
                                                                                                                              cTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGc                      >  1:108456/1‑150 (MQ=255)
                                                                                                                                  cACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                  <  2:883037/150‑1 (MQ=255)
                                                                                                                                  cACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                  >  1:779466/1‑150 (MQ=255)
                                                                                                                                     ctcaaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAg               >  2:1491525/8‑150 (MQ=255)
                                                                                                                                      ggTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGc              <  2:314611/150‑1 (MQ=255)
                                                                                                                                       gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTGACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCg             <  2:1152318/150‑1 (MQ=255)
                                                                                                                                       caaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGCACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCCTCTGGCAGCCGAAAGCCATGTTGCCCGTTCTATTGACGGCCCGGCcgcg             >  1:1667535/6‑150 (MQ=255)
                                                                                                                                       caaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTg                           >  1:424536/6‑136 (MQ=255)
                                                                                                                                       caaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTg                           <  2:424536/131‑1 (MQ=255)
                                                                                                                                         aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt           >  2:1164116/4‑150 (MQ=255)
                                                                                                                                          gcgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGttt          >  2:991875/1‑150 (MQ=255)
                                                                                                                                          gcgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGttt          <  1:347072/150‑1 (MQ=255)
                                                                                                                                          aggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGttt          >  1:414876/3‑150 (MQ=255)
                                                                                                                                           cgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAACCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTa         >  2:404626/1‑150 (MQ=255)
                                                                                                                                              aaaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTg      <  1:1168142/150‑1 (MQ=255)
                                                                                                                                               aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGTGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACCGCCCGGCAGCGTTTATTGa     >  2:852007/1‑150 (MQ=255)
                                                                                                                                                  cgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAc  <  1:1666029/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAC  >  NC_000913/3972525‑3972820

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