Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,907,967 C→A 21.9% L369F (TTG→TTT prc ← carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,907,9670CA21.9% 799.0 / 184.2 301L369F (TTG→TTTprccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base C (119/116);  new base A (37/29);  total (156/145)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 304 total aligned reads displayed.

ACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTACCAGCAGG  >  CP009273/1907818‑1908113
                                                                                                                                                     |                                                                                                                                                  
   tGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCaaa                                                                                                                                                 <  1:2802435/150‑1 (MQ=255)
        aaTGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCcac                                                                                                                                            <  2:3017151/150‑1 (MQ=255)
           gggACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGCCTGTCAAACCCACATa                                                                                                                                         >  2:3545536/1‑150 (MQ=255)
            ggACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAg                                                                                                                                        <  1:738900/150‑1 (MQ=255)
            ggACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAg                                                                                                                                        >  2:3501047/1‑150 (MQ=255)
               cccGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAg                                                                                                                                     >  1:3242060/1‑150 (MQ=255)
               cccGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAg                                                                                                                                     >  2:865539/1‑150 (MQ=255)
                 cGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGcc                                                                                                                                   <  2:3155631/150‑1 (MQ=255)
                 cGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGcc                                                                                                                                   >  2:3192635/1‑150 (MQ=255)
                  gCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGccc                                                                                                                                  >  1:1946901/1‑150 (MQ=255)
                   cAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCg                                                                                                                                 <  1:893228/150‑1 (MQ=255)
                   cAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCg                                                                                                                                 >  2:1836823/1‑150 (MQ=255)
                    aGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCgg                                                                                                                                <  2:2423040/150‑1 (MQ=255)
                      gAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGaa                                                                                                                              <  2:2649970/150‑1 (MQ=255)
                        aTAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAata                                                                                                                            <  2:3862133/150‑1 (MQ=255)
                        aTAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAata                                                                                                                            <  1:3142098/150‑1 (MQ=255)
                        aTAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAata                                                                                                                            <  2:590662/150‑1 (MQ=255)
                         tAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAatat                                                                                                                           >  1:711028/1‑150 (MQ=255)
                             cAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAg                                                                                                                       >  1:1741165/1‑150 (MQ=255)
                              aGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGc                                                                                                                      >  1:578707/1‑150 (MQ=255)
                              aGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGc                                                                                                                      <  2:4005557/150‑1 (MQ=255)
                               gACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCa                                                                                                                     >  2:1611239/1‑150 (MQ=255)
                               gACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTCCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCa                                                                                                                     <  2:3133948/150‑1 (MQ=255)
                                aCCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCAc                                                                                                                    <  1:2973057/150‑1 (MQ=255)
                                  cGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGc                                                                                                                  >  1:2688553/1‑150 (MQ=255)
                                   ggAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGcc                                                                                                                 <  1:1067070/150‑1 (MQ=255)
                                   ggAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGcc                                                                                                                 >  2:2370070/1‑150 (MQ=255)
                                   ggAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGcc                                                                                                                 >  2:2503952/1‑150 (MQ=255)
                                    gagaGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCg                                                                                                                >  1:4161717/1‑150 (MQ=255)
                                       agCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACt                                                                                                             >  1:1145032/1‑150 (MQ=255)
                                         cGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACttt                                                                                                           <  1:4664815/150‑1 (MQ=255)
                                         cGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACttt                                                                                                           <  2:4474200/150‑1 (MQ=255)
                                          gATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTtc                                                                                                          >  2:1820959/1‑150 (MQ=255)
                                          gATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTtc                                                                                                          >  2:4644887/1‑150 (MQ=255)
                                           aTACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTtct                                                                                                         >  2:3669517/1‑150 (MQ=255)
                                              cGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCtt                                                                                                      <  1:4384467/150‑1 (MQ=255)
                                               ggCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCttt                                                                                                     <  1:378754/150‑1 (MQ=255)
                                                gCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTa                                                                                                    >  1:4280740/1‑150 (MQ=255)
                                                gCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTa                                                                                                    >  1:1806298/1‑150 (MQ=255)
                                                   tCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCg                                                                                                 <  1:4763854/150‑1 (MQ=255)
                                                    cAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGa                                                                                                >  1:1841500/1‑150 (MQ=255)
                                                     aGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGAc                                                                                               <  2:4287455/150‑1 (MQ=255)
                                                              cccACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCAc                                                                                      <  1:2494615/150‑1 (MQ=255)
                                                              cccACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCAc                                                                                      <  1:2493514/150‑1 (MQ=255)
                                                              cccACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCAc                                                                                      >  2:211494/1‑150 (MQ=255)
                                                                  ccGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGAc                                                                                  <  2:3949013/150‑1 (MQ=255)
                                                                     ccATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACAtt                                                                               <  1:4446478/150‑1 (MQ=255)
                                                                     ccATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACAtt                                                                               <  1:4617729/150‑1 (MQ=255)
                                                                     ccATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACAtt                                                                               <  1:3282501/150‑1 (MQ=255)
                                                                      cATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACAttt                                                                              <  2:4398078/150‑1 (MQ=255)
                                                                       aTTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACAtttt                                                                             >  1:2465756/1‑150 (MQ=255)
                                                                            tACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACcg                                                                        <  1:1293699/150‑1 (MQ=255)
                                                                            tACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACcg                                                                        <  1:1149045/150‑1 (MQ=255)
                                                                            tACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACcg                                                                        <  1:2753031/150‑1 (MQ=255)
                                                                                  ggTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGt                                                                  <  2:2266333/150‑1 (MQ=255)
                                                                                    tCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCt                                                                >  2:2197927/1‑150 (MQ=255)
                                                                                     cGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCtt                                                               <  1:4826170/150‑1 (MQ=255)
                                                                                     cGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCtt                                                               <  1:3433393/150‑1 (MQ=255)
                                                                                      gatgatGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTc                                                              <  2:4767989/150‑1 (MQ=255)
                                                                                       atgatgACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCg                                                             >  1:3575657/1‑150 (MQ=255)
                                                                                         gatgaCGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCgag                                                           <  1:2180565/150‑1 (MQ=255)
                                                                                           tgaCGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGACTATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGAc                                                         >  2:4025656/1‑150 (MQ=255)
                                                                                                        aTTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACg                                            >  2:4067900/1‑150 (MQ=255)
                                                                                                             gTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCa                                       <  2:2779148/150‑1 (MQ=255)
                                                                                                                  tCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGtt                                  <  1:2869828/150‑1 (MQ=255)
                                                                                                                   ccAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTa                                 <  1:1608824/150‑1 (MQ=255)
                                                                                                                    cAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTAc                                <  1:21565/150‑1 (MQ=255)
                                                                                                                      gTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAg                              >  2:811060/1‑150 (MQ=255)
                                                                                                                      gTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAg                              >  2:2098073/1‑150 (MQ=255)
                                                                                                                      gTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAg                              >  2:671937/1‑150 (MQ=255)
                                                                                                                       tttCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGt                             >  1:3311769/1‑150 (MQ=255)
                                                                                                                       tttCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGt                             >  1:304951/1‑150 (MQ=255)
                                                                                                                          cTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACg                          <  1:3917732/150‑1 (MQ=255)
                                                                                                                          cTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACg                          <  2:1046658/150‑1 (MQ=255)
                                                                                                                           tGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACgg                         <  1:3501047/150‑1 (MQ=255)
                                                                                                                            gCAGTTGCACTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACggg                        <  1:2236321/150‑1 (MQ=255)
                                                                                                                             cAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGt                       <  2:1436288/150‑1 (MQ=255)
                                                                                                                              aGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTc                      <  1:4110986/150‑1 (MQ=255)
                                                                                                                              aGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTc                      <  1:1920589/150‑1 (MQ=255)
                                                                                                                              aGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTc                      >  1:1437409/1‑150 (MQ=255)
                                                                                                                               gTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCt                     <  2:4144507/150‑1 (MQ=255)
                                                                                                                               gTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCt                     >  1:551174/1‑150 (MQ=255)
                                                                                                                               gTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCt                     >  2:2027710/1‑150 (MQ=255)
                                                                                                                                 tGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTg                   <  1:666962/150‑1 (MQ=255)
                                                                                                                                   cACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGt                 >  1:4306996/1‑150 (MQ=255)
                                                                                                                                    aCTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTc                >  1:2795642/1‑150 (MQ=255)
                                                                                                                                    aCTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTc                >  1:3819277/1‑150 (MQ=255)
                                                                                                                                    aCTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTc                >  1:365403/1‑150 (MQ=255)
                                                                                                                                    aCTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTc                <  1:406826/150‑1 (MQ=255)
                                                                                                                                    aCTTTGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTc                >  1:2269881/1‑150 (MQ=255)
                                                                                                                                       ttGACATCGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCt             <  2:4853023/150‑1 (MQ=255)
                                                                                                                                        tGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCtt            <  2:1946901/150‑1 (MQ=255)
                                                                                                                                          aCATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTa          <  1:865539/150‑1 (MQ=255)
                                                                                                                                           cATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTAc         >  1:3013879/1‑150 (MQ=255)
                                                                                                                                           cATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTAc         >  1:526347/1‑150 (MQ=255)
                                                                                                                                           cATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTAc         <  2:1917455/150‑1 (MQ=255)
                                                                                                                                            aTCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTAcc        <  1:2557620/150‑1 (MQ=255)
                                                                                                                                              cGTCTGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTGTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTACcag      <  2:3048226/150‑1 (MQ=255)
                                                                                                                                                 cTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTACcagcag   <  2:3107674/150‑1 (MQ=255)
                                                                                                                                                  tGTAAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTACCAGCAgg  <  1:3192635/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCCATTGCTACGCAGGTCGATGATGACGCTGCTGACATTCTGTTTTTCCAGTTTCTGCAGTTGCACTTTGACATCGTCTGTCAAACCCACATAGAAGCCCGGAATATCCAGCACGCCGACTTTCTCTTTACCGACGGTCTTCACCGACATTTTAACCGCGCGGTCTTCGAGACGAATACGTTCACGGGTCAACGTTACAGTACGGGTCTTGGTCCCTTTACCAGCAGG  >  CP009273/1907818‑1908113

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: