Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 729,456 G→A 100% R75R (AGG→AGA ybfB → putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273729,4560GA98.6% 1155.9 / 4.0 353R75R (AGG→AGAybfBputative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/5);  new base A (128/219);  total (128/225)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 353 total aligned reads displayed.

TTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGGATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATT  >  CP009273/729312‑729600
                                                                                                                                                |                                                                                                                                                
ttCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAAtt                                                                                                                                               >  1:664238/1‑148 (MQ=255)
 tCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTg                                                                                                                                              <  1:2901266/148‑1 (MQ=255)
   aatATTAGACTCCACTGGTTATGTGAGCGAGTTAATTAGTGTGTTGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGGATTGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTTAGAATTGAt                                                                                                                                            <  1:817072/146‑1 (MQ=255)
     tATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATag                                                                                                                                          <  2:144018/148‑1 (MQ=255)
     tATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATag                                                                                                                                          >  1:1969607/1‑148 (MQ=255)
       ttAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATagta                                                                                                                                        <  1:3301926/148‑1 (MQ=255)
         aGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCAAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCTGACATTATTGGCGTAAGAATTGATagtaag                                                                                                                                      >  2:1410735/1‑148 (MQ=255)
            cTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGtat                                                                                                                                   <  2:1803698/148‑1 (MQ=255)
            cTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGtat                                                                                                                                   <  1:1393192/148‑1 (MQ=255)
              ccACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTatat                                                                                                                                 >  1:453845/1‑148 (MQ=255)
                aCTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTAtatta                                                                                                                               <  1:3190152/148‑1 (MQ=255)
                 cTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTAtattat                                                                                                                              <  1:407957/148‑1 (MQ=255)
                 cTCGTGATGTGAGAGAGTTAATTAGTCTCATGTCTTATGGAATGATGTTAATATGTTTTCCAACAGGCATATTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTTGCGTAATAATTGATAGTAAGTAtattat                                                                                                                              <  2:1038567/148‑1 (MQ=255)
                  tCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATg                                                                                                                             <  2:583360/148‑1 (MQ=255)
                       aTGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGat                                                                                                                        <  2:3405656/148‑1 (MQ=255)
                           gAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGataata                                                                                                                    <  1:511120/148‑1 (MQ=255)
                           gAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGataata                                                                                                                    >  1:164913/1‑148 (MQ=255)
                            aGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGataataa                                                                                                                   >  2:1299596/1‑148 (MQ=255)
                             gCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGataataat                                                                                                                  <  1:2246617/148‑1 (MQ=255)
                                aGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATg                                                                                                               >  1:1697175/1‑148 (MQ=255)
                                 gTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATgg                                                                                                              >  2:393738/1‑148 (MQ=255)
                                 gTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATgg                                                                                                              <  1:3580659/148‑1 (MQ=255)
                                   tAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCt                                                                                                            <  2:2307568/148‑1 (MQ=255)
                                   tAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCt                                                                                                            <  1:252691/148‑1 (MQ=255)
                                     aTTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCttt                                                                                                          <  1:1917012/148‑1 (MQ=255)
                                      ttAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTa                                                                                                         <  1:1031461/148‑1 (MQ=255)
                                       tAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTAt                                                                                                        >  2:3052410/1‑148 (MQ=255)
                                          tctcATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTAtttt                                                                                                     >  1:993651/1‑148 (MQ=255)
                                             cATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTg                                                                                                  >  1:3446895/1‑148 (MQ=255)
                                             cATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTg                                                                                                  >  2:3605231/1‑148 (MQ=255)
                                              aTGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGt                                                                                                 >  1:1154724/1‑148 (MQ=255)
                                                 tCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCag                                                                                              >  2:6466/1‑148 (MQ=255)
                                                   ttATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCagga                                                                                            <  1:2477587/148‑1 (MQ=255)
                                                    tATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCaggag                                                                                           >  2:697602/1‑148 (MQ=255)
                                                        gAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGta                                                                                       >  1:2339906/1‑148 (MQ=255)
                                                        gAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGta                                                                                       >  2:1763947/1‑148 (MQ=255)
                                                         aaTGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGGATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGtat                                                                                      <  1:913852/148‑1 (MQ=255)
                                                         aaTGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGtat                                                                                      <  2:1123900/148‑1 (MQ=255)
                                                         aaTGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGtat                                                                                      <  2:2422737/148‑1 (MQ=255)
                                                              tgGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCa                                                                                 <  1:3551159/148‑1 (MQ=255)
                                                               ggTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCaa                                                                                <  2:1612042/148‑1 (MQ=255)
                                                                 tAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATg                                                                              >  2:3758675/1‑148 (MQ=255)
                                                                 tAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATg                                                                              >  1:966754/1‑148 (MQ=255)
                                                                   atatGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGt                                                                            <  1:2969828/148‑1 (MQ=255)
                                                                   atatGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGt                                                                            <  2:664238/148‑1 (MQ=255)
                                                                   atatGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGt                                                                            <  2:1726128/148‑1 (MQ=255)
                                                                           tCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATta                                                                    <  1:2857724/148‑1 (MQ=255)
                                                                              aaCAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATtaagt                                                                 <  1:524797/148‑1 (MQ=255)
                                                                               aCAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATtaagta                                                                <  2:329610/148‑1 (MQ=255)
                                                                                cAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATtaagtaa                                                               <  2:1662680/148‑1 (MQ=255)
                                                                                 aGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATtaagtaag                                                              <  1:1790088/148‑1 (MQ=255)
                                                                                  ggCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAgc                                                             >  2:3418579/1‑148 (MQ=255)
                                                                                   gCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAgcg                                                            >  2:191260/1‑148 (MQ=255)
                                                                                    cATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAgcgc                                                           <  1:3432784/148‑1 (MQ=255)
                                                                                       actATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCAtt                                                        >  1:1367450/3‑148 (MQ=255)
                                                                                       aGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCAtt                                                        <  1:437724/148‑1 (MQ=255)
                                                                                       aGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCAtt                                                        >  2:786191/1‑148 (MQ=255)
                                                                                           ttCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATaa                                                    <  2:3132130/148‑1 (MQ=255)
                                                                                           ttCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATaa                                                    >  1:3162717/1‑148 (MQ=255)
                                                                                            tCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATaaa                                                   >  1:3541124/1‑148 (MQ=255)
                                                                                            tCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATaaa                                                   >  2:1946076/1‑148 (MQ=255)
                                                                                            tCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATaaa                                                   <  2:1994402/148‑1 (MQ=255)
                                                                                            tCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATaaa                                                   <  1:1228838/148‑1 (MQ=255)
                                                                                               ttATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACaa                                                <  2:297181/148‑1 (MQ=255)
                                                                                                   tGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACaaataa                                            <  1:2048658/148‑1 (MQ=255)
                                                                                                    gCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACaaataaa                                           >  2:2362400/1‑148 (MQ=255)
                                                                                                     ctctAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACaaataaat                                          <  2:2890399/148‑1 (MQ=255)
                                                                                                        tAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTaa                                       <  2:2571489/148‑1 (MQ=255)
                                                                                                         aaTATTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAg                                      <  1:828020/148‑1 (MQ=255)
                                                                                                            atTTATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTa                                   <  2:2810150/148‑1 (MQ=255)
                                                                                                             tttATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTaa                                  <  2:418264/148‑1 (MQ=255)
                                                                                                             tttATAGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTaa                                  <  1:835095/148‑1 (MQ=255)
                                                                                                                ataGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTg                               <  1:1101114/148‑1 (MQ=255)
                                                                                                                ataGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTg                               <  2:2845935/148‑1 (MQ=255)
                                                                                                                 taGTGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGc                              <  2:2328620/148‑1 (MQ=255)
                                                                                                                  aGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGct                             >  1:2273284/1‑148 (MQ=255)
                                                                                                                  aGGGACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGct                             <  2:447194/148‑1 (MQ=255)
                                                                                                                   gggACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGctc                            <  2:757300/148‑1 (MQ=255)
                                                                                                                   gggACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGctc                            <  1:970404/148‑1 (MQ=255)
                                                                                                                    ggACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATAATGTCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCAGTATAAAAAAATAAATTAAGTTAATTGatct                           >  1:2649299/1‑148 (MQ=255)
                                                                                                                     gACTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCtt                          <  1:2023959/148‑1 (MQ=255)
                                                                                                                      aCTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCtta                         <  2:3192001/148‑1 (MQ=255)
                                                                                                                      aCTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCtta                         <  2:1154724/148‑1 (MQ=255)
                                                                                                                      aCTGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCtta                         <  2:3229783/148‑1 (MQ=255)
                                                                                                                        tGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCttatt                       <  2:3422767/148‑1 (MQ=255)
                                                                                                                        tGTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCttatt                       <  2:1205094/148‑1 (MQ=255)
                                                                                                                         gTATCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCttatta                      <  2:721768/148‑1 (MQ=255)
                                                                                                                            tCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCttattatta                   <  1:714509/148‑1 (MQ=255)
                                                                                                                            tCAGACATTATTGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCttattatta                   <  2:2043697/148‑1 (MQ=255)
                                                                                                                                  attattGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTa             <  1:3605231/148‑1 (MQ=255)
                                                                                                                                  attattGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTa             <  2:2446923/148‑1 (MQ=255)
                                                                                                                                   ttattGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTaa            <  2:3559627/148‑1 (MQ=255)
                                                                                                                                     attGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAAcc          <  1:3126620/148‑1 (MQ=255)
                                                                                                                                      ttGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCt         >  1:371401/1‑148 (MQ=255)
                                                                                                                                       tGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTg        <  2:886922/148‑1 (MQ=255)
                                                                                                                                       tGGCGTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTg        >  1:621773/1‑148 (MQ=255)
                                                                                                                                         gCGTAAGAATTGAGAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTggg      <  2:3446910/148‑1 (MQ=255)
                                                                                                                                           gTAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCa    <  2:512054/148‑1 (MQ=255)
                                                                                                                                            tAAGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCAt   >  1:1561568/1‑148 (MQ=255)
                                                                                                                                             aaGAATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCAtt  <  1:301578/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCATGTCTTATGGAATGATGGTAATATGTTTTCCAACAGGCATAGTATTCTTTATTGCTCTAATATTTATAGGGACTGTATCAGACATTATTGGCGTAAGGATTGATAGTAAGTATATTATGGCGATAATAATATGGCTTTATTTTCTGTCAGGAGGGTATATTCAATGGTTTGTATTAAGTAAGCGCATTATAAACAAATAAATTAAGTTAATTGCTCTTATTATTATATGTAACCTGGGCATT  >  CP009273/729312‑729600

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: