Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,235,258 A→T 6.7% S104C (AGC→TGC)  PP_2834 → D‑galactarate dehydratase/altronate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,235,2580AT 6.7% 175.5 / 4.5 60S104C (AGC→TGC) PP_2834D‑galactarate dehydratase/altronate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (29/26);  new base T (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

AGCCTGTGGCGCAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTCCGCCGTGACAT  >  NC_002947/3235109‑3235403
                                                                                                                                                       |                                                                                                                                                 
aGCCTGTGGCGCAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCa                                                                                                                                                   <  2:286996/150‑1 (MQ=255)
  ccTGTGGCGCAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGc                                                                                                                                                 >  1:1311417/1‑150 (MQ=255)
        gcgcAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCa                                                                                                                                           <  2:107821/150‑1 (MQ=255)
        gcgcAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCa                                                                                                                                           <  2:281875/150‑1 (MQ=255)
          gcAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCTGCCTGGCAgg                                                                                                                                         >  2:1111657/1‑150 (MQ=255)
             aGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAg                                                                                                                                      <  2:1661848/150‑1 (MQ=255)
                 cAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTTCCAGGCAGcgcc                                                                                                                                  >  1:400118/1‑150 (MQ=255)
                     cGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCTGCCTGGCAGGCTgcgggtatc                                                                                                                              <  1:1111657/150‑7 (MQ=255)
                       ttGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGAAGCGCCGCGCAg                                                                                                                            <  2:634629/150‑1 (MQ=255)
                        tGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGt                                                                                                                           <  1:1940066/150‑1 (MQ=255)
                           gcCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCACCGCGCAGTTcc                                                                                                                        >  2:148870/1‑150 (MQ=255)
                            ccGCTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCa                                                                                                                       >  2:427876/1‑150 (MQ=255)
                                       aGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTgc                                                                                                            <  2:35142/150‑1 (MQ=255)
                                            atcGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCTGCCTGGCAGGCcgcgggtatc                                                                                                                              <  1:11686/127‑7 (MQ=255)
                                            atcGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCTGCCTGGCAGGCcgcgggtatc                                                                                                                              >  2:11686/1‑121 (MQ=255)
                                            atcGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGAGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGcc                                                                                                       >  1:1850756/1‑150 (MQ=255)
                                             tcGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCg                                                                                                      >  2:348210/1‑150 (MQ=255)
                                               ggATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACCAC‑‑CAAGGTTGCCTACCCGTGGCACGCGAGTCCGCGTTTGAGGGCGATCCACGAAGCATGACAGCCAGCTCGGCGCAGTGCCAGGGCATCGTGCGCGCCGAt                                                                                                    <  1:443219/150‑1 (MQ=255)
                                                   ttGCGTCCCAGGCGATCGACGCTGGCTACCATGTGCATGTGCACAACGTGCAGATCAGCGA‑‑TCTTGCCCGCGCCTATGCGTTTTGTGTCGATATCCCCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGcc                                                                                                <  2:867493/150‑1 (MQ=255)
                                                      cGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCgtg                                                                                             >  2:312426/1‑150 (MQ=255)
                                                       gTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCgtgt                                                                                            <  2:1764887/150‑1 (MQ=255)
                                                          ccAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGa                                                                                                     <  1:1419001/138‑1 (MQ=255)
                                                          ccAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGa                                                                                                     >  2:1419001/1‑138 (MQ=255)
                                                            aGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCa                                                                                                                             <  2:856995/112‑1 (MQ=255)
                                                            aGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCa                                                                                                                             >  1:856995/1‑112 (MQ=255)
                                                            aGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCa                                                                                       >  1:270848/1‑150 (MQ=255)
                                                            aGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCa                                                                                       <  2:829363/150‑1 (MQ=255)
                                                                 aTCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCAc                                                                                      >  1:304511/1‑146 (MQ=255)
                                                                 aTCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCAc                                                                                      <  2:304511/146‑1 (MQ=255)
                                                                  tCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGGGACTAGGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCCCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCa                                                                                 >  1:904276/1‑150 (MQ=255)
                                                                  tCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCa                                                                                                                                                   >  1:1158447/1‑84 (MQ=255)
                                                                  tCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCa                                                                                                                                                   <  2:1158447/84‑1 (MQ=255)
                                                                       gCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTAc                                                                            <  2:667622/150‑1 (MQ=255)
                                                                            cGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCgg                                                                       <  1:312426/150‑1 (MQ=255)
                                                                             gACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGc                                                                                         <  1:1933203/131‑1 (MQ=255)
                                                                             gACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGc                                                                                         >  2:1933203/1‑131 (MQ=255)
                                                                                 atgtgcatgtgcACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCg                                                                        >  2:1134704/1‑144 (MQ=255)
                                                                                 atgtgcatgtgcACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCg                                                                        <  1:1134704/144‑1 (MQ=255)
                                                                                    tgcatgtgcaCAACGTGCAGATGGGCGG‑‑GCGGGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCa                                                               >  1:50638/1‑150 (MQ=255)
                                                                                     gcatgtgcaCAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCAc                                                              >  2:393582/1‑150 (MQ=255)
                                                                                       atgtgcaCAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCt                                                            >  1:856933/1‑150 (MQ=255)
                                                                                           gcaCAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGt                                                        >  2:624961/1‑150 (MQ=255)
                                                                                           gcaCAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCACCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGt                                                        <  1:148870/150‑1 (MQ=255)
                                                                                               aaCGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAAc                                                    <  1:1388432/150‑1 (MQ=255)
                                                                                                     cAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGa                                                                                                     >  1:1702983/1‑95 (MQ=255)
                                                                                                     cAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGa                                                                                                     <  2:1702983/95‑1 (MQ=255)
                                                                                                     cAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCg                                              <  2:1702961/150‑1 (MQ=255)
                                                                                                         tGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGa                                          >  1:1697659/1‑150 (MQ=255)
                                                                                                                  ttttGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGt                                                        >  2:350579/1‑129 (MQ=255)
                                                                                                                  ttttGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGt                                                        <  1:350579/129‑1 (MQ=255)
                                                                                                                         cgcgACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGcc                            <  2:399558/150‑1 (MQ=255)
                                                                                                                                 tGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGAc                                         <  2:1880257/129‑1 (MQ=255)
                                                                                                                                 tGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGAc                                         >  1:1880257/1‑129 (MQ=255)
                                                                                                                                 tGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGa                    <  2:385489/150‑1 (MQ=255)
                                                                                                                                    gTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCa                 >  1:173229/1‑150 (MQ=255)
                                                                                                                                       tGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCAttt              <  1:624961/150‑1 (MQ=255)
                                                                                                                                         gtgtCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTcc            >  2:1619474/1‑150 (MQ=255)
                                                                                                                                            tCGATACCCGCCGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCTGTGTCGCCAGCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTccgcc         <  1:229725/150‑1 (MQ=255)
                                                                                                                                            tCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTccgcc         >  1:1308154/1‑150 (MQ=255)
                                                                                                                                            tCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTccgcc         >  1:1817768/1‑150 (MQ=255)
                                                                                                                                               aTACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTCCGCCGTg      >  2:1861373/1‑150 (MQ=255)
                                                                                                                                               aTACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCACGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGAGCGTGGCCCGGGCCGTGGCCGACCATTTCCGCCGTg      >  2:1038189/1‑150 (MQ=255)
                                                                                                                                                   ccGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTCCGCCGTGACAt  >  2:348205/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                 
AGCCTGTGGCGCAAGGCCAGCCGTTGCGCCGTTACGGGCAGATCATCGGATTTGCGTCCCAGGCGATCGACGCTGGCGACCATGTGCATGTGCACAACGTGCAGATGGGCGA‑‑TTTTGCCCGCGACTATGCGTTTGGTGTCGATACCCGCAGCCTGCCAGGCAGCGCCGCGCAGTTCCAGGGCATCGTGCGCGCCGATGGCCGTGTCGCCACCCGCAACTACATCGGCATTCTCACCTCGGTGAACTGTTCGGCGACCGTGGCCCGGGCCGTGGCCGACCATTTCCGCCGTGACAT  >  NC_002947/3235109‑3235403

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