Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 97,521 C→A 100% A90A (GCC→GCA murC → UDP‑N‑acetylmuramate:L‑alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927397,5210CA96.9% 876.0 / 15.9 261A90A (GCC→GCAmurCUDP‑N‑acetylmuramate:L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base A (121/132);  total (126/135)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 261 total aligned reads displayed.

TGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCA  >  CP009273/97374‑97666
                                                                                                                                                   |                                                                                                                                                 
tGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGa                                                                                                                                                 >  1:688918/1‑150 (MQ=255)
     gTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAcc                                                                                                                                            >  1:4170038/1‑150 (MQ=255)
      ttATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAccc                                                                                                                                           <  2:4229998/150‑1 (MQ=255)
      ttATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAAccc                                                                                                                                           <  2:2660626/150‑1 (MQ=255)
             atcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAtt                                                                                                                                    <  1:4068425/150‑1 (MQ=255)
             atcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAtt                                                                                                                                    <  2:3703155/150‑1 (MQ=255)
             atcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAAtt                                                                                                                                    <  2:262603/150‑1 (MQ=255)
              tcagTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTg                                                                                                                                   <  1:3746259/150‑1 (MQ=255)
                    gTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgccg                                                                                                                             <  1:1307011/150‑1 (MQ=255)
                    gTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTcgccg                                                                                                                             >  1:3223380/1‑150 (MQ=255)
                         gATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAt                                                                                                                        >  1:269577/1‑150 (MQ=255)
                         gATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCAt                                                                                                                        >  1:3052282/1‑150 (MQ=255)
                              aGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAgc                                                                                                                   >  2:1202591/1‑150 (MQ=255)
                              aGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAgc                                                                                                                   <  2:3127914/150‑1 (MQ=255)
                               gcgcCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAgcg                                                                                                                  <  1:2505274/150‑1 (MQ=255)
                                   cAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTa                                                                                                              >  1:3480537/1‑150 (MQ=255)
                                    aaaTCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAt                                                                                                             >  1:469962/1‑150 (MQ=255)
                                       tCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTcc                                                                                                          <  2:1908840/150‑1 (MQ=255)
                                       tCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTcc                                                                                                          <  2:3976278/150‑1 (MQ=255)
                                       tCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTcc                                                                                                          >  1:1485663/1‑150 (MQ=255)
                                        ccGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCg                                                                                                         <  1:2886780/150‑1 (MQ=255)
                                        ccGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCg                                                                                                         <  1:3814756/150‑1 (MQ=255)
                                        ccGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCg                                                                                                         <  2:4163881/150‑1 (MQ=255)
                                          ggTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGt                                                                                                       >  2:4283943/1‑150 (MQ=255)
                                           gTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTg                                                                                                      >  2:4355836/1‑150 (MQ=255)
                                            tCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGa                                                                                                     <  2:2090513/150‑1 (MQ=255)
                                             cACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGAt                                                                                                    <  1:721207/150‑1 (MQ=255)
                                               cGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATcc                                                                                                  <  2:138675/150‑1 (MQ=255)
                                                gcagcaGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcg                                                                                                 <  1:3967568/150‑1 (MQ=255)
                                                  agcagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtc                                                                                               >  1:17383/1‑150 (MQ=255)
                                                    cagTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCcgtcgt                                                                                             >  2:1780634/1‑150 (MQ=255)
                                                     agTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTg                                                                                            >  2:2054254/1‑150 (MQ=255)
                                                       ttAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTgtc                                                                                          >  2:3362107/1‑148 (MQ=255)
                                                       ttAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGcc                                                                                          >  2:3742251/1‑150 (MQ=255)
                                                         aaTGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGa                                                                                        <  1:2966491/150‑1 (MQ=255)
                                                          aTGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGaa                                                                                       >  1:984038/1‑150 (MQ=255)
                                                           tGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGaaa                                                                                      >  1:4531335/1‑150 (MQ=255)
                                                            gAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAAt                                                                                     >  2:3324126/1‑150 (MQ=255)
                                                             aaTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATg                                                                                    >  1:50319/1‑150 (MQ=255)
                                                             aaTCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATg                                                                                    <  2:3376142/150‑1 (MQ=255)
                                                               tCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgct                                                                                  <  1:623158/150‑1 (MQ=255)
                                                                cTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctg                                                                                 <  2:560576/150‑1 (MQ=255)
                                                                 tGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctgg                                                                                <  2:3947336/150‑1 (MQ=255)
                                                                 tGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctgg                                                                                >  1:1864870/1‑150 (MQ=255)
                                                                  gggTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggc                                                                               >  2:43779/1‑150 (MQ=255)
                                                                    gTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggctg                                                                             <  2:3842830/150‑1 (MQ=255)
                                                                    gTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATgctggctg                                                                             <  2:1366441/150‑1 (MQ=255)
                                                                      gCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGGGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAg                                                                           <  1:3560611/150‑1 (MQ=255)
                                                                        gacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                         <  1:3263066/150‑1 (MQ=255)
                                                                        gacgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGtt                                                                         <  2:3863991/150‑1 (MQ=255)
                                                                         acgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTa                                                                        <  2:3784684/150‑1 (MQ=255)
                                                                          cgaTTTATTTCAACCCTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTaa                                                                       <  1:1560693/150‑1 (MQ=255)
                                                                          cgaTTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACACGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTaa                                                                       <  1:1253535/150‑1 (MQ=255)
                                                                             tttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGc                                                                    >  1:283840/1‑150 (MQ=255)
                                                                              ttatttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCg                                                                   <  1:3303503/150‑1 (MQ=255)
                                                                                atttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtt                                                                 >  2:149302/1‑150 (MQ=255)
                                                                                 tttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGttt                                                                >  2:1956078/1‑150 (MQ=255)
                                                                                  ttCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGtttt                                                               <  2:2504703/150‑1 (MQ=255)
                                                                                    cAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCg                                                             <  1:1336947/150‑1 (MQ=255)
                                                                                        cATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCAt                                                         >  1:817923/1‑150 (MQ=255)
                                                                                         aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                        >  1:1505688/1‑150 (MQ=255)
                                                                                         aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                        <  2:2728556/150‑1 (MQ=255)
                                                                                         aTCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATg                                                        <  1:2247993/150‑1 (MQ=255)
                                                                                           cGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGc                                                      <  1:3201600/150‑1 (MQ=255)
                                                                                            gCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCa                                                     <  2:2777793/150‑1 (MQ=255)
                                                                                               cGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCg                                                  <  1:4223910/150‑1 (MQ=255)
                                                                                                  aaaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                               >  2:2360278/1‑150 (MQ=255)
                                                                                                  aaaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                               >  2:4439684/1‑150 (MQ=255)
                                                                                                  aaaaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCa                                               >  2:4198269/1‑150 (MQ=255)
                                                                                                    accgtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAtt                                             >  1:2278893/3‑150 (MQ=255)
                                                                                                    aaCGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCAtt                                             <  2:582957/150‑1 (MQ=255)
                                                                                                     acgtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTg                                            >  1:490030/1‑150 (MQ=255)
                                                                                                      cgtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGc                                           <  1:2053643/150‑1 (MQ=255)
                                                                                                       gtacgtGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGAGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGcc                                          >  2:701593/1‑150 (MQ=255)
                                                                                                            tGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAAc                                     <  1:2533114/150‑1 (MQ=255)
                                                                                                             gATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACg                                    >  1:3272412/1‑150 (MQ=255)
                                                                                                              aTGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGc                                   >  1:889697/1‑150 (MQ=255)
                                                                                                                gCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCAc                                 >  1:2718090/1‑150 (MQ=255)
                                                                                                                gCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCAc                                 >  2:1315014/1‑150 (MQ=255)
                                                                                                                gCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTAATCCGGTGATCCGTCGTGCCCAAAAGCTGGCAGAGTTAATGCGTTTTCGTCATGGCATCGCCATAGCCGGAACGCAc                                 >  2:4255625/1‑150 (MQ=255)
                                                                                                                gCCAGCGTGGTCGTTGTTTCCAGCGCGAGTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCAc                                 <  2:1636221/150‑1 (MQ=255)
                                                                                                                  cAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACgg                               <  2:4502866/150‑1 (MQ=255)
                                                                                                                  cAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACgg                               >  2:3267720/1‑150 (MQ=255)
                                                                                                                   aGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGc                              <  2:1720477/150‑1 (MQ=255)
                                                                                                                       tGGTCGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCaaaa                          >  2:3484332/1‑150 (MQ=255)
                                                                                                                           cGTTGTTTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGAc                      <  1:1082650/150‑1 (MQ=255)
                                                                                                                            gttgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACa                     >  1:1729310/1‑150 (MQ=255)
                                                                                                                             ttgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACaa                    >  2:4255268/1‑150 (MQ=255)
                                                                                                                              tgttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAac                   >  1:4363640/1‑150 (MQ=255)
                                                                                                                               gttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacc                  <  1:2326156/150‑1 (MQ=255)
                                                                                                                               gttTCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacc                  >  2:3020215/1‑150 (MQ=255)
                                                                                                                                tttCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacca                 <  1:4254445/150‑1 (MQ=255)
                                                                                                                                tttCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAacca                 >  2:4439502/1‑150 (MQ=255)
                                                                                                                                  tCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc               >  2:2024387/1‑150 (MQ=255)
                                                                                                                                  tCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc               >  1:1986387/1‑150 (MQ=255)
                                                                                                                                  tCCAGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTAATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAaccacc               <  2:1588080/150‑1 (MQ=255)
                                                                                                                                     aGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgcg            >  1:514255/1‑150 (MQ=255)
                                                                                                                                     aGCGCGATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACcgcg            >  2:2485110/1‑150 (MQ=255)
                                                                                                                                       cgcgATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGAt          <  2:2870547/150‑1 (MQ=255)
                                                                                                                                         cgATTTCTGCAGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATgg        >  2:4072559/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TGAAGGTTATCAGATCAGTGGTTCCGATTTAGCGCCAAATCCGGTCACGCAGCAGTTAATGAATCTGGGTGCGACGATTTATTTCAACCATCGCCCGGAAAACGTACGTGATGCCAGCGTGGTCGTTGTTTCCAGCGCGATTTCTGCCGATAACCCGGAAATTGTCGCCGCTCATGAAGCGCGTATTCCGGTGATCCGTCGTGCCGAAATGCTGGCTGAGTTAATGCGTTTTCGTCATGGCATCGCCATTGCCGGAACGCACGGCAAAACGACAACCACCGCGATGGTTTCCA  >  CP009273/97374‑97666

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: