Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,127,022 T→G 100% E188D (GAA→GAC yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,127,0220TG99.6% 921.3 / ‑5.4 228E188D (GAA→GACyijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (127/99);  total (128/100)
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: E-value score below prediction cutoff.
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 231 total aligned reads displayed.

TTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTC  >  CP009273/4126877‑4127166
                                                                                                                                                 |                                                                                                                                                
            tGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGcc                                                                                                                                    <  2:692108/148‑1 (MQ=255)
                   ggTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCa                                                                                                                             >  1:1118044/1‑148 (MQ=255)
                       aGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGt                                                                                                                         <  2:2285404/148‑1 (MQ=255)
                        gCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                        >  2:1610998/1‑148 (MQ=255)
                        gCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                        <  2:1795871/148‑1 (MQ=255)
                        gCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                        <  1:1512962/148‑1 (MQ=255)
                          cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTAt                                                                                                                      <  1:46079/148‑1 (MQ=255)
                             aaGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTa                                                                                                                   >  2:1591904/1‑148 (MQ=255)
                             aaGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTa                                                                                                                   <  2:2398327/148‑1 (MQ=255)
                             aaGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTa                                                                                                                   >  2:372173/1‑148 (MQ=255)
                              aGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTAt                                                                                                                  <  2:1152356/148‑1 (MQ=255)
                                 aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAc                                                                                                               >  1:2294276/1‑148 (MQ=255)
                                 aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAc                                                                                                               >  1:501821/1‑148 (MQ=255)
                                  aaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAcc                                                                                                              >  1:1920358/1‑148 (MQ=255)
                                   aaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCa                                                                                                             >  1:140300/1‑148 (MQ=255)
                                   aaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACCATAATCCCCCACAATCCGCCAGTTAAACAGAACATCGTCTTCCTGCTCCCCACCGACAATCTTATGTCTCACCa                                                                                                             >  2:1533977/1‑148 (MQ=255)
                                    aagaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAg                                                                                                            <  1:2696733/148‑1 (MQ=255)
                                    aagaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAg                                                                                                            >  2:3091491/1‑148 (MQ=255)
                                     agaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgc                                                                                                           <  1:2587645/148‑1 (MQ=255)
                                      gaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                          >  1:2595662/1‑148 (MQ=255)
                                      gaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                          >  1:1457080/1‑148 (MQ=255)
                                      gaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcg                                                                                                          >  1:1234128/1‑148 (MQ=255)
                                         gCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcggcg                                                                                                       >  2:1730217/1‑148 (MQ=255)
                                          ccGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGt                                                                                                      >  2:1515404/1‑148 (MQ=255)
                                            gTCAGTGTAAACAGTATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTAc                                                                                                    <  2:2502003/148‑1 (MQ=255)
                                            gTCAGTGTAAAAATGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTAc                                                                                                    <  1:677970/148‑1 (MQ=255)
                                             tCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCACTGTTATGTATCACCAGCGGCGTAcc                                                                                                   >  2:535852/1‑148 (MQ=255)
                                              cAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCg                                                                                                  >  2:3236459/1‑148 (MQ=255)
                                                 tgtAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcg                                                                                               <  1:2429223/148‑1 (MQ=255)
                                                   tAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcgcg                                                                                             >  2:2166376/1‑148 (MQ=255)
                                                      aaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                          <  1:3086302/148‑1 (MQ=255)
                                                       aaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCg                                                                                         >  1:3011832/1‑148 (MQ=255)
                                                       aaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCg                                                                                         >  1:1287316/1‑148 (MQ=255)
                                                          gATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAg                                                                                      >  1:2835373/1‑148 (MQ=255)
                                                          gATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAg                                                                                      <  1:2366907/148‑1 (MQ=255)
                                                          gATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAg                                                                                      <  2:1447837/148‑1 (MQ=255)
                                                          gATCTGAATTTTTTCTTCATGTGGGGCGATCACATATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAg                                                                                      >  1:1704044/1‑148 (MQ=255)
                                                            tCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGcc                                                                                    <  2:1029742/148‑1 (MQ=255)
                                                             cTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCa                                                                                   >  2:2320461/1‑148 (MQ=255)
                                                              tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                  >  2:3013399/1‑148 (MQ=255)
                                                              tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                  <  1:862566/148‑1 (MQ=255)
                                                              tGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCAt                                                                                  <  2:2985028/148‑1 (MQ=255)
                                                               gAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATc                                                                                 >  1:2133573/1‑148 (MQ=255)
                                                                aaTTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCt                                                                                >  2:2411719/1‑148 (MQ=255)
                                                                       tcttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacca                                                                         >  2:85031/1‑148 (MQ=255)
                                                                       tcttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacca                                                                         >  2:246533/1‑148 (MQ=255)
                                                                        cttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAaccac                                                                        >  2:414094/1‑148 (MQ=255)
                                                                          tcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAccc                                                                      <  1:1008740/148‑1 (MQ=255)
                                                                          tcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAccc                                                                      <  1:1008650/148‑1 (MQ=255)
                                                                           cATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAcccc                                                                     <  2:1118044/148‑1 (MQ=255)
                                                                             tgtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCaa                                                                   >  1:1197469/1‑148 (MQ=255)
                                                                              gtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAat                                                                  <  1:790208/148‑1 (MQ=255)
                                                                              gtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAat                                                                  <  1:544490/148‑1 (MQ=255)
                                                                               tgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAata                                                                 <  2:501821/148‑1 (MQ=255)
                                                                               tgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAata                                                                 <  2:240098/148‑1 (MQ=255)
                                                                               tgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAata                                                                 <  2:2133573/148‑1 (MQ=255)
                                                                                 gggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtg                                                               >  1:2178906/1‑148 (MQ=255)
                                                                                 gggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtg                                                               <  2:1988304/148‑1 (MQ=255)
                                                                                 gggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtg                                                               >  1:383125/1‑148 (MQ=255)
                                                                                 gggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtg                                                               <  1:367324/148‑1 (MQ=255)
                                                                                   gCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtgtg                                                             >  1:719088/1‑148 (MQ=255)
                                                                                    cGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGc                                                            >  2:351172/1‑148 (MQ=255)
                                                                                    cGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGc                                                            <  1:2518492/148‑1 (MQ=255)
                                                                                      aTCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGACa                                                          >  1:2637120/1‑148 (MQ=255)
                                                                                        ctctTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGc                                                        >  1:3244699/1‑148 (MQ=255)
                                                                                            tATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGt                                                    >  2:3275895/1‑148 (MQ=255)
                                                                                                    aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTGCCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGt                                            >  2:1324199/1‑148 (MQ=255)
                                                                                                      aaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                          >  2:2112657/1‑148 (MQ=255)
                                                                                                         aaCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCa                                       <  1:2995642/148‑1 (MQ=255)
                                                                                                         aaCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCa                                       <  2:1234128/148‑1 (MQ=255)
                                                                                                         aaCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCa                                       >  2:2457338/1‑148 (MQ=255)
                                                                                                          aCGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAg                                      >  2:106775/1‑148 (MQ=255)
                                                                                                            gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGga                                    <  1:2329289/148‑1 (MQ=255)
                                                                                                            gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGga                                    >  2:2631616/1‑148 (MQ=255)
                                                                                                            gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGACAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGAGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGga                                    >  2:330314/1‑148 (MQ=255)
                                                                                                             aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGgag                                   >  2:1098166/1‑148 (MQ=255)
                                                                                                               aaTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGAc                                 >  1:1861639/1‑148 (MQ=255)
                                                                                                                aTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                                <  1:423549/148‑1 (MQ=255)
                                                                                                                 tGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                               <  2:59717/148‑1 (MQ=255)
                                                                                                                  gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAat                              <  1:173660/148‑1 (MQ=255)
                                                                                                                  gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAat                              <  2:2813213/148‑1 (MQ=255)
                                                                                                                         caTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                       >  2:2332796/1‑148 (MQ=255)
                                                                                                                          aTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                      >  1:2916867/1‑148 (MQ=255)
                                                                                                                           tCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                     <  2:311302/148‑1 (MQ=255)
                                                                                                                           tCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                     <  2:654984/148‑1 (MQ=255)
                                                                                                                           tCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                     >  2:2334675/1‑148 (MQ=255)
                                                                                                                           tCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                     <  2:2443564/148‑1 (MQ=255)
                                                                                                                               ccAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                 <  1:1965298/148‑1 (MQ=255)
                                                                                                                               ccAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGa                 >  1:2048494/1‑148 (MQ=255)
                                                                                                                                cAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGAt                <  1:414094/148‑1 (MQ=255)
                                                                                                                                 aGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa               <  1:3236459/148‑1 (MQ=255)
                                                                                                                                   ttAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGt             <  1:1098166/148‑1 (MQ=255)
                                                                                                                                     aaaCAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCa           <  2:1920358/148‑1 (MQ=255)
                                                                                                                                     aaaCAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCa           >  1:1578122/1‑148 (MQ=255)
                                                                                                                                         aGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgg       <  1:3527155/148‑1 (MQ=255)
                                                                                                                                         aGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgg       >  1:790763/1‑148 (MQ=255)
                                                                                                                                          gCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTggg      >  1:2404972/1‑148 (MQ=255)
                                                                                                                                           cacaTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg     <  1:3190705/148‑1 (MQ=255)
                                                                                                                                           cacaTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgggg     >  1:140725/1‑148 (MQ=255)
                                                                                                                                             caTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGtt   >  1:245257/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTC  >  CP009273/4126877‑4127166

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: