Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,809,935 G→T 5.4% A604A (GCG→GCT katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,809,9350GT 5.4% 835.2 / 13.8 219A604A (GCG→GCTkatEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base G (115/91);  new base T (11/2);  total (126/93)
Fisher's exact test for biased strand distribution p-value = 4.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 222 total aligned reads displayed.

CCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCT  >  CP009273/1809788‑1810083
                                                                                                                                                   |                                                                                                                                                    
 cAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAt                                                                                                                                                    >  1:1078687/1‑149 (MQ=255)
   ggCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAttt                                                                                                                                                  >  1:3479082/1‑149 (MQ=255)
   ggCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAttt                                                                                                                                                  <  2:1032056/149‑1 (MQ=255)
    gCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAtttt                                                                                                                                                 >  2:2748987/1‑149 (MQ=255)
       gTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACt                                                                                                                                              <  2:1742841/149‑1 (MQ=255)
          gCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTaa                                                                                                                                           <  1:1406065/149‑1 (MQ=255)
           cGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                                                                                                                          <  2:3176316/149‑1 (MQ=255)
             aaaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatga                                                                                                                                        <  2:3704016/149‑1 (MQ=255)
              aaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgat                                                                                                                                       >  1:2485193/1‑149 (MQ=255)
              aaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgat                                                                                                                                       >  2:3389930/1‑149 (MQ=255)
              aaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgat                                                                                                                                       <  1:3654742/149‑1 (MQ=255)
                 aTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGaa                                                                                                                                    >  1:344890/1‑149 (MQ=255)
                  tctcGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAg                                                                                                                                   <  2:1798370/149‑1 (MQ=255)
                    tcGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTg                                                                                                                                 <  2:917937/149‑1 (MQ=255)
                     cGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTga                                                                                                                                <  1:867128/149‑1 (MQ=255)
                       gTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgaga                                                                                                                              >  2:978922/1‑149 (MQ=255)
                            gAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGc                                                                                                                         <  2:61497/149‑1 (MQ=255)
                              actgactgACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAg                                                                                                                       >  1:772174/1‑149 (MQ=255)
                                tgactgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAc                                                                                                                     >  2:1868365/1‑149 (MQ=255)
                                 gactgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                                                                                                                    <  2:3041673/149‑1 (MQ=255)
                                 gactgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGAcc                                                                                                                    >  1:3025675/1‑149 (MQ=255)
                                  actgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCt                                                                                                                   <  1:3227482/149‑1 (MQ=255)
                                   ctgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCtt                                                                                                                  >  1:1527300/1‑149 (MQ=255)
                                    tgacGACCAGCTGAATATCACCCCACCTCCGGACGTCATCGGTCTGAAAAAGGATCCATCCTTTAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTc                                                                                                                 <  1:3140124/149‑1 (MQ=255)
                                     gacgacCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                                                                                                                >  2:3037488/1‑149 (MQ=255)
                                     gacgacCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                                                                                                                <  2:1482998/149‑1 (MQ=255)
                                         acCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGcc                                                                                                            <  2:1447321/149‑1 (MQ=255)
                                           cAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCAt                                                                                                          <  2:3591227/149‑1 (MQ=255)
                                             gCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtc                                                                                                        <  2:287893/149‑1 (MQ=255)
                                                 aaTATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCaa                                                                                                    >  2:859154/1‑149 (MQ=255)
                                                   tatCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAgg                                                                                                  <  1:2956946/149‑1 (MQ=255)
                                                   tatCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAgg                                                                                                  <  2:620477/149‑1 (MQ=255)
                                                       aCCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCt                                                                                              >  1:3554130/1‑149 (MQ=255)
                                                        ccccACCTCCGGACGTCAGCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg                                                                                             <  2:2028346/149‑1 (MQ=255)
                                                        ccccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg                                                                                             >  1:2648201/1‑149 (MQ=255)
                                                        ccccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg                                                                                             >  1:3856069/1‑149 (MQ=255)
                                                            aCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAgg                                                                                         <  2:2143492/149‑1 (MQ=255)
                                                                ccGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCaa                                                                                     >  2:2192713/1‑149 (MQ=255)
                                                                  ggACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAg                                                                                   >  1:3611647/1‑149 (MQ=255)
                                                                   gACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAgg                                                                                  <  1:2709635/149‑1 (MQ=255)
                                                                    aCGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGc                                                                                 <  2:3530396/149‑1 (MQ=255)
                                                                        cAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTc                                                                             >  2:3206617/1‑149 (MQ=255)
                                                                        cAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTc                                                                             >  2:2776967/1‑149 (MQ=255)
                                                                         aaCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCa                                                                            >  1:2794575/1‑149 (MQ=255)
                                                                         aaCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCa                                                                            <  2:474349/149‑1 (MQ=255)
                                                                         aaCGGTCTGAAAAAGGATCCAGCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCa                                                                            <  2:1895031/149‑1 (MQ=255)
                                                                          aCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCAt                                                                           >  2:2440676/1‑149 (MQ=255)
                                                                            ggTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGc                                                                         <  1:1232646/149‑1 (MQ=255)
                                                                            ggTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGc                                                                         >  1:326493/1‑149 (MQ=255)
                                                                             gTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGcc                                                                        >  2:2541600/1‑149 (MQ=255)
                                                                              tCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCa                                                                       <  1:3395826/149‑1 (MQ=255)
                                                                              tCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCa                                                                       >  1:1329381/1‑149 (MQ=255)
                                                                               cTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCaa                                                                      >  1:117356/1‑149 (MQ=255)
                                                                                 gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAAc                                                                    >  1:2141533/1‑149 (MQ=255)
                                                                                   aaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTg                                                                  >  1:1753439/1‑149 (MQ=255)
                                                                                      aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctc                                                               >  2:90635/1‑149 (MQ=255)
                                                                                       ggATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctct                                                              >  1:2047753/1‑149 (MQ=255)
                                                                                         atccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTAc                                                            >  2:3792515/1‑149 (MQ=255)
                                                                                             atccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCTATTTTACTTAATGATGACGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTccc                                                        >  2:3140776/1‑149 (MQ=255)
                                                                                             atccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTccc                                                        >  2:2209130/1‑149 (MQ=255)
                                                                                              tccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCg                                                       >  1:2228818/1‑149 (MQ=255)
                                                                                                   aaGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATgg                                                  >  2:1242066/1‑149 (MQ=255)
                                                                                                   aaGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATgg                                                  >  2:2934047/1‑149 (MQ=255)
                                                                                                     gTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGt                                                >  2:3616872/1‑149 (MQ=255)
                                                                                                      tttGTACGCCATTCCTGACGGTTATGTGAAATGTAGCGTGGTAGCTATTTTACTTAATGATGAAGTGAGATCGGCAGACATTCTGGCCATTCTCAAGGCGCTGAAGGCAAATGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGtt                                               <  2:504650/149‑2 (MQ=255)
                                                                                                      tttGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTg                                               >  2:2580291/1‑149 (MQ=255)
                                                                                                       gtgcaCGCCAGACCTGACGGAGATGTGAAAGGTCGCGTGGTAGCGCTTTTCCTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGa                                              <  2:1544078/145‑1 (MQ=255)
                                                                                                        tGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGaa                                             <  1:3954530/149‑1 (MQ=255)
                                                                                                         gTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAg                                            <  2:2485193/149‑1 (MQ=255)
                                                                                                         gTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAg                                            >  1:2827512/1‑149 (MQ=255)
                                                                                                          tACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGt                                           >  1:686465/1‑149 (MQ=255)
                                                                                                           aCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTg                                          <  2:2990836/149‑1 (MQ=255)
                                                                                                            cGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTg                                          >  1:3820093/1‑148 (MQ=255)
                                                                                                            cGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTg                                          <  2:3820093/148‑1 (MQ=255)
                                                                                                              ccATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACt                                       <  1:3149913/149‑1 (MQ=255)
                                                                                                               cATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTg                                      >  1:3054368/1‑149 (MQ=255)
                                                                                                                  tCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCgg                                   >  1:4048837/1‑149 (MQ=255)
                                                                                                                   ccTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGa                                  <  1:3571908/149‑1 (MQ=255)
                                                                                                                   ccTGACGGTGAAGTGAAAGGACGCGTGGTAGCTATATTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATACTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTTACTGAGGa                                  >  2:663167/1‑149 (MQ=255)
                                                                                                                    cTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGAt                                 >  1:1774232/1‑149 (MQ=255)
                                                                                                                      gACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGa                               <  1:3126072/149‑1 (MQ=255)
                                                                                                                      gACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATNCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGa                               <  1:2086050/149‑1 (MQ=255)
                                                                                                                         ggTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCGAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacgg                            >  2:2098563/1‑149 (MQ=255)
                                                                                                                          gTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggt                           >  2:2357863/1‑149 (MQ=255)
                                                                                                                           tGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggta                          <  2:448943/149‑1 (MQ=255)
                                                                                                                            gATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGacggtac                         >  1:1985270/1‑149 (MQ=255)
                                                                                                                                tgAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGgtg                     >  1:1582335/1‑149 (MQ=255)
                                                                                                                                tgAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGgtg                     >  2:114926/1‑149 (MQ=255)
                                                                                                                                  aaaGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGtt                   >  1:2592752/1‑149 (MQ=255)
                                                                                                                                   aaGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTg                  >  1:937893/1‑149 (MQ=255)
                                                                                                                                     ggTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCAGTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGcc                <  1:2765560/149‑1 (MQ=255)
                                                                                                                                       tCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCta              <  2:4012294/149‑1 (MQ=255)
                                                                                                                                         gcgTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCtata            >  1:1729725/1‑149 (MQ=255)
                                                                                                                                           gTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGc          <  2:1440865/149‑1 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTa     >  2:3650923/1‑149 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTa     >  1:386873/1‑149 (MQ=255)
                                                                                                                                                aGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGAGTTGCCTATAGCCGCTa     >  2:386720/1‑149 (MQ=255)
                                                                                                                                                 gCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCTAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTAc    >  1:942663/1‑149 (MQ=255)
                                                                                                                                                  cGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTAcc   <  2:3023687/149‑1 (MQ=255)
                                                                                                                                                   gATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCt  >  2:1268783/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTGACTGCGGATGACGGTACGGTGTTGCCTATAGCCGCTACCT  >  CP009273/1809788‑1810083

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: