Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,259,122 T→G 6.3% E176D (GAA→GAC) ‡ tdcG ← L‑serine deaminase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,259,1220TG 6.3% 409.5 / 15.9 140E176D (GAA→GAC) ‡tdcGL‑serine deaminase III
Reads supporting (aligned to +/- strand):  ref base T (68/63);  new base G (3/6);  total (71/69)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 141 total aligned reads displayed.

GTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGC  >  NC_000913/3258973‑3259269
                                                                                                                                                     |                                                                                                                                                   
gtgtTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGt                                                                                                                                                     <  2:1907055/150‑1 (MQ=255)
     cATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCAcc                                                                                                                                                <  2:691629/150‑1 (MQ=255)
       tCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctg                                                                                                                                              >  1:874754/1‑150 (MQ=255)
         ccACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctgct                                                                                                                                            >  1:778152/1‑150 (MQ=255)
          cACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctgctg                                                                                                                                           <  2:1063683/150‑1 (MQ=255)
            cGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaa                                                                                                                                         <  2:1151216/150‑2 (MQ=255)
            cGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaa                                                                                                                                         <  2:2225884/150‑2 (MQ=255)
            cGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaa                                                                                                                                         <  2:2418073/150‑2 (MQ=255)
             gTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGGTGAACTGCTGaaa                                                                                                                                        <  1:1748700/150‑3 (MQ=255)
              ttCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTg                                                                                                                                       <  2:1977082/150‑1 (MQ=255)
               tCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTgg                                                                                                                                      <  1:1253767/150‑1 (MQ=255)
                cAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGa                                                                                                                                     >  1:1452424/1‑150 (MQ=255)
                cAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGa                                                                                                                                     >  2:683384/1‑150 (MQ=255)
                   tACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACctgctg                                                                                                                                           >  1:847740/1‑141 (MQ=255)
                    aCCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATc                                                                                                                                 <  1:1227489/150‑1 (MQ=255)
                        tCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTaa                                                                                                                             >  1:2352210/1‑150 (MQ=255)
                        tCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTaa                                                                                                                             <  1:405591/150‑1 (MQ=255)
                        gcATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTaa                                                                                                                             <  2:2067299/149‑1 (MQ=255)
                           tGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGt                                                                                                                          >  1:484640/1‑150 (MQ=255)
                            gCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTa                                                                                                                         <  2:543077/150‑1 (MQ=255)
                            gCATCACTTCCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTa                                                                                                                         <  1:1276732/150‑1 (MQ=255)
                                cACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGa                                                                                                                     >  2:2110934/1‑150 (MQ=255)
                                     gCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGttt                                                                                                                >  2:850181/1‑150 (MQ=255)
                                       cAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCg                                                                                                              <  2:897717/150‑1 (MQ=255)
                                          aTACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCTCACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACa                                                                                                           <  2:831283/150‑1 (MQ=255)
                                           tACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACAt                                                                                                          >  1:456641/1‑150 (MQ=255)
                                              gggCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGt                                                                                                       >  1:2267077/1‑150 (MQ=255)
                                               ggCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTg                                                                                                      >  2:1847836/1‑150 (MQ=255)
                                               ggCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTg                                                                                                      >  1:2288652/1‑150 (MQ=255)
                                               ggCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTg                                                                                                      >  1:1308482/1‑150 (MQ=255)
                                               ggCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTg                                                                                                      >  1:1685691/1‑150 (MQ=255)
                                                gCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAAATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGc                                                                                                     <  1:418580/150‑1 (MQ=255)
                                                 cAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCg                                                                                                    >  1:614246/1‑150 (MQ=255)
                                                  aaaaCCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGa                                                                                                   >  1:2459098/1‑150 (MQ=255)
                                                    aaCCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACa                                                                                                 >  1:2368133/1‑150 (MQ=255)
                                                      ccGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAgg                                                                                               >  2:631599/1‑150 (MQ=255)
                                                       cGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGc                                                                                              <  2:2267077/150‑1 (MQ=255)
                                                        ggCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGcc                                                                                             <  2:1485205/150‑1 (MQ=255)
                                                         gCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCg                                                                                            <  1:2200365/150‑1 (MQ=255)
                                                          cGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGa                                                                                           <  1:2290687/150‑1 (MQ=255)
                                                           gTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGaa                                                                                          >  2:1017750/1‑150 (MQ=255)
                                                            tCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAg                                                                                         >  2:1563741/1‑150 (MQ=255)
                                                             cAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAgt                                                                                        <  1:804174/150‑1 (MQ=255)
                                                              aaTATCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGGATCATCAGACCAGATATAGTCAGGCGGTTGTAATGACACATTTTCAGCAGTTCAGCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAgtg                                                                                       <  2:2071798/150‑1 (MQ=255)
                                                                  tCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttct                                                                                   <  2:1197389/150‑1 (MQ=255)
                                                                      cTTTGCTGCGCAGCGCTAGCTCGTTATGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttctt                                                                               >  1:1305170/1‑150 (MQ=255)
                                                                       tttGCTTCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttcttc                                                                              <  1:1694700/150‑1 (MQ=255)
                                                                       tttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGttcttcttc                                                                              <  2:484640/150‑1 (MQ=255)
                                                                        ttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCg                                                                             <  1:610449/150‑1 (MQ=255)
                                                                        ttGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCg                                                                             <  2:110090/150‑1 (MQ=255)
                                                                          gctgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGcc                                                                                             <  1:1029393/132‑1 (MQ=255)
                                                                          gctgcGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGAc                                                                           >  2:1508053/1‑150 (MQ=255)
                                                                                cAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATaaa                                                                     >  2:373685/1‑150 (MQ=255)
                                                                                    gcTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACccg                                                                 >  1:398164/1‑150 (MQ=255)
                                                                                       aGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACccgccg                                                              >  2:967898/1‑150 (MQ=255)
                                                                                         cTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACccgccgcc                                                            >  1:905182/1‑150 (MQ=255)
                                                                                            gTTGTGCATCATCAGACGAGATATAGACTGTCCGTTGTAATCACACATTTTCAGCAGTTCACCTGGTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGGTCTTCTTCGACAATAAACCCGCCGCCGAc                                                         <  2:1898768/150‑1 (MQ=255)
                                                                                            gTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGAc                                                         >  1:90906/1‑150 (MQ=255)
                                                                                             ttGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACa                                                        <  1:683384/150‑1 (MQ=255)
                                                                                              tgtgCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACag                                                       <  2:107224/150‑1 (MQ=255)
                                                                                              tgtgCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACag                                                       <  1:1472588/150‑1 (MQ=255)
                                                                                               gtgCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACaga                                                      >  2:1927761/1‑150 (MQ=255)
                                                                                               gtgCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGATTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACtga                                                      >  1:2449321/1‑150 (MQ=255)
                                                                                                tgCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACagag                                                     >  1:969784/1‑150 (MQ=255)
                                                                                                 gCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGt                                                    <  1:1182904/150‑1 (MQ=255)
                                                                                                   atcatcAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTaa                                                  <  1:1940987/150‑1 (MQ=255)
                                                                                                    tcatcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATACACCCGCGGCCGACAGAGTAAt                                                 >  2:235795/1‑150 (MQ=255)
                                                                                                     catcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGc                                                                                                     <  1:396266/97‑1 (MQ=255)
                                                                                                     catcaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATa                                                >  2:1939027/1‑150 (MQ=255)
                                                                                                        aaGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGt                                             <  1:1847836/149‑1 (MQ=255)
                                                                                                         aGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGtt                                            >  1:2306391/1‑150 (MQ=255)
                                                                                                          gACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGttt                                           >  2:260491/1‑150 (MQ=255)
                                                                                                          gACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGttt                                           >  2:251776/1‑150 (MQ=255)
                                                                                                              aGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGTCGCCGACAGAGTAATAGGTTTTACt                                       <  1:1133367/150‑1 (MQ=255)
                                                                                                                    aGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAg                                 >  2:2249532/1‑150 (MQ=255)
                                                                                                                      aCAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGct                               <  2:543889/150‑1 (MQ=255)
                                                                                                                       cAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGctc                              <  1:2110934/150‑1 (MQ=255)
                                                                                                                         ggCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCtt                            >  2:298012/1‑150 (MQ=255)
                                                                                                                            cGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCt                         >  1:938723/1‑150 (MQ=255)
                                                                                                                                gTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTcc                     <  2:778152/150‑1 (MQ=255)
                                                                                                                                 tAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGAAAGAGTAATAGGTTTTACTTAATAGATCTTCCTGTcac                    >  2:2023411/1‑148 (MQ=255)
                                                                                                                                   aTCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCtt                  >  2:2160681/1‑150 (MQ=255)
                                                                                                                                     cacacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTcc                >  1:574066/1‑150 (MQ=255)
                                                                                                                                     cacacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTcc                >  2:78214/1‑150 (MQ=255)
                                                                                                                                     cacacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTcc                >  1:1038428/1‑150 (MQ=255)
                                                                                                                                     cacacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTcc                >  1:2232867/1‑150 (MQ=255)
                                                                                                                                     cacacaTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTcc                <  1:2132428/150‑1 (MQ=255)
                                                                                                                                       tttccactCAGCAGGTGAACTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAg              >  2:1700495/8‑150 (MQ=255)
                                                                                                                                         caTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGc            <  1:1927761/150‑1 (MQ=255)
                                                                                                                                         caTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGc            >  1:1437230/1‑150 (MQ=255)
                                                                                                                                          aTTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCa           >  2:935154/1‑150 (MQ=255)
                                                                                                                                           ttttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAg          <  2:1685691/150‑1 (MQ=255)
                                                                                                                                           cactCAGCAGGTGAACTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAg          >  2:2138279/4‑150 (MQ=255)
                                                                                                                                            tttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGt         <  1:1017750/150‑1 (MQ=255)
                                                                                                                                            tttCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCATCCAGGCAGt         <  2:367121/150‑1 (MQ=255)
                                                                                                                                            actCAGCAGGTGAACTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGt         >  1:2183106/3‑150 (MQ=255)
                                                                                                                                               aagcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGAt      <  1:2459903/149‑1 (MQ=255)
                                                                                                                                                agcagTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATc     >  1:1944143/1‑150 (MQ=255)
                                                                                                                                                   agTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGc  >  2:534606/1‑150 (MQ=255)
                                                                                                                                                   agTTCACATGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGc  >  1:2201984/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTGTTCATCCCACGTTCAATACCGTCATGCATCACTTGCCAGATACGGGCAAAACCGGCGTCAATTTCCGCTTTGCTGCGCAGCGCTAGCTCGTTGTGCATCATCAGACCAGATATAGACAGGCCGTTGTAATCACACATTTTCAGCAGTTCACCTGCTGAGTGGAAATCGTAAGGTACGGACGTTTCGACATCGTGCGACAGGCCGAAGTGTTCTTCTTCGACAATAAACCCGCCGCCGACAGAGTAATAGGTTTTACTTAATAGCTCTTCCTGTCCCTTCCAGGCAGTGATCCGC  >  NC_000913/3258973‑3259269

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