Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,181,504 G→T 100% D69Y (GAT→TAT)  pepT → peptidase T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,181,5040GT98.1% 988.2 / 7.4 310D69Y (GAT→TAT) pepTpeptidase T
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (148/156);  total (152/158)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 313 total aligned reads displayed.

GCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGTTC  >  CP009273/1181357‑1181650
                                                                                                                                                   |                                                                                                                                                  
 cAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCtata                                                                                                                                                 >  2:739511/1‑150 (MQ=255)
     ggTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATccc                                                                                                                                             <  2:2440340/150‑1 (MQ=255)
       tGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCgg                                                                                                                                           >  2:4719107/1‑150 (MQ=255)
        gagaCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGc                                                                                                                                          >  1:1468161/1‑150 (MQ=255)
          gaCAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGa                                                                                                                                        >  2:5023530/1‑150 (MQ=255)
            cAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAtt                                                                                                                                      <  2:2225081/150‑1 (MQ=255)
             aGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTg                                                                                                                                     >  2:1441742/1‑150 (MQ=255)
               gTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGc                                                                                                                                   >  2:279565/1‑150 (MQ=255)
                ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGGGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCt                                                                                                                                  >  1:3016968/1‑150 (MQ=255)
                ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCt                                                                                                                                  >  1:4356910/1‑150 (MQ=255)
                ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCt                                                                                                                                  >  2:498297/1‑150 (MQ=255)
                ttCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCt                                                                                                                                  >  2:2390357/1‑150 (MQ=255)
                  cccAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCttt                                                                                                                                >  2:3538455/1‑150 (MQ=255)
                    cAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAt                                                                                                                              <  2:3472548/150‑1 (MQ=255)
                    cAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAt                                                                                                                              >  2:1546465/1‑150 (MQ=255)
                     aGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAtt                                                                                                                             <  1:2056121/150‑1 (MQ=255)
                     aGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAtt                                                                                                                             >  1:2146840/1‑150 (MQ=255)
                      gCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTAttt                                                                                                                            <  2:1697528/150‑1 (MQ=255)
                       cACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtc                                                                                                                           >  1:3852933/1‑150 (MQ=255)
                         cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtctc                                                                                                                         >  2:3770971/1‑150 (MQ=255)
                         cGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTtctc                                                                                                                         <  2:2857150/150‑1 (MQ=255)
                           gAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAt                                                                                                                       <  1:2390203/150‑1 (MQ=255)
                            aaGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCTGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtg                                                                                                                      >  1:3593514/1‑150 (MQ=255)
                              ggCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAtgtg                                                                                                                    <  1:3278212/150‑1 (MQ=255)
                                   aTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAc                                                                                                               <  2:4263786/150‑1 (MQ=255)
                                   aTGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCATAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAc                                                                                                               >  1:1854740/1‑150 (MQ=255)
                                     ggAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCt                                                                                                             <  2:2151179/150‑1 (MQ=255)
                                     ggAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTAATGTGGATACCt                                                                                                             >  1:1020660/1‑150 (MQ=255)
                                      gAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTc                                                                                                            <  2:4192062/150‑1 (MQ=255)
                                       aaGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCa                                                                                                           <  2:851525/150‑1 (MQ=255)
                                        aGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCAc                                                                                                          >  2:2940761/1‑150 (MQ=255)
                                          ttATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAGTGGCTTTATTTCTCATGTGGATACCTCACCg                                                                                                        >  1:94754/1‑150 (MQ=255)
                                           tATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCgg                                                                                                       >  1:5130202/1‑150 (MQ=255)
                                                cATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGc                                                                                                  >  1:4411296/1‑150 (MQ=255)
                                                  tCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAg                                                                                                >  1:3685610/1‑150 (MQ=255)
                                                    tgctgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCg                                                                                              <  2:539871/150‑1 (MQ=255)
                                                       tgACAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCa                                                                                           <  2:4658374/150‑1 (MQ=255)
                                                       tgAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCa                                                                                           >  1:4490629/1‑150 (MQ=255)
                                                        gAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaa                                                                                          >  2:2651091/1‑150 (MQ=255)
                                                           agagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCaaaaa                                                                                       >  1:403728/1‑150 (MQ=255)
                                                            gagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAt                                                                                      >  1:2849018/1‑150 (MQ=255)
                                                            gagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAt                                                                                      <  2:4801275/150‑1 (MQ=255)
                                                            gagCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCATTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAt                                                                                      <  1:1281504/150‑1 (MQ=255)
                                                             agcagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtg                                                                                     >  1:2724606/1‑150 (MQ=255)
                                                             agcagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtg                                                                                     >  1:1581049/1‑150 (MQ=255)
                                                               cagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgtg                                                                                   <  1:1959938/150‑1 (MQ=255)
                                                               cagcTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAAtgtg                                                                                   <  1:3965173/150‑1 (MQ=255)
                                                                  cTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAAt                                                                                <  2:1237514/150‑1 (MQ=255)
                                                                    cGAAGAGATGGGGCTTATCAATGTGACCTTAAGTTAGAATGGCACTTTGATTGCGACGTTACCTGCTAAAGTCCCTGGATATATCCCGTCGATTGGCTTTATTTATCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                              <  2:3768191/150‑1 (MQ=255)
                                                                    cGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATcc                                                                              <  2:3638760/150‑1 (MQ=255)
                                                                     gAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCCAAAATGTGAATCCg                                                                             >  1:703301/1‑150 (MQ=255)
                                                                      aaGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGc                                                                            >  2:2111768/1‑150 (MQ=255)
                                                                       agagATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCa                                                                           >  2:1973247/1‑150 (MQ=255)
                                                                            tGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttg                                                                      <  1:4966729/150‑1 (MQ=255)
                                                                            tGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAACTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttg                                                                      >  2:1090895/1‑150 (MQ=255)
                                                                               ggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgttg                                                                   >  1:4559537/1‑150 (MQ=255)
                                                                               ggCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAAttgttg                                                                   <  1:3664745/150‑1 (MQ=255)
                                                                                gCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGa                                                                  <  2:5192370/150‑1 (MQ=255)
                                                                                  ttATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                                <  2:3346183/150‑1 (MQ=255)
                                                                                  ttATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGaaa                                                                <  2:3244532/150‑1 (MQ=255)
                                                                                    aTCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAAc                                                              <  2:389344/150‑1 (MQ=255)
                                                                                     tCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACt                                                             <  1:5191421/150‑1 (MQ=255)
                                                                                     tCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACt                                                             <  1:2464981/150‑1 (MQ=255)
                                                                                       aaTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTAt                                                           <  1:2894907/150‑1 (MQ=255)
                                                                                       aaTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGAGTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCCAAAATGTGAATCCGCAAATTGTTGAAAACAAt                                                           >  1:1199196/1‑150 (MQ=255)
                                                                                        aTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATc                                                          >  2:4077527/1‑150 (MQ=255)
                                                                                        aTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATc                                                          >  2:4015849/1‑150 (MQ=255)
                                                                                        aTGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATc                                                          <  1:3384334/150‑1 (MQ=255)
                                                                                         tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                                         <  1:4204183/150‑1 (MQ=255)
                                                                                         tgtgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcg                                                         <  1:54831/150‑1 (MQ=255)
                                                                                           tgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATcgcg                                                       >  1:1124625/1‑150 (MQ=255)
                                                                                           tgACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATACGCAAATTGTTGAAAACTATcgcg                                                       >  1:1372565/1‑150 (MQ=255)
                                                                                            gACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCgg                                                      >  2:2073427/1‑150 (MQ=255)
                                                                                               cTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTgg                                                   >  1:339705/1‑150 (MQ=255)
                                                                                                 tAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCg                                                 >  1:1905706/1‑150 (MQ=255)
                                                                                                      gagaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATAtt                                            <  1:4605886/150‑1 (MQ=255)
                                                                                                      gagaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATAtt                                            >  2:3164805/1‑150 (MQ=255)
                                                                                                       agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTg                                           >  1:4731638/1‑150 (MQ=255)
                                                                                                       agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTg                                           >  2:5256749/1‑150 (MQ=255)
                                                                                                       agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTg                                           >  2:2040606/1‑150 (MQ=255)
                                                                                                       agaAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTg                                           <  1:2085656/150‑1 (MQ=255)
                                                                                                           gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCt                                       <  2:77410/150‑1 (MQ=255)
                                                                                                           gggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCt                                       <  1:1222624/150‑1 (MQ=255)
                                                                                                            ggCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTg                                      >  1:3120587/1‑150 (MQ=255)
                                                                                                               ccTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGt                                   >  2:5102719/2‑150 (MQ=255)
                                                                                                               aCTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGt                                   <  1:2741746/150‑1 (MQ=255)
                                                                                                                    gATGGCGACGTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTTCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCGGGGTATCgg                              >  1:1767649/1‑150 (MQ=255)
                                                                                                                             gTTACCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATAACTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGAGATTGCGCTGGGTATCGGCGAGGAAgg                     >  1:4250782/1‑149 (MQ=255)
                                                                                                                                aCCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt                  <  2:4518089/150‑1 (MQ=255)
                                                                                                                                aCCGGCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGtttt                  <  2:3005835/150‑1 (MQ=255)
                                                                                                                                    gCTAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCAAGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCCAATTGTTGAAAACTATCGCGGTGGCGAGATTGCGCTGGGTATCGGCGATGAAGTTTTATCa              >  1:1997719/1‑150 (MQ=255)
                                                                                                                                      tAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCAcc            >  2:3798551/1‑150 (MQ=255)
                                                                                                                                      tAACGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCAcc            <  2:1854740/150‑1 (MQ=255)
                                                                                                                                       aaCGTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCg           <  1:4108063/150‑1 (MQ=255)
                                                                                                                                          gTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGtt        <  2:2702767/150‑1 (MQ=255)
                                                                                                                                          gTCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGtt        <  1:3928848/150‑1 (MQ=255)
                                                                                                                                           tCCCTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTa       <  1:88303/150‑1 (MQ=255)
                                                                                                                                            cccTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTAt      >  1:4823385/1‑150 (MQ=255)
                                                                                                                                              cTGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGt    >  2:4946597/1‑150 (MQ=255)
                                                                                                                                               tGGCTATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGtt   >  2:1778515/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCAGGGGTGAGACAGGTTCCCAGCACGGAAGGCCAATGGAAGTTATTGCATCTGCTGAAAGAGCAGCTCGAAGAGATGGGGCTTATCAATGTGACCTTAAGTGAGAAGGGCACTTTGATGGCGACGTTACCGGCTAACGTCCCTGGCGATATCCCGGCGATTGGCTTTATTTCTCATGTGGATACCTCACCGGATTGCAGCGGCAAAAATGTGAATCCGCAAATTGTTGAAAACTATCGCGGTGGCGATATTGCGCTGGGTATCGGCGATGAAGTTTTATCACCGGTTATGTTC  >  CP009273/1181357‑1181650

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: ‑