Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 987,229 C→A 10.5% G384G (GGC→GGA pepN → aminopeptidase N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273987,2290CA10.5% 1164.5 / 67.1 328G384G (GGC→GGApepNaminopeptidase N
Reads supporting (aligned to +/- strand):  ref base C (136/157);  new base A (5/30);  total (141/187)
Fisher's exact test for biased strand distribution p-value = 2.21e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 336 total aligned reads displayed.

TTCCCGCGCAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTC  >  CP009273/987082‑987378
                                                                                                                                                   |                                                                                                                                                     
  cccGCGCAGTTAACCGCATCAATCATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCgg                                                                                                                                                   <  2:3187922/150‑1 (MQ=255)
   ccGCGCAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGa                                                                                                                                                  <  1:808257/150‑1 (MQ=255)
   ccGCGCAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGa                                                                                                                                                  <  2:676489/150‑1 (MQ=255)
   ccGCGCAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGa                                                                                                                                                  >  2:3383096/1‑150 (MQ=255)
    cgcgCAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGaa                                                                                                                                                 >  2:4429368/1‑150 (MQ=255)
      cgcAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGt                                                                                                                                               >  2:942001/1‑150 (MQ=255)
        cAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGa                                                                                                                                             >  2:217859/1‑150 (MQ=255)
        cAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGa                                                                                                                                             >  2:2515057/1‑150 (MQ=255)
          gTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGAtt                                                                                                                                           <  2:3779264/150‑1 (MQ=255)
           ttAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTc                                                                                                                                          <  2:1931890/150‑1 (MQ=255)
               ccGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCAt                                                                                                                                      <  2:1693843/150‑1 (MQ=255)
                cGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATg                                                                                                                                     >  2:230731/1‑150 (MQ=255)
                   aTCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATc                                                                                                                                  <  1:392822/150‑1 (MQ=255)
                    tCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATcc                                                                                                                                 >  1:1494897/1‑150 (MQ=255)
                     cAATAATGTACGCACCCTGCGCGGATTGAAGTTTGCAGAAGACGCCTGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCTTTGAGATGAACTACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGTTCca                                                                                                                                <  1:2966742/150‑1 (MQ=255)
                     cAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCca                                                                                                                                <  2:388095/150‑1 (MQ=255)
                     cAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCca                                                                                                                                <  2:2054820/150‑1 (MQ=255)
                        taatGTACGCACCCTGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCacac                                                                                                                             <  1:3267890/150‑1 (MQ=255)
                          atGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACAccc                                                                                                                           <  1:1713369/150‑1 (MQ=255)
                          atGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACAccc                                                                                                                           >  1:3652130/1‑150 (MQ=255)
                               cGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCtt                                                                                                                      >  2:923230/1‑150 (MQ=255)
                               cGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCtt                                                                                                                      <  1:2204444/150‑1 (MQ=255)
                               cGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCtt                                                                                                                      <  1:1986995/150‑1 (MQ=255)
                                 cACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTgg                                                                                                                    <  1:1531118/150‑1 (MQ=255)
                                  aCCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGc                                                                                                                   >  2:4165720/1‑150 (MQ=255)
                                   ccATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCg                                                                                                                  >  1:381286/1‑150 (MQ=255)
                                   ccATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCg                                                                                                                  >  1:161665/1‑150 (MQ=255)
                                    cATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCga                                                                                                                 >  2:781772/1‑150 (MQ=255)
                                     aTGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCgaa                                                                                                                <  2:4060718/150‑1 (MQ=255)
                                      tGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCgaag                                                                                                               <  1:3204381/150‑1 (MQ=255)
                                      tGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCgaag                                                                                                               <  2:299434/150‑1 (MQ=255)
                                         gcgGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGaaa                                                                                                            <  1:2263144/150‑1 (MQ=255)
                                          cgGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGaaaa                                                                                                           <  1:73021/150‑1 (MQ=255)
                                            gATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACt                                                                                                         >  2:858929/1‑150 (MQ=255)
                                            gATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACt                                                                                                         <  1:3106922/150‑1 (MQ=255)
                                            gATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACt                                                                                                         <  1:3106431/150‑1 (MQ=255)
                                            gATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACt                                                                                                         <  1:888625/150‑1 (MQ=255)
                                                  aGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGa                                                                                                   >  1:3830907/1‑150 (MQ=255)
                                                    tttGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGaaa                                                                                                 <  1:4479638/150‑1 (MQ=255)
                                                      tGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAgg                                                                                               <  2:3877761/150‑1 (MQ=255)
                                                      tGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAgg                                                                                               <  1:1051537/150‑1 (MQ=255)
                                                      tGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTGTACGAGAAGGGCGCGGACGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAgg                                                                                               >  2:2641320/1‑150 (MQ=255)
                                                       gCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAggg                                                                                              >  1:3454103/1‑150 (MQ=255)
                                                          gaagaCGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATg                                                                                           <  2:4457388/150‑1 (MQ=255)
                                                            agaCGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCa                                                                                         <  1:2124167/150‑1 (MQ=255)
                                                             gaCGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAg                                                                                        <  1:4402715/150‑1 (MQ=255)
                                                              aCGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGc                                                                                       >  1:2637389/1‑150 (MQ=255)
                                                               cGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCt                                                                                      >  1:82392/1‑150 (MQ=255)
                                                                gCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCtt                                                                                     >  1:2431233/1‑150 (MQ=255)
                                                                gCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCtt                                                                                     >  1:128204/1‑150 (MQ=255)
                                                                   aGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTAt                                                                                  <  1:2701636/150‑1 (MQ=255)
                                                                        gATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGa                                                                             >  2:4320688/1‑150 (MQ=255)
                                                                          tGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTGCGAGc                                                                           >  2:2855396/1‑150 (MQ=255)
                                                                           ggCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCg                                                                          <  1:206299/150‑1 (MQ=255)
                                                                            gcgcACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGt                                                                         <  1:1564461/150‑1 (MQ=255)
                                                                              gcACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCa                                                                       >  1:4037869/1‑150 (MQ=255)
                                                                              gcACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCa                                                                       <  1:82319/150‑1 (MQ=255)
                                                                               cACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCat                                                                      <  2:4251200/150‑1 (MQ=255)
                                                                               cACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCat                                                                      <  1:2817397/150‑1 (MQ=255)
                                                                                   cGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatgatg                                                                  <  2:1433438/150‑1 (MQ=255)
                                                                                   cGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatgatg                                                                  >  2:1568759/1‑150 (MQ=255)
                                                                                    gATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATgg                                                                 <  2:554927/150‑1 (MQ=255)
                                                                                     aTCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGt                                                                <  2:3088685/150‑1 (MQ=255)
                                                                                       ccGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAg                                                              <  2:3912486/150‑1 (MQ=255)
                                                                                        cGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGt                                                             >  1:1436280/1‑150 (MQ=255)
                                                                                         gCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTg                                                            >  1:685479/1‑150 (MQ=255)
                                                                                              gATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCg                                                       <  2:1505821/150‑1 (MQ=255)
                                                                                               atatGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGa                                                      >  1:3523205/1‑150 (MQ=255)
                                                                                                  tGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCt                                                   <  1:853003/150‑1 (MQ=255)
                                                                                                  tGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCt                                                   >  2:3157638/1‑150 (MQ=255)
                                                                                                  tGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCt                                                   <  2:2569451/150‑1 (MQ=255)
                                                                                                      cATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                                               <  1:3137851/150‑1 (MQ=255)
                                                                                                        ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGGGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacg                                             >  1:4004023/1‑150 (MQ=255)
                                                                                                         tGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacga                                            <  2:1890474/150‑1 (MQ=255)
                                                                                                         tGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacga                                            <  2:3658713/150‑1 (MQ=255)
                                                                                                           agaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACt                                          <  1:4516558/150‑1 (MQ=255)
                                                                                                           agaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACt                                          >  1:4435545/1‑150 (MQ=255)
                                                                                                           agaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACt                                          >  2:2445369/1‑150 (MQ=255)
                                                                                                           agaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACt                                          <  1:1658233/150‑1 (MQ=255)
                                                                                                            gaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACtt                                         <  2:1066824/150‑1 (MQ=255)
                                                                                                             aTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACttt                                        <  1:3198285/150‑1 (MQ=255)
                                                                                                                 acaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGc                                    >  2:356822/1‑150 (MQ=255)
                                                                                                                 acaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGc                                    <  1:1873496/150‑1 (MQ=255)
                                                                                                                  caacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCa                                   <  1:1741938/150‑1 (MQ=255)
                                                                                                                  caacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCa                                   <  2:4040817/150‑1 (MQ=255)
                                                                                                                   aacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAg                                  <  1:1385109/150‑1 (MQ=255)
                                                                                                                    acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAgg                                 >  2:3745521/1‑150 (MQ=255)
                                                                                                                    acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAgg                                 >  2:2426665/1‑150 (MQ=255)
                                                                                                                         tACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATg                            >  1:1149833/1‑150 (MQ=255)
                                                                                                                            acCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa                         >  1:550775/1‑150 (MQ=255)
                                                                                                                            acCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTAGTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGaa                         >  1:714319/1‑150 (MQ=255)
                                                                                                                                   ccGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGt                  >  1:1346912/1‑150 (MQ=255)
                                                                                                                                       ttACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGaa              <  2:1335616/150‑1 (MQ=255)
                                                                                                                                        tACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAAt             >  2:3360318/1‑150 (MQ=255)
                                                                                                                                          cGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGt           <  1:3352075/150‑1 (MQ=255)
                                                                                                                                          cGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGt           >  2:2242172/1‑150 (MQ=255)
                                                                                                                                          cGAGAAGGGAGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGt           <  1:2832984/150‑1 (MQ=255)
                                                                                                                                           gagaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTCGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTc          >  2:1112511/1‑150 (MQ=255)
                                                                                                                                           gagaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTc          <  1:4100615/150‑1 (MQ=255)
                                                                                                                                             gaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGa        <  1:923230/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTCCCGCGCAGTTAACCGCATCAATAATGTACGCACCATGCGCGGATTGCAGTTTGCAGAAGACGCCAGCCCGATGGCGCACCCGATCCGCCCGGATATGGTCATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGTCGAATGTCGATCTCTC  >  CP009273/987082‑987378

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: