Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,186,799 G→T 5.3% Q666H (CAG→CAT purL → phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,186,7990GT 5.3% 297.9 / 4.7 94Q666H (CAG→CATpurLphosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (51/38);  new base T (2/3);  total (53/41)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

GAAGCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCGGCGA  >  NC_002947/1186650‑1186948
                                                                                                                                                     |                                                                                                                                                     
gAAGCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGAAAAGCTCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAg                                                                                                                                                       <  1:1953929/150‑1 (MQ=255)
   gCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATg                                                                                                                                                    <  2:192748/150‑1 (MQ=255)
   gCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATg                                                                                                                                                    <  2:1155809/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)
                                                                      accaccCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAg                                                                                                                                                       <  2:821417/80‑1 (MQ=255)
                                                                      accaccCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAg                                                                                                                                                       >  1:821417/1‑80 (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                                                                    >  1:872837/1‑150 (MQ=255)
                                                                                   ggCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggc                                                                    <  2:1022940/150‑1 (MQ=255)
                                                                                   ggCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggc                                                                    >  1:567083/1‑150 (MQ=255)
                                                                                     ccAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggcca                                                                  <  1:605821/150‑1 (MQ=255)
                                                                                      cAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggccat                                                                 >  2:614057/1‑150 (MQ=255)
                                                                                      cAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggccat                                                                 >  2:2006547/1‑150 (MQ=255)
                                                                                      cAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAgccatggccat                                                                 >  1:999686/1‑150 (MQ=255)
                                                                                            gagCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGt                                                                               >  1:1624568/1‑130 (MQ=255)
                                                                                            gagCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGt                                                                               <  2:1624568/130‑1 (MQ=255)
                                                                                            gagCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGa                                                           <  1:245117/150‑1 (MQ=255)
                                                                                               cTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGTAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCg                                                        >  2:1182669/1‑150 (MQ=255)
                                                                                                ttCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCATCTGGTCgcgg                                                                                                                                             <  2:1946088/64‑4 (MQ=255)
                                                                                                ttCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCATCTGGTCgcgg                                                                                                                                             >  1:1946088/1‑61 (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                                                                                                            >  2:2230309/1‑74 (MQ=255)
                                                                                                                       cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTa                                                                                                            <  1:2230309/74‑1 (MQ=255)
                                                                                                                       cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGc                                >  1:1284869/1‑150 (MQ=255)
                                                                                                                       cACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGc                                >  1:2266428/1‑150 (MQ=255)
                                                                                                                        aCCCTCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTCGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGcc                               >  2:906254/1‑150 (MQ=255)
                                                                                                                        aCCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGcc                               >  1:12286/1‑150 (MQ=255)
                                                                                                                          cATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGa                                                           <  2:1544531/120‑1 (MQ=255)
                                                                                                                          cATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGa                                                           >  1:1544531/1‑120 (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                                               <  1:359140/119‑1 (MQ=255)
                                                                                                                                       gTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCt                                               >  2:359140/1‑119 (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)
                                                                                                                                                     |                                                                                                                                                     
GAAGCCGAGCTGGGCGATGACTTCGACCCGAGCGAGCTGGACCTGGACAGCGCCGTCCAGCGCGTGCTGAACCACCCGGCCGTGGCCAGCAAGAGCTTCCTGATCACCATTGGCGACCGCACCATCACCGGCCTGGTTGCCCGCGACCAGATGGTCGGCCCGTGGCAGGTTCCGGTGGCCGACTGCGCCGTTACCGCCACCAGCTTCGACGTCTACACCGGTGAAGCCATGGCCATGGGCGAGCGTACGCCGCTGGCGCTGCTGGATGCCCCGGCGTCCGGCCGCATGGCCATCGGCGA  >  NC_002947/1186650‑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: