Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,808 G→T 11.7% Q97K (CAG→AAG)  ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,589,8080GT11.7% 1061.4 / 72.2 309Q97K (CAG→AAG) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (129/142);  new base T (17/20);  total (147/162)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 329 total aligned reads displayed.

TAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCT  >  CP009273/2589660‑2589956
                                                                                                                                                    |                                                                                                                                                    
tAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtg                                                                                                                                                      >  2:1605150/1‑149 (MQ=255)
  gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgtt                                                                                                                                                    >  2:305145/1‑149 (MQ=255)
   aaGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttg                                                                                                                                                   <  2:3584238/149‑1 (MQ=255)
    aGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttgt                                                                                                                                                  >  1:2756136/1‑149 (MQ=255)
     ggAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTc                                                                                                                                                 <  2:799512/149‑1 (MQ=255)
       agagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTccc                                                                                                                                               <  2:919119/149‑1 (MQ=255)
       agagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTccc                                                                                                                                               <  2:2963373/149‑1 (MQ=255)
        gagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcccc                                                                                                                                              <  2:1124067/149‑1 (MQ=255)
        gagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcccc                                                                                                                                              >  1:4356791/1‑149 (MQ=255)
            tCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgt                                                                                                                                          <  2:3479712/149‑1 (MQ=255)
             cGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgtg                                                                                                                                         >  2:2575191/1‑149 (MQ=255)
               atatAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                       >  2:1073436/1‑149 (MQ=255)
               atatAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                       >  2:50535/1‑149 (MQ=255)
               atatAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                       <  1:3281319/149‑1 (MQ=255)
                 ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCtt                                                                                                                                     >  2:4714732/1‑149 (MQ=255)
                 ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCtt                                                                                                                                     <  1:4007626/149‑1 (MQ=255)
                    aaCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCg                                                                                                                                  >  1:2561156/1‑149 (MQ=255)
                     aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGt                                                                                                                                 <  2:586861/149‑1 (MQ=255)
                       cGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGACCTGTTGTCCCCAGGTGTCTTCCGTgg                                                                                                                               <  2:2469218/149‑1 (MQ=255)
                        gTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGt                                                                                                                              <  2:1797158/149‑1 (MQ=255)
                         tGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGtt                                                                                                                             >  2:1287338/1‑149 (MQ=255)
                          gCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTg                                                                                                                            >  2:4014170/1‑149 (MQ=255)
                            cATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGcc                                                                                                                          <  2:3693807/149‑1 (MQ=255)
                               ccgccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGa                                                                                                                       >  2:67030/1‑149 (MQ=255)
                               ccgccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGa                                                                                                                       >  2:2440196/1‑149 (MQ=255)
                                 gccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAAt                                                                                                                     >  1:4757207/1‑149 (MQ=255)
                                     ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCAcc                                                                                                                 >  2:1496621/1‑149 (MQ=255)
                                          tAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGt                                                                                                            >  1:2843605/1‑149 (MQ=255)
                                          tAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGt                                                                                                            >  1:394632/1‑149 (MQ=255)
                                          tAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGt                                                                                                            >  2:4312632/1‑149 (MQ=255)
                                           aGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTg                                                                                                           <  2:1450835/149‑1 (MQ=255)
                                              aCGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAt                                                                                                        >  1:421659/1‑149 (MQ=255)
                                                ggCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAttt                                                                                                      <  2:3104439/149‑1 (MQ=255)
                                                 gCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAtttt                                                                                                     <  2:2985251/149‑1 (MQ=255)
                                                  ccccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTg                                                                                                    >  1:1363935/1‑149 (MQ=255)
                                                   cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCAGTCCACGGTACATGACCAGCTTCGGTTGCTGATCGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGc                                                                                                   <  2:3237799/149‑1 (MQ=255)
                                                      aTTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGc                                                                                                <  1:304734/149‑1 (MQ=255)
                                                         aTGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTc                                                                                             >  1:2382485/1‑149 (MQ=255)
                                                              cTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCtt                                                                                        <  2:4398343/149‑1 (MQ=255)
                                                              cTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCtt                                                                                        <  1:2511597/149‑1 (MQ=255)
                                                              cTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCtt                                                                                        <  1:4779345/149‑1 (MQ=255)
                                                               tGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTc                                                                                       <  1:2259560/149‑1 (MQ=255)
                                                                 gcccgcccCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAt                                                                                     <  1:2424644/149‑1 (MQ=255)
                                                                 gcccgcccCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAt                                                                                     >  2:2129472/1‑149 (MQ=255)
                                                                  cccgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAtt                                                                                    >  1:154230/1‑149 (MQ=255)
                                                                     gccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTgg                                                                                 <  2:3532020/149‑1 (MQ=255)
                                                                               cGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAc                                                                       <  1:687552/149‑1 (MQ=255)
                                                                                 gTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcg                                                                     >  1:2571592/1‑149 (MQ=255)
                                                                                 gTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcg                                                                     >  2:4297890/1‑149 (MQ=255)
                                                                                   aCGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgcg                                                                   >  1:1010123/1‑149 (MQ=255)
                                                                                   aCGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgcg                                                                   >  1:4122490/1‑149 (MQ=255)
                                                                                     gcgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGtt                                                                 <  1:2603516/149‑1 (MQ=255)
                                                                                     gcgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGtt                                                                 >  1:3355257/1‑149 (MQ=255)
                                                                                     gcgTCACTCCACTGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCATGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTTAATTTTGCGGCTTCGTATTCATTTGGGCTAATTGACCGCGtt                                                                 <  1:3613807/149‑1 (MQ=255)
                                                                                      cgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTg                                                                <  1:2059585/149‑1 (MQ=255)
                                                                                       gTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                               >  1:47237/1‑149 (MQ=255)
                                                                                         cATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAtt                                                             <  2:719252/149‑1 (MQ=255)
                                                                                          aTTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttg                                                            <  2:2091480/149‑1 (MQ=255)
                                                                                           ttCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgt                                                           >  1:3763281/1‑149 (MQ=255)
                                                                                           ttCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgt                                                           <  1:1061992/149‑1 (MQ=255)
                                                                                           ttCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgt                                                           <  1:955345/149‑1 (MQ=255)
                                                                                            tCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgtt                                                          <  1:1113237/149‑1 (MQ=255)
                                                                                             ccACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttg                                                         <  2:2210001/149‑1 (MQ=255)
                                                                                               aCGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgtt                                                       >  1:3676598/1‑149 (MQ=255)
                                                                                                    aCATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGaa                                                  >  1:2930288/1‑149 (MQ=255)
                                                                                                     cATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGaaa                                                 <  1:2089266/149‑1 (MQ=255)
                                                                                                       tGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAcc                                               >  1:174322/1‑149 (MQ=255)
                                                                                                        gACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCg                                              <  2:4224703/149‑1 (MQ=255)
                                                                                                         aCCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCgg                                             <  2:465236/149‑1 (MQ=255)
                                                                                                           cAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCt                                           <  1:1880261/149‑1 (MQ=255)
                                                                                                            aGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTg                                          >  1:3437834/1‑149 (MQ=255)
                                                                                                            aGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTg                                          <  2:460556/149‑1 (MQ=255)
                                                                                                             gCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGc                                         >  2:4000489/1‑149 (MQ=255)
                                                                                                               ttCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                       >  2:3173249/1‑149 (MQ=255)
                                                                                                               ttCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                       <  2:1331375/149‑1 (MQ=255)
                                                                                                                 cGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCgg                                     <  2:1249892/149‑1 (MQ=255)
                                                                                                                    ttGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGACACCGGCTGCCCGGtca                                  >  1:613955/1‑149 (MQ=255)
                                                                                                                    ttGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtca                                  >  1:1480901/1‑149 (MQ=255)
                                                                                                                    ttGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtca                                  <  1:2013753/149‑1 (MQ=255)
                                                                                                                    ttGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtca                                  <  2:3847264/149‑1 (MQ=255)
                                                                                                                      gCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcatc                                <  1:3084831/149‑1 (MQ=255)
                                                                                                                         gATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAc                             >  2:4679312/1‑149 (MQ=255)
                                                                                                                           tAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAccc                           >  2:776870/1‑149 (MQ=255)
                                                                                                                            aGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCg                          >  1:2793303/1‑149 (MQ=255)
                                                                                                                                cTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTa                      >  1:1162295/1‑149 (MQ=255)
                                                                                                                                   aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   <  1:4322330/149‑1 (MQ=255)
                                                                                                                                   aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                   <  1:4490432/149‑1 (MQ=255)
                                                                                                                                      cATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCaa                >  1:1705476/1‑149 (MQ=255)
                                                                                                                                          ttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACAcc            >  2:1979796/1‑149 (MQ=255)
                                                                                                                                           tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCa           <  2:4356791/149‑1 (MQ=255)
                                                                                                                                           tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCa           <  2:3367305/149‑1 (MQ=255)
                                                                                                                                           tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCa           <  2:2723048/149‑1 (MQ=255)
                                                                                                                                              cGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag        >  1:277451/1‑149 (MQ=255)
                                                                                                                                               gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt       <  2:2628615/149‑1 (MQ=255)
                                                                                                                                               gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt       <  2:1010123/149‑1 (MQ=255)
                                                                                                                                                aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta      >  1:730888/1‑149 (MQ=255)
                                                                                                                                                 aCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag     <  2:1182149/149‑1 (MQ=255)
                                                                                                                                                  cTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGc    >  1:3139574/1‑149 (MQ=255)
                                                                                                                                                   tgttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGcc   <  2:837638/149‑1 (MQ=255)
                                                                                                                                                    gttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCt  <  1:1605150/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCT  >  CP009273/2589660‑2589956

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