Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,322,609 C→A 100% W223L (TGG→TTG)  yfaQ ← tandem DUF2300 domain protein, putative host defense protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,322,6090CA98.0% 822.4 / 4.2 262W223L (TGG→TTG) yfaQtandem DUF2300 domain protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base A (117/139);  total (120/142)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 269 total aligned reads displayed.

GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACA  >  CP009273/2322462‑2322755
                                                                                                                                                   |                                                                                                                                                  
gCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCTCATAACGTGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaac                                                                                                                                                  >  1:4985210/1‑150 (MQ=255)
gCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaaca                                                                                                                                                 <  1:988685/150‑1 (MQ=255)
 cTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCaaca                                                                                                                                                 <  1:5097788/150‑1 (MQ=255)
  tgtgATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAg                                                                                                                                                <  2:700614/150‑1 (MQ=255)
      aTTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgt                                                                                                                                            >  1:46504/1‑150 (MQ=255)
      aTTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgt                                                                                                                                            <  2:3029884/150‑1 (MQ=255)
       ttttCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                           >  1:4506221/1‑150 (MQ=255)
       ttttCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGgtg                                                                                                                                           >  2:535509/1‑150 (MQ=255)
          tcttcAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTAc                                                                                                                                        <  2:2778483/150‑1 (MQ=255)
                gCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCgg                                                                                                                                  <  2:3082143/150‑1 (MQ=255)
                 cGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGt                                                                                                                                 >  2:220291/1‑150 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCACCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTCACCAACAGCGGTGTACCATCGGtt                                                                                                                                >  2:5186196/1‑150 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGtt                                                                                                                                <  2:1019941/150‑1 (MQ=255)
                  ggATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGtt                                                                                                                                >  2:2608790/1‑150 (MQ=255)
                   gATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTa                                                                                                                               <  1:3298520/150‑1 (MQ=255)
                      aGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAcc                                                                                                                            <  1:3785828/150‑1 (MQ=255)
                       gCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAcca                                                                                                                           <  2:812980/150‑1 (MQ=255)
                           gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                       <  2:1817879/150‑1 (MQ=255)
                           gCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAAccagcca                                                                                                                       >  1:2163027/1‑150 (MQ=255)
                             aaaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                     >  2:2006216/1‑150 (MQ=255)
                             aaaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                     >  1:827729/1‑150 (MQ=255)
                             aaaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                     >  1:55930/1‑150 (MQ=255)
                             aaaCAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACcagccagc                                                                                                                     >  1:1571789/1‑150 (MQ=255)
                              aaCAGCTCCACCTCTACACATAGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCa                                                                                                                    >  2:808890/1‑150 (MQ=255)
                                cAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCaaa                                                                                                                  >  2:3880387/1‑150 (MQ=255)
                                  gCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAac                                                                                                                <  1:4796128/150‑1 (MQ=255)
                                   cTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAacc                                                                                                               <  2:1011923/150‑1 (MQ=255)
                                   cTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAacc                                                                                                               >  2:4684085/1‑150 (MQ=255)
                                     ccACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccac                                                                                                             <  1:2854229/150‑1 (MQ=255)
                                     ccACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccac                                                                                                             >  2:426619/1‑150 (MQ=255)
                                      cACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAaccacc                                                                                                            >  2:306692/1‑150 (MQ=255)
                                       aCCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCAccc                                                                                                           <  2:2543978/150‑1 (MQ=255)
                                        ccTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCt                                                                                                          <  1:484807/150‑1 (MQ=255)
                                          tctACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGa                                                                                                        >  2:4060801/1‑150 (MQ=255)
                                              cacaTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGa                                                                                                    <  1:1812634/150‑1 (MQ=255)
                                                caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAtt                                                                                                  <  1:4163747/150‑1 (MQ=255)
                                                caTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAtt                                                                                                  <  2:287179/150‑1 (MQ=255)
                                                 aTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAttt                                                                                                 <  2:2304277/150‑1 (MQ=255)
                                                 aTTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGAttt                                                                                                 <  1:1890874/150‑1 (MQ=255)
                                                   tGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTcc                                                                                               <  2:4919548/150‑1 (MQ=255)
                                                     cccACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTg                                                                                             <  2:3520661/150‑1 (MQ=255)
                                                      ccACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGt                                                                                            >  2:3400852/1‑150 (MQ=255)
                                                        aCTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGttt                                                                                          <  1:862353/150‑1 (MQ=255)
                                                        aCTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGttt                                                                                          <  2:585642/150‑1 (MQ=255)
                                                         cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGTTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTa                                                                                         >  1:4255672/1‑150 (MQ=255)
                                                         cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTa                                                                                         >  2:671314/1‑150 (MQ=255)
                                                         cTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGGTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTa                                                                                         >  2:1509679/1‑150 (MQ=255)
                                                          tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTat                                                                                        <  2:3779139/150‑1 (MQ=255)
                                                          tGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTat                                                                                        <  2:4443560/150‑1 (MQ=255)
                                                             cAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcat                                                                                     <  2:1735843/150‑1 (MQ=255)
                                                              aaCCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTatcatc                                                                                    >  1:5291332/1‑150 (MQ=255)
                                                                ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCgg                                                                                  <  2:95008/150‑1 (MQ=255)
                                                                ccTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCgg                                                                                  >  2:4673588/1‑150 (MQ=255)
                                                                 cTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATACCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGc                                                                                 <  2:188624/150‑1 (MQ=255)
                                                                        ggCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAg                                                                          >  1:3033578/1‑150 (MQ=255)
                                                                         gCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGc                                                                         >  1:276100/1‑150 (MQ=255)
                                                                          cACTGGCAAAACCCGATTCAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCt                                                                        <  2:3342925/150‑1 (MQ=255)
                                                                          cACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCt                                                                        <  2:3628781/150‑1 (MQ=255)
                                                                          cACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCt                                                                        <  1:2495464/150‑1 (MQ=255)
                                                                             tGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCa                                                                     <  1:3658945/150‑1 (MQ=255)
                                                                              ggCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAg                                                                    <  2:4725166/150‑1 (MQ=255)
                                                                               gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAgg                                                                   <  1:2552895/150‑1 (MQ=255)
                                                                               gCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAgg                                                                   <  2:2438818/150‑1 (MQ=255)
                                                                                 aaaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtt                                                                 <  1:2819196/150‑1 (MQ=255)
                                                                                 aaaaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTACAGGtt                                                                 >  2:4306148/1‑150 (MQ=255)
                                                                                   aaCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGtttt                                                               >  1:1205065/1‑150 (MQ=255)
                                                                                    aCCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTc                                                              <  2:1259355/150‑1 (MQ=255)
                                                                                     cccGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCa                                                             <  2:4200336/150‑1 (MQ=255)
                                                                                        gATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCCTCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCCAACCCGCTCCCGGGGTTCAccc                                                          >  1:483945/1‑150 (MQ=255)
                                                                                           tAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTa                                                       <  1:2965729/150‑1 (MQ=255)
                                                                                                     cGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTTACCAAGAGCGGTGTACCATCGGTTATCCAGCGTGCAATACCACCCTGACGGATTTGCTGTTTATCATCGGCTAACCAGCTCGAGGTTTTCACCCGTACCCGACTAcc                                             <  2:85929/150‑1 (MQ=255)
                                                                                                        aTAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                          <  2:3154978/150‑1 (MQ=255)
                                                                                                        aTAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                          <  1:3897668/150‑1 (MQ=255)
                                                                                                        aTAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                          >  2:4411860/1‑150 (MQ=255)
                                                                                                        aTAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAg                                          >  1:798064/1‑150 (MQ=255)
                                                                                                           aCGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCAt                                       <  1:339640/150‑1 (MQ=255)
                                                                                                           aCGGGTTAGTACCGTTTTGCTGGAGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCAt                                       <  2:5051567/150‑1 (MQ=255)
                                                                                                            cGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATg                                      >  1:2345994/1‑150 (MQ=255)
                                                                                                            cGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATg                                      <  1:923472/150‑1 (MQ=255)
                                                                                                            cGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGAGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCAta                                      <  1:1675604/150‑2 (MQ=255)
                                                                                                             gggTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATgg                                     <  2:2750272/150‑1 (MQ=255)
                                                                                                                  aGTACCGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCg                                >  2:3958699/1‑150 (MQ=255)
                                                                                                                       cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                           <  1:1004330/150‑1 (MQ=255)
                                                                                                                       cGTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACAcc                           <  2:1466496/150‑1 (MQ=255)
                                                                                                                        gTTTTGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCg                          <  1:1980361/150‑1 (MQ=255)
                                                                                                                         ttttGCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGa                         >  1:1481004/1‑150 (MQ=255)
                                                                                                                             gCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAtttt                     <  2:2469777/150‑1 (MQ=255)
                                                                                                                             gCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAtttt                     <  2:1636988/150‑1 (MQ=255)
                                                                                                                             gCTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGAtttt                     <  2:5115446/150‑1 (MQ=255)
                                                                                                                              cTGGTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTc                    <  1:3880777/150‑1 (MQ=255)
                                                                                                                                ggTGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGc                  >  1:2660595/1‑150 (MQ=255)
                                                                                                                                  tGCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCt                <  2:46504/150‑1 (MQ=255)
                                                                                                                                   gCCACTGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTc               >  1:1626080/1‑150 (MQ=255)
                                                                                                                                        tGCCAGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCa          <  1:1400394/150‑1 (MQ=255)
                                                                                                                                            aGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac      >  1:5044091/1‑150 (MQ=255)
                                                                                                                                            aGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac      >  1:2065925/1‑150 (MQ=255)
                                                                                                                                            aGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac      >  2:3989564/1‑150 (MQ=255)
                                                                                                                                            aGTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcac      >  1:595089/1‑150 (MQ=255)
                                                                                                                                             gTAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcacc     >  2:4962582/1‑150 (MQ=255)
                                                                                                                                              tAACCAACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGcacca    >  2:4830133/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTGTGATTTTCTTCAGCGGATAGCGGGCAAACAGCTCCACCTCTACACATTGCCCACTGGCAACCTGCGTCGGCACTGGCAAAACCCGATTAAGCACCGTCGCATAACGGGTTAGTACCGTTTTGCTGGTGCCACTGCCAGTAACCCACAGCGGTGTACCATCGGTTAACCAGCCAGCAAAACCACCCTGACGGATTTCCTGTTTATCATCGGCAAACCAGCTCCAGGTTTTCACCCGTACCCGACTACCGAGCATGGATGCGACACCGATTTTCGCCTCGTCCAGCACCACA  >  CP009273/2322462‑2322755

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: