Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,186,800 A→C 5.3% M667L (ATG→CTG)  purL → phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,186,8000AC 5.3% 279.3 / 5.8 95M667L (ATG→CTG) purLphosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base A (52/38);  new base C (2/3);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCGGCGA  >  NC_002947/1186653‑1186948
                                                                                                                                                   |                                                                                                                                                    
gCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATg                                                                                                                                                    <  2:1155809/150‑1 (MQ=255)
gCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATg                                                                                                                                                    <  2:192748/150‑1 (MQ=255)
   gAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTc                                                                                                                                                 >  2:1927439/1‑150 (MQ=255)
    aGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCg                                                                                                                                                <  1:553515/150‑1 (MQ=255)
       tGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCATCTGGTCgcgg                                                                                                                                             <  1:342618/150‑4 (MQ=255)
       tGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCATCTGGTCgcgg                                                                                                                                             <  1:1319607/150‑4 (MQ=255)
       tGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGcc                                                                                                                                             <  1:756429/150‑1 (MQ=255)
          gCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAATTGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGt                                                                                                                                          <  2:1666088/150‑1 (MQ=255)
               gACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAg                                                                                                                                     <  2:1224581/150‑1 (MQ=255)
                  ttCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGtt                                                                                                                                  >  1:98260/1‑150 (MQ=255)
                  ttCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGtt                                                                                                                                  >  1:504718/1‑150 (MQ=255)
                  ttCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGtt                                                                                                                                  >  1:687947/1‑150 (MQ=255)
                  ttCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGtt                                                                                                                                  >  1:72353/1‑150 (MQ=255)
                      aCCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCgg                                                                                                                              <  2:1431525/150‑1 (MQ=255)
                           agcgagcTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGTCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGcc                                                                                                                         <  1:1364558/150‑1 (MQ=255)
                                    gacctggacAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGcc                                                                                                                <  1:2397656/150‑1 (MQ=255)
                                       ctggacAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGtt                                                                                                             >  1:1388819/1‑150 (MQ=255)
                                        tggacAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTa                                                                                                            >  2:245117/1‑150 (MQ=255)
                                         ggacAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTAc                                                                                                           <  2:504718/150‑1 (MQ=255)
                                                   gTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATAACTATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGCCCGACTGCGCCGTTACCGCCACCAGc                                                                                                 <  1:1132174/150‑1 (MQ=255)
                                                    tCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCt                                                                                                <  1:1786128/150‑1 (MQ=255)
                                                        gcgcGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGAcca                                                                                                                    <  1:1328390/126‑3 (MQ=255)
                                                        gcgcGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGAcca                                                                                                                    >  2:1328390/1‑124 (MQ=255)
                                                        gcgcGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGccacca                                                                                                   >  1:1196726/1‑143 (MQ=255)
                                                        gcgcGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGccacca                                                                                                   <  2:1196726/143‑1 (MQ=255)
                                                          gcgTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTAAGTACCCACTGTTCACGACCGCACCATCACCGGCGTTGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACg                                                                                          >  1:1022940/1‑150 (MQ=255)
                                                             tGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCt                                                                                       >  2:286399/1‑150 (MQ=255)
                                                                tGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTaca                                                                                    >  1:1540207/1‑150 (MQ=255)
                                                                    ccaccCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGTTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCgg                                                                                <  1:903967/150‑1 (MQ=255)
                                                                      accCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCTCCCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGGCTACCCCGGTg                                                                              >  2:162619/1‑150 (MQ=255)
                                                                      accCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTg                                                                              <  1:299236/150‑1 (MQ=255)
                                                                        ccGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGTTCTGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTCCCGCCACCAGCTTCGACGTCTACACCGGTGaa                                                                            >  1:1621058/1‑150 (MQ=255)
                                                                        ccGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGaa                                                                            >  1:325629/1‑150 (MQ=255)
                                                                             cGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccat                                                                       <  1:621676/150‑1 (MQ=255)
                                                                                ggCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggc                                                                    <  2:1022940/150‑1 (MQ=255)
                                                                                ggCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggc                                                                    >  1:872837/1‑150 (MQ=255)
                                                                                ggCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggc                                                                    >  1:567083/1‑150 (MQ=255)
                                                                                  ccAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggcca                                                                  <  1:605821/150‑1 (MQ=255)
                                                                                   cAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggccat                                                                 >  1:999686/1‑150 (MQ=255)
                                                                                   cAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggccat                                                                 >  2:614057/1‑150 (MQ=255)
                                                                                   cAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggccat                                                                 >  2:2006547/1‑150 (MQ=255)
                                                                                         gagCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGt                                                                               <  2:1624568/130‑1 (MQ=255)
                                                                                         gagCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGt                                                                               >  1:1624568/1‑130 (MQ=255)
                                                                                         gagCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGa                                                           <  1:245117/150‑1 (MQ=255)
                                                                                            cTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGTAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCg                                                        >  2:1182669/1‑150 (MQ=255)
                                                                                             ttCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCATCTGGTCgcgg                                                                                                                                             >  1:1946088/1‑61 (MQ=255)
                                                                                             ttCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCATCTGGTCgcgg                                                                                                                                             <  2:1946088/64‑4 (MQ=255)
                                                                                              tCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGCTCCTGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTa                                                      <  1:162619/150‑1 (MQ=255)
                                                                                                 tGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTACGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTAcgc                                                   >  1:1094348/1‑150 (MQ=255)
                                                                                                     cACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCt                                               >  2:1142510/1‑150 (MQ=255)
                                                                                                      aCCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGCCCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTg                                              >  1:2294655/1‑150 (MQ=255)
                                                                                                        cATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGgc                                            >  2:195520/1‑150 (MQ=255)
                                                                                                        cATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGgc                                            >  1:1838149/1‑150 (MQ=255)
                                                                                                         aTTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGgcg                                           >  2:771868/1‑150 (MQ=255)
                                                                                                          ttGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGgcgc                                          <  1:1927439/150‑1 (MQ=255)
                                                                                                           tGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCgct                                         <  2:2172111/150‑1 (MQ=255)
                                                                                                             gCGACCGCACCATCACCGGCCTGGTTGCCCGCGGGCCGATGGTCCGTCCGTGGCAGGTTCCGGTCGCTGAATGTGCCAGTACCGCCGCGAGCTTCGATGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCgctgc                                       <  2:2077452/150‑1 (MQ=255)
                                                                                                             gCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGccatggccatgg                                                               <  1:831613/126‑1 (MQ=255)
                                                                                                             gCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGccatggccatgg                                                               >  2:831613/1‑126 (MQ=255)
                                                                                                              cGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCgctgct                                      <  1:1221968/150‑1 (MQ=255)
                                                                                                              cGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCgctgct                                      >  2:1743574/1‑150 (MQ=255)
                                                                                                                  cGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGAt                                  >  2:1971556/1‑150 (MQ=255)
                                                                                                                    cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTa                                                                                                            <  1:2230309/74‑1 (MQ=255)
                                                                                                                    cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTa                                                                                                            >  2:2230309/1‑74 (MQ=255)
                                                                                                                    cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGc                                >  1:2266428/1‑150 (MQ=255)
                                                                                                                    cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGc                                >  1:1284869/1‑150 (MQ=255)
                                                                                                                     aCCCTCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTCGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGcc                               >  2:906254/1‑150 (MQ=255)
                                                                                                                     aCCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGcc                               >  1:12286/1‑150 (MQ=255)
                                                                                                                       cATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGa                                                           >  1:1544531/1‑120 (MQ=255)
                                                                                                                       cATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGa                                                           <  2:1544531/120‑1 (MQ=255)
                                                                                                                           aCCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCg                         >  1:1438345/1‑150 (MQ=255)
                                                                                                                           aCCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCg                         >  2:2148024/1‑150 (MQ=255)
                                                                                                                            ccGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGt                        >  1:2099083/1‑150 (MQ=255)
                                                                                                                               gCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCg                     <  1:286399/150‑1 (MQ=255)
                                                                                                                                  tGGTCGCCCGCGGCCAGAGGGTCGGCTCGTGGCAGGTTCCGGTGGCCGACTGCCCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGCATGCCCCGGCGTCCGGcc                  <  2:72353/150‑1 (MQ=255)
                                                                                                                                   ggTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCg                 >  2:1928215/1‑150 (MQ=255)
                                                                                                                                    gTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCt                                               >  2:359140/1‑119 (MQ=255)
                                                                                                                                    gTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCt                                               <  1:359140/119‑1 (MQ=255)
                                                                                                                                    gTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGc                >  1:2409748/1‑150 (MQ=255)
                                                                                                                                      tGCCCGCGACCAGATGGCCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCAt              <  2:98260/150‑1 (MQ=255)
                                                                                                                                        cccGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATgg            <  2:1838149/150‑1 (MQ=255)
                                                                                                                                         ccGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGc           >  2:1732255/1‑150 (MQ=255)
                                                                                                                                         ccGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGc           >  2:1521326/1‑150 (MQ=255)
                                                                                                                                         ccGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGc           >  1:95469/1‑150 (MQ=255)
                                                                                                                                         ccGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGc           >  1:945874/1‑150 (MQ=255)
                                                                                                                                         ccGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGc           >  1:2360511/1‑150 (MQ=255)
                                                                                                                                          gccgACCATCTGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGcc          >  2:101469/3‑150 (MQ=255)
                                                                                                                                          cgcgACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGc                   <  2:382350/141‑1 (MQ=255)
                                                                                                                                          cgcgACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGc                   >  1:382350/1‑141 (MQ=255)
                                                                                                                                             gACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCACCGGCGTTCGGCCGCATGGCCATc       <  2:1419075/150‑1 (MQ=255)
                                                                                                                                              aCCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCg      <  2:2294655/150‑1 (MQ=255)
                                                                                                                                                cAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGcca                                                                                                      <  2:258090/52‑1 (MQ=255)
                                                                                                                                                cAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGcca                                                                                                      >  1:258090/1‑52 (MQ=255)
                                                                                                                                                cAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCGGc    >  2:218287/1‑150 (MQ=255)
                                                                                                                                                  gATGGTCGGCCCGTGGCAGGTTCCGGTGGCAGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCGGCGa  >  1:936232/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCGGCGA  >  NC_002947/1186653‑1186948

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