Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,165,425 C→A 5.4% D136Y (GAC→TAC)  ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,165,4250CA 5.4% 1189.8 / 26.7 296D136Y (GAC→TAC) ldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (118/162);  new base A (5/11);  total (123/173)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 304 total aligned reads displayed.

GCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGC  >  CP009273/1165276‑1165574
                                                                                                                                                     |                                                                                                                                                     
   cATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTcacc                                                                                                                                                    <  1:1348950/150‑1 (MQ=255)
      tGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCaccacc                                                                                                                                                 <  1:3192120/150‑1 (MQ=255)
      tGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAaccacc                                                                                                                                                 <  1:1703572/150‑1 (MQ=255)
      tGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAaccacc                                                                                                                                                 <  1:2623199/150‑1 (MQ=255)
       gggTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCt                                                                                                                                                >  1:43589/1‑150 (MQ=255)
       gggTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCt                                                                                                                                                >  2:4090020/1‑150 (MQ=255)
          ttATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAAc                                                                                                                                             <  1:635166/150‑1 (MQ=255)
           tATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACt                                                                                                                                            <  1:1442091/150‑1 (MQ=255)
            aTCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTg                                                                                                                                           <  1:4945765/150‑1 (MQ=255)
            aTCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTg                                                                                                                                           <  2:1230140/150‑1 (MQ=255)
            aTCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTg                                                                                                                                           <  1:3014160/150‑1 (MQ=255)
             tccagtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGa                                                                                                                                          <  1:4507209/150‑1 (MQ=255)
              ccagtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGAc                                                                                                                                         >  1:1567490/1‑150 (MQ=255)
              ccagtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGAc                                                                                                                                         <  1:993596/150‑1 (MQ=255)
                agtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCa                                                                                                                                       >  1:1008106/1‑150 (MQ=255)
                agtccagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCa                                                                                                                                       <  1:1721043/150‑1 (MQ=255)
                    cagCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATAc                                                                                                                                   >  1:4608873/1‑150 (MQ=255)
                       ccGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTa                                                                                                                                >  1:3269439/1‑150 (MQ=255)
                        cGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTAt                                                                                                                               <  1:1146627/150‑1 (MQ=255)
                        cGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTAt                                                                                                                               >  2:1427827/1‑150 (MQ=255)
                          gCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTg                                                                                                                             <  2:3026830/150‑1 (MQ=255)
                          gCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTg                                                                                                                             <  2:1385894/150‑1 (MQ=255)
                           cACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTgg                                                                                                                            <  2:3457651/150‑1 (MQ=255)
                           cACTACCGCAGGCAACTCAAGTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTgg                                                                                                                            <  1:2553880/150‑1 (MQ=255)
                            aCTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTCTTGGa                                                                                                                           >  2:1594246/1‑150 (MQ=255)
                              tACCGCAGGCAACTCAAGTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATa                                                                                                                         <  2:3719636/150‑1 (MQ=255)
                                ccGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATaca                                                                                                                       <  1:2893839/150‑1 (MQ=255)
                                 cGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATacac                                                                                                                      <  1:338195/150‑1 (MQ=255)
                                   cAGGCAACTCAAGTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACgg                                                                                                                    <  2:3473776/150‑1 (MQ=255)
                                     ggCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGtt                                                                                                                  <  1:2590548/150‑1 (MQ=255)
                                       cAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTAc                                                                                                                <  2:2103916/150‑1 (MQ=255)
                                         aCTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCg                                                                                                              >  1:1832305/1‑150 (MQ=255)
                                         aCTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCg                                                                                                              >  1:443125/1‑150 (MQ=255)
                                           tCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGaa                                                                                                            >  1:3619259/1‑150 (MQ=255)
                                            cAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAt                                                                                                           <  2:892778/150‑1 (MQ=255)
                                             aaTTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAtt                                                                                                          <  1:4686706/150‑1 (MQ=255)
                                                tCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAttttt                                                                                                       <  2:1992243/150‑1 (MQ=255)
                                                tCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAAttttt                                                                                                       <  1:4815076/150‑1 (MQ=255)
                                                  ccTGTGCTTTATAACGGGCGCGGATGTTTGCCTTTGGCGTCCAGTTTTGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTAc                                                                                                     <  1:885859/150‑1 (MQ=255)
                                                     gtgCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcg                                                                                                  <  2:77894/150‑1 (MQ=255)
                                                       gCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggc                                                                                                <  2:3832926/150‑1 (MQ=255)
                                                        cTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggcg                                                                                               >  2:234333/1‑150 (MQ=255)
                                                        cTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCcggcg                                                                                               <  1:1120195/150‑1 (MQ=255)
                                                          ttATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCggg                                                                                             >  2:662160/1‑150 (MQ=255)
                                                           tataACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgt                                                                                            >  1:480328/1‑150 (MQ=255)
                                                           tataACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGGTGGGTTTGCACGTTTGTATGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGTATACACGGTTACCGAATTTTTACCCGGCGGgt                                                                                            >  1:2417683/1‑150 (MQ=255)
                                                            ataACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgta                                                                                           >  2:4306801/1‑150 (MQ=255)
                                                            ataACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGgta                                                                                           <  2:2189151/150‑1 (MQ=255)
                                                                cGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTaa                                                                                       >  1:3620516/1‑150 (MQ=255)
                                                                cGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTaa                                                                                       >  2:2224358/1‑150 (MQ=255)
                                                                   gcgcGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATag                                                                                    <  2:2783835/150‑1 (MQ=255)
                                                                    cgcgGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATaga                                                                                   >  1:850693/1‑150 (MQ=255)
                                                                    cgcgGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATaga                                                                                   <  2:1512620/150‑1 (MQ=255)
                                                                       ggATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTCGTAATAGAGAc                                                                                >  2:606728/1‑150 (MQ=255)
                                                                       ggATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGAc                                                                                <  2:1090008/150‑1 (MQ=255)
                                                                          tGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCa                                                                             <  2:4387216/150‑1 (MQ=255)
                                                                            tttGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGc                                                                           <  1:504714/150‑1 (MQ=255)
                                                                               gCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTcc                                                                        >  1:2405565/1‑150 (MQ=255)
                                                                               gCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTcc                                                                        >  2:3489937/1‑150 (MQ=255)
                                                                               gCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTcc                                                                        <  2:4205841/150‑1 (MQ=255)
                                                                               gCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTcc                                                                        <  1:1310212/150‑1 (MQ=255)
                                                                                    tGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAAt                                                                   <  2:4596461/150‑1 (MQ=255)
                                                                                       cGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGtt                                                                <  2:183156/150‑1 (MQ=255)
                                                                                        gTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTg                                                               >  2:2268504/1‑150 (MQ=255)
                                                                                           cAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATc                                                            >  1:3767564/1‑150 (MQ=255)
                                                                                           cAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATc                                                            <  2:3562626/150‑1 (MQ=255)
                                                                                                tGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTAACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAAt                                                       <  1:3575816/150‑1 (MQ=255)
                                                                                                    tttGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCt                                                   <  1:593848/150‑1 (MQ=255)
                                                                                                      tGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTc                                                 >  1:2491970/1‑150 (MQ=255)
                                                                                                      tGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTc                                                 >  1:2336634/1‑150 (MQ=255)
                                                                                                        cACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTcgc                                               <  1:2765019/150‑1 (MQ=255)
                                                                                                            tttGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCgcgg                                           >  2:2887165/1‑150 (MQ=255)
                                                                                                              tGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCgcggcg                                         <  2:4349525/150‑1 (MQ=255)
                                                                                                              tGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCgcggcg                                         >  1:3489224/1‑150 (MQ=255)
                                                                                                                 cTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCAt                                      >  1:3880447/1‑150 (MQ=255)
                                                                                                                  tGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATc                                     >  2:243081/1‑150 (MQ=255)
                                                                                                                   gAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCt                                    >  2:2727324/1‑150 (MQ=255)
                                                                                                                   gAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCt                                    <  2:1531287/150‑1 (MQ=255)
                                                                                                                    aaaCGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTg                                   <  2:3495487/150‑1 (MQ=255)
                                                                                                                        ggtggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAg                               <  2:514658/150‑1 (MQ=255)
                                                                                                                          tggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTa                             >  1:2011869/1‑150 (MQ=255)
                                                                                                                           ggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAg                            <  1:2264892/150‑1 (MQ=255)
                                                                                                                           ggtCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAg                            <  2:1684041/150‑1 (MQ=255)
                                                                                                                             tCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAggg                          >  1:1472187/1‑150 (MQ=255)
                                                                                                                             tCACCATTGTCGGTGACAGAGTGTCAACACATAACTTACCAATACCTATTGGATACACGGTTACCTAATTTTTAACCGGCGGGTAGTAATAGAGAAGCAGCTCCGCAATGTTGATAAAAAATACTTCGCGCGGCGCATCTTGAAGTAggg                          >  2:4038881/1‑150 (MQ=255)
                                                                                                                              cACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGt                         <  2:76168/150‑1 (MQ=255)
                                                                                                                                 cATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTg                      >  2:3433133/1‑150 (MQ=255)
                                                                                                                                  aTTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGc                     >  2:3934859/1‑150 (MQ=255)
                                                                                                                                  aTTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGc                     >  2:5049033/1‑150 (MQ=255)
                                                                                                                                    tGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCaa                   >  1:1735736/1‑150 (MQ=255)
                                                                                                                                    tGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCaa                   <  1:3956173/150‑1 (MQ=255)
                                                                                                                                    tGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATATTTACCCGGCGGGTAGTAATATAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCaa                   >  2:4742315/1‑150 (MQ=255)
                                                                                                                                        ggTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACggg               >  2:3288787/1‑150 (MQ=255)
                                                                                                                                         gtgtCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGa              <  1:4849198/150‑1 (MQ=255)
                                                                                                                                           gtCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATc            <  1:2224358/150‑1 (MQ=255)
                                                                                                                                           gtCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATc            <  2:2571743/150‑1 (MQ=255)
                                                                                                                                            tCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCg           <  1:4272968/150‑1 (MQ=255)
                                                                                                                                                  tgtCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAaca     <  2:123567/150‑1 (MQ=255)
                                                                                                                                                    tCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACg   >  1:4235217/1‑150 (MQ=255)
                                                                                                                                                     caccacCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGc  >  1:3286777/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCCATTGGGTTATCCAGTCCAGCCGGCACTACCGCAGGCAACTCAATTCCCTGTGCTTTATAACGGGCGCGGATGTTTGCCGTTGGCGTCCAGGTTGGGTTTGCACGTTTGTCTGAAACGGTGGTCACCATTGTCGGTGTCAGCGTGTCACCACCTAACTGACCAATACCTATTGGATACACGGTTACCGAATTTTTACCCGGCGGGTAGTAATAGAGACGCAGCTCCGCAATGTTGATCACAATGCCTTCGCGCGGCGCATCTGGAAGTAGGGTTTGCAACGGGATCGTTAACACGC  >  CP009273/1165276‑1165574

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: