Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,527,466 T→G 5.2% Y311D (TAC→GAC)  PP_4860 → mannose‑6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,527,4660TG 5.2% 256.5 / 4.6 88Y311D (TAC→GAC) PP_4860mannose‑6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base G (3/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGCGTGACACCCCGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCAT  >  NC_002947/5527317‑5527610
                                                                                                                                                       |                                                                                                                                                
tgctgcGTGACACCCCGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGt                                                                                                                                                  >  2:1398629/1‑150 (MQ=255)
        gACACCCCGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCgg                                                                                                                                          <  2:158954/150‑1 (MQ=255)
         gtcatccctgctcaccgcgccgtcaccccctcggcactcataccagccggcgtcacgcaggaaccgcgctgaaagggcgcccaactgccccgcgagctttgcttcactgccttcacgcaatgtaagcgccCGCAGGCCGAGTACCTGCgg                                                                                                                                          >  1:964658‑M1/131‑150 (MQ=255)
            ccccGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGgcgc                                                                                                                                      <  1:1090982/150‑1 (MQ=255)
             cccGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCt                                                                                                                                     <  1:561061/150‑1 (MQ=255)
               cGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTa                                                                                                                                   >  2:771510/1‑150 (MQ=255)
                     aTGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGc                                                                                                                             <  1:865102/150‑1 (MQ=255)
                       ccgcgccgtcaccccctcggcactcataccagccggcgtcacgcaggaaccgcgctgaaagggcgcccaactgccccgcgagctttgcttcactgccttcacgcaatgtaagcgccCGCAGGCCGAGTACCTGCg                                                                                                                                           >  2:865877‑M1/117‑135 (MQ=255)
                       ccgcgccgtcaccccctcggcactcataccagccggcgtcacgcaggaaccgcgctgaaagggcgcccaactgccccgcgagctttgcttcactgccttcacgcaatgtaagcgccCGCAGGCCGAGTACCTGCg                                                                                                                                           <  1:865877‑M1/19‑1 (MQ=255)
                         gTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGt                                                                                                                           >  2:1580465/1‑148 (MQ=255)
                         gTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGt                                                                                                                           <  1:1580465/148‑1 (MQ=255)
                            aaTAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAgg                                                                                                                      <  1:1702045/150‑1 (MQ=255)
                             aTAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAgg                                                                                                                      <  2:453198/149‑1 (MQ=255)
                             aTAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAgg                                                                                                                      >  1:453198/1‑149 (MQ=255)
                               aGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGTGGCACGGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAg                                                                                                                   <  1:118848/150‑1 (MQ=255)
                                   cGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGc                                                                                                                     <  1:1636302/144‑1 (MQ=255)
                                   cGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGc                                                                                                                     >  2:1636302/1‑144 (MQ=255)
                                   cGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGaa                                                                                                               <  1:713031/150‑1 (MQ=255)
                                    gggCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAg                                                                                                              >  1:739734/1‑150 (MQ=255)
                                     ggCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGc                                                                                                             <  1:1234541/150‑1 (MQ=255)
                                        tttCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCaaa                                                                                                          >  2:1203865/1‑150 (MQ=255)
                                               tttGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTcgcg                                                                                                   <  2:1272467/150‑1 (MQ=255)
                                                ttGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCgg                                                                                                  <  2:1730851/150‑1 (MQ=255)
                                                 tGCCGAGCAATTCGGGGTGAAGGATTCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCggg                                                                                                 <  1:8023/150‑1 (MQ=255)
                                                  gCCGAGCAATTCGGGGTGAACGATGCGGCCCTAATGGCAAGCTGTAATTTGGATTGCCGTGTGCTTGATGACTCTCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCgggg                                                                                                <  1:1398629/150‑1 (MQ=255)
                                                         aaTTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGATCAAATTCGAGCAAAGCTCGCGGGGCAGTTgg                                                                                         <  1:1021371/150‑1 (MQ=255)
                                                          aTTCGGGGTGAAGGATGCGTCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCGGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTggg                                                                                        <  1:1688589/150‑1 (MQ=255)
                                                               gggTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGTCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGccc                                                                                   >  1:425094/1‑150 (MQ=255)
                                                               gggTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGccc                                                                                   <  2:739734/150‑1 (MQ=255)
                                                                    aaGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCa                                                                              >  2:361041/1‑150 (MQ=255)
                                                                     aGGATGCGGCCGTACTGGCCATGATGAATGTGGCTGCCGGTGTCTTTGAGACCGCCCTGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAg                                                                             <  1:1398452/150‑1 (MQ=255)
                                                                      ggATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAgc                                                                            >  2:1762482/1‑150 (MQ=255)
                                                                       gATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGGATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAgcg                                                                           >  1:1267931/1‑150 (MQ=255)
                                                                        gcatggCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAgcgc                                                                          >  2:1334196/5‑150 (MQ=255)
                                                                          gCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCgg                                                                        <  1:1117852/150‑1 (MQ=255)
                                                                            ggCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGtt                                                                      >  1:684667/1‑150 (MQ=255)
                                                                                 tACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGTCTCCCTGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGc                                                                 <  2:425094/150‑1 (MQ=255)
                                                                                   cTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGt                                                               <  1:1334196/150‑1 (MQ=255)
                                                                                     ggCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGa                                                             <  2:684667/150‑1 (MQ=255)
                                                                                      gCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCgg                                                                                                                                          <  1:1883156/72‑1 (MQ=255)
                                                                                      gCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCgg                                                                                                                                          >  2:1883156/1‑72 (MQ=255)
                                                                                      gCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGgcgc                                                                                                                                      >  2:487905/1‑76 (MQ=255)
                                                                                      gCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGgcgc                                                                                                                                      <  1:487905/76‑1 (MQ=255)
                                                                                       ccATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACg                                                           >  1:1390629/1‑150 (MQ=255)
                                                                                        cATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGc                                                          <  1:361041/150‑1 (MQ=255)
                                                                                             tGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGc                                                                                               <  1:397407/108‑1 (MQ=255)
                                                                                             tGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGc                                                                                               >  2:397407/1‑108 (MQ=255)
                                                                                              gAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTg                                                    >  2:578040/1‑150 (MQ=255)
                                                                                                aaTGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTgg                                                   >  2:1830827/1‑150 (MQ=255)
                                                                                                aaTGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGCACAGATGGGCGCCCTGTCATCGCGGTTCCTGCGTGAAGAAGGCTgg                                                   >  2:1287856/1‑150 (MQ=255)
                                                                                                  tgtgGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTa                                                 <  1:771510/150‑1 (MQ=255)
                                                                                                    tgGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGgc                                                                                       <  2:1090558/110‑1 (MQ=255)
                                                                                                    tgGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGGGAAGCAAAGCTCGCGGGGCAGTTGGgc                                                                                       >  1:1090558/1‑110 (MQ=255)
                                                                                                      gATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGt                                                               <  1:167141/132‑1 (MQ=255)
                                                                                                      gATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGt                                                               >  2:167141/1‑132 (MQ=255)
                                                                                                        tGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGCATGAGTg                                           >  2:894990/1‑150 (MQ=255)
                                                                                                           ccGGGTGCTTGATGCCTCCCCGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCg                                        >  2:320968/1‑150 (MQ=255)
                                                                                                              ggTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAgg                                     >  1:1694298/1‑150 (MQ=255)
                                                                                                                tGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAgggg                                   >  2:1891149/1‑150 (MQ=255)
                                                                                                                  cTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGt                                 >  1:1349307/1‑150 (MQ=255)
                                                                                                                     gccttcacgcaatgtaagcgccCGCAGGTACTCGGCCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGa                                              >  1:1646325‑M1/23‑134 (MQ=255)
                                                                                                                     gccttcacgcaatgtaagcgccCGCAGGTACTCGGCCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGa                                              <  2:1646325‑M1/112‑1 (MQ=255)
                                                                                                                          cTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGgcgc                         <  2:1694298/150‑1 (MQ=255)
                                                                                                                           tggcAGTGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGcgcg                        >  2:520133/4‑150 (MQ=255)
                                                                                                                           tCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGc                                          >  1:342109/1‑132 (MQ=255)
                                                                                                                           tCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGc                                          <  2:342109/132‑1 (MQ=255)
                                                                                                                           tCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAgggg                                   >  2:1146687/1‑139 (MQ=255)
                                                                                                                           tCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAgggg                                   <  1:1146687/139‑1 (MQ=255)
                                                                                                                             gcaatgtaagcgccCGCAGGTACTCGGCCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCa                                                                              <  2:24456‑M1/80‑1 (MQ=255)
                                                                                                                             gcaatgtaagcgccCGCAGGTACTCGGCCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCa                                                                              >  1:24456‑M1/15‑94 (MQ=255)
                                                                                                                             ccAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGt                      <  1:364458/150‑1 (MQ=255)
                                                                                                                             ccAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGt                      >  2:686063/1‑150 (MQ=255)
                                                                                                                              cAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTGCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTg                     >  1:1163992/1‑150 (MQ=255)
                                                                                                                              cAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTg                     >  2:1124197/1‑150 (MQ=255)
                                                                                                                               aGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTgg                    <  1:1830827/150‑2 (MQ=255)
                                                                                                                               aGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGa                    >  1:276036/1‑150 (MQ=255)
                                                                                                                                gCGTATCTGGGCGCAGGCCGAGGACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAg                   >  1:1441072/1‑150 (MQ=255)
                                                                                                                                     tCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAggg              <  1:1891149/150‑1 (MQ=255)
                                                                                                                                     tCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAggg              <  1:191341/150‑1 (MQ=255)
                                                                                                                                     tCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAggg              >  2:736514/1‑150 (MQ=255)
                                                                                                                                       tGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGAt            <  1:686063/150‑1 (MQ=255)
                                                                                                                                        gggCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAgcgc                                                                          <  1:1631191/87‑1 (MQ=255)
                                                                                                                                        gggCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAgcgc                                                                          >  2:1631191/1‑87 (MQ=255)
                                                                                                                                          gcgcAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGAc         >  2:463963/1‑150 (MQ=255)
                                                                                                                                           CGCAGGCCGAGTACCTGGGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACa        >  2:1267931‑M1/1‑150 (MQ=255)
                                                                                                                                               GGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTGTGAGTGCCGAGAGGGGGACGGCGCGGTGAGCAGGGATGACATGcc    <  2:148837‑M1/150‑1 (MQ=255)
                                                                                                                                               GGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGcc    >  1:1124197‑M1/1‑150 (MQ=255)
                                                                                                                                                GCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCa   <  2:98873‑M1/150‑1 (MQ=255)
                                                                                                                                                 ccGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCAt  >  1:1849358/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                
TGCTGCGTGACACCCCGCTGCATGCGTCAATAGACCGGGCTTTCGCCTTTGCCGAGCAATTCGGGGTGAAGGATGCGGCCGTACTGGCCATGCTGAATGTGGATGGCCGGGTGCTTGATGCCTCCCAGCGTATCTGGGCGCAGGCCGAGTACCTGCGGGCGCTTACATTGCGTGAAGGCAGTGAAGCAAAGCTCGCGGGGCAGTTGGGCGCCCTTTCAGCGCGGTTCCTGCGTGACGCCGGCTGGTATGAGTGCCGAGGGGGTGACGGCGCGGTGAGCAGGGATGACATGCCAT  >  NC_002947/5527317‑5527610

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