Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,808 G→T 7.9% 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,8080GT 7.9% 1188.3 / 50.8 314Q97K (CAG→AAG) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (132/157);  new base T (12/13);  total (144/170)
Fisher's exact test for biased strand distribution p-value = 8.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 325 total aligned reads displayed.

ATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTG  >  CP009273/2589659‑2589957
                                                                                                                                                     |                                                                                                                                                     
  aGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgtt                                                                                                                                                     <  2:3472304/150‑1 (MQ=255)
  aGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgtt                                                                                                                                                     >  2:3761372/1‑150 (MQ=255)
   gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttg                                                                                                                                                    >  2:4907849/1‑150 (MQ=255)
   gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttg                                                                                                                                                    <  2:2868681/150‑1 (MQ=255)
   gAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACtgttg                                                                                                                                                    <  2:3848341/150‑1 (MQ=255)
     aGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTc                                                                                                                                                  >  1:455629/1‑150 (MQ=255)
      ggAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcc                                                                                                                                                 <  2:4827315/150‑1 (MQ=255)
       gagagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTccc                                                                                                                                                >  1:211305/1‑150 (MQ=255)
        agagaTCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTcccc                                                                                                                                               <  1:3322761/150‑1 (MQ=255)
              cGATATAAACCGTGCCCTCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgtgt                                                                                                                                         >  1:1913813/1‑150 (MQ=255)
              cGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGgtgt                                                                                                                                         >  2:2669501/1‑150 (MQ=255)
               gATATCAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCAGTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                        <  2:3818643/150‑1 (MQ=255)
               gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                        >  1:2072020/1‑150 (MQ=255)
               gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                        >  1:3448576/1‑150 (MQ=255)
               gATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTc                                                                                                                                        >  2:1386417/1‑150 (MQ=255)
                atatAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCt                                                                                                                                       >  2:3322732/1‑150 (MQ=255)
                  ataAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTc                                                                                                                                     <  1:4945830/150‑1 (MQ=255)
                   taAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTcc                                                                                                                                    <  1:1060111/150‑1 (MQ=255)
                     aaCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGt                                                                                                                                  >  1:4705591/1‑150 (MQ=255)
                      aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                 >  2:3216873/1‑150 (MQ=255)
                      aCCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTg                                                                                                                                 >  1:4204070/1‑150 (MQ=255)
                       ccGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTgg                                                                                                                                >  2:417033/1‑150 (MQ=255)
                         gTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGtt                                                                                                                              >  1:1444640/1‑150 (MQ=255)
                         gTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGtt                                                                                                                              <  2:1564040/150‑1 (MQ=255)
                              aTCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAg                                                                                                                         >  2:2248098/1‑150 (MQ=255)
                                 cgccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAAt                                                                                                                      <  1:1712934/150‑1 (MQ=255)
                                   ccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCa                                                                                                                    <  1:1748806/150‑1 (MQ=255)
                                   ccgGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCa                                                                                                                    >  1:714033/1‑150 (MQ=255)
                                      ggCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCg                                                                                                                 <  2:60336/150‑1 (MQ=255)
                                        cactAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCTGAATCACCGAg                                                                                                               <  2:1814994/147‑1 (MQ=255)
                                              aaCGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAAt                                                                                                         <  2:5066028/150‑1 (MQ=255)
                                              aaCGGCCACATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTAACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTAGCCAGAATCACCGAGGTGAAt                                                                                                         >  1:437484/1‑150 (MQ=255)
                                                   ccccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGc                                                                                                    >  2:1141856/1‑150 (MQ=255)
                                                   ccccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGc                                                                                                    <  2:2205094/150‑1 (MQ=255)
                                                    cccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCg                                                                                                   <  2:1730110/150‑1 (MQ=255)
                                                     ccATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCgg                                                                                                  <  1:3949835/150‑1 (MQ=255)
                                                        ttATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCtt                                                                                               >  2:1618334/1‑150 (MQ=255)
                                                        ttATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCtt                                                                                               >  1:1151581/1‑150 (MQ=255)
                                                         tATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTc                                                                                              <  1:934738/150‑1 (MQ=255)
                                                         tATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTc                                                                                              >  2:3133284/1‑150 (MQ=255)
                                                         tATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCAGAATCACCGAGGTGAATTTTGCGGCTTCg                                                                                             >  2:5101602/1‑150 (MQ=255)
                                                           tGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGt                                                                                            >  1:1217123/1‑150 (MQ=255)
                                                            ggACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTc                                                                                           >  2:529868/1‑150 (MQ=255)
                                                               cTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTc                                                                                        <  2:4100495/150‑1 (MQ=255)
                                                               cTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTc                                                                                        <  1:1933390/150‑1 (MQ=255)
                                                                tGGCCCGCCCCGCAGCCGGTACGCGTCATTCCTCGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCa                                                                                       <  2:3651462/150‑1 (MQ=255)
                                                                 ggCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAt                                                                                      <  2:2624155/150‑1 (MQ=255)
                                                                  gcccgcccCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCAtt                                                                                     <  2:2585472/150‑1 (MQ=255)
                                                                    ccgccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTg                                                                                   <  1:4817570/150‑1 (MQ=255)
                                                                      gccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTggg                                                                                 >  1:864225/1‑150 (MQ=255)
                                                                      gccccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTggg                                                                                 <  1:2629618/150‑1 (MQ=255)
                                                                         ccgcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTa                                                                              >  1:4749932/1‑150 (MQ=255)
                                                                           gcAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCGGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAAt                                                                            <  1:4237043/150‑1 (MQ=255)
                                                                             aGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTg                                                                          <  1:2315489/150‑1 (MQ=255)
                                                                             aGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTg                                                                          <  2:2348335/150‑1 (MQ=255)
                                                                             aGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTg                                                                          >  2:2491309/1‑150 (MQ=255)
                                                                              gCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGa                                                                         <  1:4217824/150‑1 (MQ=255)
                                                                              gCCGGTACGCGTCAGTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGa                                                                         <  1:3474163/150‑1 (MQ=255)
                                                                               ccGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAc                                                                        <  1:2543072/150‑1 (MQ=255)
                                                                                cGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAcc                                                                       <  1:4667368/150‑1 (MQ=255)
                                                                                cGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAcc                                                                       <  1:393033/150‑1 (MQ=255)
                                                                                  gTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                                     >  1:1369755/1‑150 (MQ=255)
                                                                                      gcgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTg                                                                 <  1:4447885/150‑1 (MQ=255)
                                                                                       cgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                                <  1:823435/150‑1 (MQ=255)
                                                                                       cgTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                                <  1:3425740/150‑1 (MQ=255)
                                                                                         tCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAtt                                                              <  1:2518909/150‑1 (MQ=255)
                                                                                           aTTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgt                                                            <  1:610929/150‑1 (MQ=255)
                                                                                             tCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttg                                                          <  2:1276871/150‑1 (MQ=255)
                                                                                                  ggTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTggg                                                     >  2:2751900/1‑150 (MQ=255)
                                                                                                     aCATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGaaa                                                  <  2:342383/150‑1 (MQ=255)
                                                                                                      cATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAc                                                 <  1:1084134/150‑1 (MQ=255)
                                                                                                       aTGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAcc                                                <  1:3074546/150‑1 (MQ=255)
                                                                                                       aTGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAcc                                                <  2:1188310/150‑1 (MQ=255)
                                                                                                       aTGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAAcc                                                <  1:4134340/150‑1 (MQ=255)
                                                                                                         gACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCgg                                              <  2:2469947/150‑1 (MQ=255)
                                                                                                           ccAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCt                                            <  2:4055883/150‑1 (MQ=255)
                                                                                                               cTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                                        >  2:4391110/1‑150 (MQ=255)
                                                                                                                ttCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCg                                       >  2:5087682/1‑150 (MQ=255)
                                                                                                                 tCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCgg                                      >  2:812259/1‑150 (MQ=255)
                                                                                                                   ggTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                    <  1:4907849/150‑1 (MQ=255)
                                                                                                                   ggTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtc                                    <  1:2171071/150‑1 (MQ=255)
                                                                                                                     ttGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                                  >  1:1870735/1‑150 (MQ=255)
                                                                                                                        cTGATAGGTCTTACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCaa                               <  1:2248022/150‑1 (MQ=255)
                                                                                                                         tGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAc                              >  2:1488760/1‑150 (MQ=255)
                                                                                                                         tGATAGGTCTTACCCATCTACTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGAGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCTGTCATCAAc                              >  2:3082705/1‑150 (MQ=255)
                                                                                                                          gATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAAcc                             >  2:250784/1‑150 (MQ=255)
                                                                                                                             aGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCga                          <  2:4705591/150‑2 (MQ=255)
                                                                                                                              ggTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCgat                         >  1:5101714/1‑148 (MQ=255)
                                                                                                                               gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                        <  1:3652777/150‑1 (MQ=255)
                                                                                                                                tCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTa                       <  2:2059737/150‑1 (MQ=255)
                                                                                                                                 cTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTaa                      <  2:4782741/150‑1 (MQ=255)
                                                                                                                                    aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATc                   <  2:4610650/150‑1 (MQ=255)
                                                                                                                                    aCCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATc                   <  1:527922/150‑1 (MQ=255)
                                                                                                                                      ccATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCaa                 <  2:69349/150‑1 (MQ=255)
                                                                                                                                       cATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac                <  2:3351234/150‑1 (MQ=255)
                                                                                                                                        aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAaca               >  1:1670322/1‑150 (MQ=255)
                                                                                                                                          cttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACAcc             <  1:2943153/150‑1 (MQ=255)
                                                                                                                                            tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAt           <  1:3254287/150‑1 (MQ=255)
                                                                                                                                              tcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag         <  1:2440500/150‑1 (MQ=255)
                                                                                                                                                 aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGTTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag      >  1:5135693/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATAGAAGGAGAGATCGATATAAACCGTGCCATCCGCCGGGCAATAGAACGGCCCCATTATGGACTGGCCCGCCCCGCAGCCGGTACGCGTCATTCCACGGTACATGACCAGCTTCGGTTGCTGATAGGTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTG  >  CP009273/2589659‑2589957

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: