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,0220TG98.7% 1256.8 / ‑0.1 313E188D (GAA→GACyijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (174/135);  total (176/137)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 322 total aligned reads displayed.

GATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTT  >  CP009273/4126875‑4127168
                                                                                                                                                   |                                                                                                                                                  
 attgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACAtcgtct                                                                                                                                                 <  1:188922/150‑1 (MQ=255)
  ttgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTTTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATcgtctt                                                                                                                                                <  1:2016513/150‑1 (MQ=255)
   tgtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCgtcttc                                                                                                                                               <  1:3536596/150‑1 (MQ=255)
    gtgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTcc                                                                                                                                              >  2:838642/1‑150 (MQ=255)
     tgctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCt                                                                                                                                             <  2:3184843/150‑1 (MQ=255)
       ctgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTgc                                                                                                                                           >  2:3560227/1‑150 (MQ=255)
        tgGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTgcg                                                                                                                                          >  1:564147/1‑150 (MQ=255)
            aCTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGcc                                                                                                                                      <  1:1105094/150‑1 (MQ=255)
              tGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTg                                                                                                                                    <  1:3044347/150‑1 (MQ=255)
               ggCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTgc                                                                                                                                   <  1:4041207/150‑1 (MQ=255)
                gcgcTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTgcg                                                                                                                                  <  2:1797923/150‑1 (MQ=255)
                     ggTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAAt                                                                                                                             <  1:3057593/150‑1 (MQ=255)
                        gAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGtt                                                                                                                          >  1:3145379/1‑150 (MQ=255)
                         aGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTa                                                                                                                         >  2:2972551/1‑150 (MQ=255)
                           ccGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACTATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATg                                                                                                                       >  1:2155797/1‑150 (MQ=255)
                            cGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGt                                                                                                                      >  1:1819892/1‑150 (MQ=255)
                               aaGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATc                                                                                                                   >  1:2530370/1‑150 (MQ=255)
                                 gAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAAAAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAc                                                                                                                 >  2:4385404/1‑150 (MQ=255)
                                  aaaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCAcc                                                                                                                <  2:4180619/150‑1 (MQ=255)
                                   aaaaaGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCa                                                                                                               >  2:1283048/1‑150 (MQ=255)
                                       agaagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcgg                                                                                                           <  2:2393176/150‑1 (MQ=255)
                                         aagCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAgcggcg                                                                                                         <  2:44458/150‑1 (MQ=255)
                                           gCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTa                                                                                                       <  1:2299583/150‑1 (MQ=255)
                                            ccGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTAc                                                                                                      >  2:300520/1‑150 (MQ=255)
                                             cGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTAcc                                                                                                     <  1:638041/150‑1 (MQ=255)
                                             cGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTAcc                                                                                                     >  1:4523241/1‑150 (MQ=255)
                                              gTCAGTTTAAAAAGGATCTGAATTTTTTCTTCATGTAGGGCGATCTCTTATTTAACAAAATAACTATAATGCCCCACCATCCGCCAGATAAACAGCACATCGTCTTCCTGCGCGACAGCGCCAATTTTAAGAATCACAAGCGGCGTACCg                                                                                                    >  2:1380969/1‑150 (MQ=255)
                                                 aGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcg                                                                                                 >  1:970382/1‑150 (MQ=255)
                                                 aGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcg                                                                                                 >  2:793435/1‑150 (MQ=255)
                                                  gtgtAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTcgc                                                                                                >  2:391063/1‑150 (MQ=255)
                                                     tAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCggg                                                                                             >  1:759026/1‑150 (MQ=255)
                                                      aaaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                            <  1:3587437/150‑1 (MQ=255)
                                                      aaaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                            <  2:1617341/150‑1 (MQ=255)
                                                      aaaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                            >  1:2584419/1‑150 (MQ=255)
                                                      aaaaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGc                                                                                            <  2:3277983/150‑1 (MQ=255)
                                                        aaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGa                                                                                          >  2:1925993/1‑150 (MQ=255)
                                                        aaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGa                                                                                          <  1:2093544/150‑1 (MQ=255)
                                                        aaaGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGa                                                                                          >  2:1073621/1‑150 (MQ=255)
                                                         aaGGATCTTAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTTCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGaa                                                                                         <  1:1701016/150‑1 (MQ=255)
                                                          aGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAg                                                                                        >  2:3385323/1‑150 (MQ=255)
                                                            gATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGcc                                                                                      >  1:1135409/1‑150 (MQ=255)
                                                             aTCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCa                                                                                     <  1:1433698/150‑1 (MQ=255)
                                                                 gAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTg                                                                                 >  2:1165281/1‑150 (MQ=255)
                                                                  aaTTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGa                                                                                >  1:2901643/1‑150 (MQ=255)
                                                                       tttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacca                                                                           >  2:1884973/1‑150 (MQ=255)
                                                                       tttCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAacca                                                                           >  1:1294746/1‑150 (MQ=255)
                                                                          cttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAccc                                                                        <  2:10353/150‑1 (MQ=255)
                                                                          cttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAccc                                                                        >  1:2730879/1‑150 (MQ=255)
                                                                          cttcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCAccc                                                                        >  1:3622356/1‑150 (MQ=255)
                                                                            tcATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCa                                                                      >  2:2190700/1‑150 (MQ=255)
                                                                                gtgGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAatat                                                                  <  1:4604709/150‑1 (MQ=255)
                                                                                  ggggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATAtgt                                                                <  1:872404/150‑1 (MQ=255)
                                                                                    ggCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGc                                                              >  1:4157047/1‑150 (MQ=255)
                                                                                       gATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                           >  2:2909990/1‑150 (MQ=255)
                                                                                       gATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAg                                                           >  2:3870408/1‑150 (MQ=255)
                                                                                         tctcTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGcc                                                         >  2:1446395/1‑150 (MQ=255)
                                                                                              tATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTa                                                    <  2:677644/150‑1 (MQ=255)
                                                                                               aTTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTAt                                                   >  1:4261191/1‑150 (MQ=255)
                                                                                                  tAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCa                                                >  2:38220/1‑150 (MQ=255)
                                                                                                  tAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCa                                                >  2:1309745/1‑150 (MQ=255)
                                                                                                    aCAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGt                                              <  1:2419721/150‑1 (MQ=255)
                                                                                                     cAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTc                                             <  1:1568503/150‑1 (MQ=255)
                                                                                                     cAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTc                                             >  2:4636130/1‑150 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACTAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            <  2:1969786/150‑1 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            >  2:3096740/1‑150 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            >  2:2734186/1‑150 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            <  1:3314804/150‑1 (MQ=255)
                                                                                                      aaaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCg                                            >  2:2853629/1‑150 (MQ=255)
                                                                                                       aaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGc                                           >  2:3674132/1‑150 (MQ=255)
                                                                                                       aaaTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGc                                           >  1:2647564/1‑150 (MQ=255)
                                                                                                         aTAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCa                                         <  1:3448857/150‑1 (MQ=255)
                                                                                                              gATAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGGTAGGTCTCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATCTGTGCCAGCCCGTTATCCAGTCGCCAGgaga                                    <  2:533603/150‑1 (MQ=255)
                                                                                                               aTAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGAc                                   >  2:818743/1‑150 (MQ=255)
                                                                                                                tAATGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACa                                  >  1:4373205/1‑150 (MQ=255)
                                                                                                                 aaTGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACaa                                 >  2:1175790/1‑150 (MQ=255)
                                                                                                                   tGCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAata                               <  2:1151210/150‑1 (MQ=255)
                                                                                                                    gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAatat                              <  2:4052301/150‑1 (MQ=255)
                                                                                                                    gCCCCACCATCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAatat                              >  2:556565/1‑150 (MQ=255)
                                                                                                                       ccaccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGc                           >  2:508699/1‑150 (MQ=255)
                                                                                                                        caccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                          >  1:1891930/1‑150 (MQ=255)
                                                                                                                        caccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGcc                          <  2:1560902/150‑1 (MQ=255)
                                                                                                                         accaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                         >  2:2163826/1‑150 (MQ=255)
                                                                                                                         accaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGccc                         >  2:2769031/1‑150 (MQ=255)
                                                                                                                          ccaTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCg                        <  1:1356707/150‑1 (MQ=255)
                                                                                                                           caTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGc                       >  2:3469720/1‑150 (MQ=255)
                                                                                                                            aTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCt                      <  2:3145379/150‑1 (MQ=255)
                                                                                                                            aTCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCt                      >  2:3979510/1‑150 (MQ=255)
                                                                                                                             tCCGCCAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCTTTATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCtt                     <  2:2220630/150‑1 (MQ=255)
                                                                                                                                 ccAGTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATa                 >  2:1229415/1‑150 (MQ=255)
                                                                                                                                    gTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTc              <  2:2530370/150‑1 (MQ=255)
                                                                                                                                    gTTAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTc              <  1:2210628/150‑1 (MQ=255)
                                                                                                                                     ttAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCa             >  2:2607921/1‑150 (MQ=255)
                                                                                                                                     ttAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCa             >  1:3265090/1‑150 (MQ=255)
                                                                                                                                      tAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCAc            <  2:3231528/150‑1 (MQ=255)
                                                                                                                                      tAAACAGCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCAc            <  2:256903/150‑1 (MQ=255)
                                                                                                                                         aCAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTgg         >  2:4243475/1‑150 (MQ=255)
                                                                                                                                            gCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGATACAATATCGCCCGCTTGATAGTCACTGGGGt      >  2:3327538/1‑150 (MQ=255)
                                                                                                                                            gCACATCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGt      <  1:296249/150‑1 (MQ=255)
                                                                                                                                               caTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCt   <  1:3780013/150‑1 (MQ=255)
                                                                                                                                                aTCGTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCtt  <  1:4193645/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GATTGTGCTGGCACTGGCGCTGGTGAGCCGTAAGAAAAAAGAAGCCGTCAGTGTAAAAAGGATCTGAATTTTTTCTTCATGTGGGGCGATCTCTTATTTAACAAAATAACGATAATGCCCCACCATCCGCCAGTTAAACAGCACATCTTCTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTT  >  CP009273/4126875‑4127168

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: