Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 54,650 C→T 8.2% S18N (AGT→AAT)  surA ← peptidyl‑prolyl cis‑trans isomerase (PPIase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927354,6500CT 8.2% 722.4 / 37.6 233S18N (AGT→AAT) surApeptidyl‑prolyl cis‑trans isomerase (PPIase)
Reads supporting (aligned to +/- strand):  ref base C (91/122);  new base T (7/12);  total (99/134)
Fisher's exact test for biased strand distribution p-value = 8.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 236 total aligned reads displayed.

GGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACGCAGCATCTC  >  CP009273/54501‑54799
                                                                                                                                                     |                                                                                                                                                     
ggCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAAc                                                                                                                                                       >  1:2643012/1‑150 (MQ=255)
  cgcAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTg                                                                                                                                                     <  2:1234417/150‑1 (MQ=255)
  cgcAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTATCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTg                                                                                                                                                     >  2:2339057/1‑150 (MQ=255)
     aGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTa                                                                                                                                                  >  1:789463/1‑150 (MQ=255)
      gcgtcgcgtcATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTAt                                                                                                                                                 <  1:1536599/150‑1 (MQ=255)
       cgtcgcgtcATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTAtt                                                                                                                                                >  1:234092/1‑150 (MQ=255)
       cgtcgcgtcATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTAtt                                                                                                                                                <  2:4529952/150‑1 (MQ=255)
        gtcgcgtcATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTc                                                                                                                                               >  1:2693993/1‑150 (MQ=255)
        gtcgcgtcATCAGGAAGTTGCTGCCGTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTc                                                                                                                                               <  2:3956083/150‑1 (MQ=255)
             gtcATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGAt                                                                                                                                          <  1:3635490/150‑1 (MQ=255)
             gtcATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGAt                                                                                                                                          <  1:1574404/150‑1 (MQ=255)
              tcatcaGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATc                                                                                                                                         <  1:3749610/150‑1 (MQ=255)
                atcaGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCAt                                                                                                                                       >  2:1575667/1‑150 (MQ=255)
                atcaGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCATTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCAt                                                                                                                                       >  2:305835/1‑150 (MQ=255)
                     gAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGa                                                                                                                                  >  1:2736461/1‑150 (MQ=255)
                      aaGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGAt                                                                                                                                 >  1:2787684/1‑150 (MQ=255)
                      aaGGTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGAt                                                                                                                                 <  2:57688/150‑1 (MQ=255)
                         ttGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATAcc                                                                                                                              <  2:3169008/150‑1 (MQ=255)
                          tgctgcCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTACTCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCg                                                                                                                             <  2:607668/150‑1 (MQ=255)
                          tgctgcCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCg                                                                                                                             >  1:1962032/1‑150 (MQ=255)
                           gctgcCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCga                                                                                                                            >  1:1217303/1‑150 (MQ=255)
                             tgccttgcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCgaga                                                                                                                          <  2:4523043/150‑1 (MQ=255)
                             tgccttgcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCgaga                                                                                                                          >  2:348120/1‑150 (MQ=255)
                               ccttgcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAag                                                                                                                        <  1:4333445/150‑1 (MQ=255)
                               ccttgcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAag                                                                                                                        <  1:2046197/150‑1 (MQ=255)
                                cttgcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAagc                                                                                                                       >  1:604446/1‑150 (MQ=255)
                                  tgcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAagcag                                                                                                                     <  2:2258151/150‑1 (MQ=255)
                                   gcctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAagcagc                                                                                                                    <  2:2039723/150‑1 (MQ=255)
                                    cctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGACGCAGCg                                                                                                                   >  2:4662725/1‑150 (MQ=255)
                                    cctGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCg                                                                                                                   <  1:3867170/150‑1 (MQ=255)
                                      tGAGCAGCGTTCAGTTTTCCCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGtt                                                                                                                 <  1:890315/150‑1 (MQ=255)
                                      tGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGtt                                                                                                                 <  2:3734672/150‑1 (MQ=255)
                                          cagcGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTcc                                                                                                             >  2:629517/1‑150 (MQ=255)
                                          cagcGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTcc                                                                                                             >  1:588038/1‑150 (MQ=255)
                                          cagcGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTcc                                                                                                             <  2:787644/150‑1 (MQ=255)
                                           agcGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCa                                                                                                            >  2:2005568/1‑150 (MQ=255)
                                            gcGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAg                                                                                                           <  1:1608006/150‑1 (MQ=255)
                                              gTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGtt                                                                                                         <  1:2684547/150‑1 (MQ=255)
                                                tCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGttct                                                                                                       <  1:4375622/150‑1 (MQ=255)
                                                 cAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGttctt                                                                                                      <  1:3154730/150‑1 (MQ=255)
                                                 cAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGttctt                                                                                                      >  2:4873875/1‑150 (MQ=255)
                                                  aGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGttcttc                                                                                                     <  2:3602563/150‑1 (MQ=255)
                                                          cGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACtttt                                                                                             <  1:4106106/150‑1 (MQ=255)
                                                               gCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCAt                                                                                        <  1:4989527/150‑1 (MQ=255)
                                                                 aTTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCAttt                                                                                      >  2:1239405/1‑150 (MQ=255)
                                                                 aTTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCAttt                                                                                      >  2:1506128/1‑150 (MQ=255)
                                                                  ttAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTc                                                                                     <  1:1765898/150‑1 (MQ=255)
                                                                  ttAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTc                                                                                     <  2:3567297/150‑1 (MQ=255)
                                                                  ttAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTc                                                                                     <  2:1107155/150‑1 (MQ=255)
                                                                      tCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAAtt                                                                                 >  1:4975834/1‑150 (MQ=255)
                                                                      tCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAAtt                                                                                 >  1:2604281/1‑150 (MQ=255)
                                                                      tCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAAtt                                                                                 <  2:1817114/150‑1 (MQ=255)
                                                                       ccATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATTGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTa                                                                                >  1:337070/1‑150 (MQ=255)
                                                                          tCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAAcc                                                                             <  1:1877685/150‑1 (MQ=255)
                                                                            aaCGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGc                                                                           >  1:3406485/1‑150 (MQ=255)
                                                                            aaCGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTATCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGc                                                                           >  2:4471394/1‑150 (MQ=255)
                                                                                   cTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCgg                                                                    >  2:1040611/1‑150 (MQ=255)
                                                                                    tttCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGa                                                                   <  1:3243534/150‑1 (MQ=255)
                                                                                        cAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTAc                                                               >  1:3978845/1‑150 (MQ=255)
                                                                                          gCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGt                                                             >  1:868917/1‑150 (MQ=255)
                                                                                             cgacgCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGt                                                          >  1:4686954/1‑150 (MQ=255)
                                                                                             cgacgCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGt                                                          >  2:2408390/1‑150 (MQ=255)
                                                                                              gacgCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTa                                                         >  2:1746885/1‑150 (MQ=255)
                                                                                                    gTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATcaa                                                   >  2:4329579/1‑150 (MQ=255)
                                                                                                      tATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGAGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATcaaca                                                 >  2:355224/1‑150 (MQ=255)
                                                                                                       aTTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATcaacaa                                                <  2:680702/150‑1 (MQ=255)
                                                                                                          gacgacGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATc                                             <  1:688375/150‑1 (MQ=255)
                                                                                                             gacgGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATcaca                                          <  2:4206685/150‑1 (MQ=255)
                                                                                                               cgGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACaaa                                        <  1:4288459/150‑1 (MQ=255)
                                                                                                               cgGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACaaa                                        <  1:1121153/150‑1 (MQ=255)
                                                                                                                ggCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAg                                       >  1:4067417/1‑150 (MQ=255)
                                                                                                                 gCTGCGACTTTATCGACTACCTGGGGGGCAGCGAATCTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAgt                                      <  2:5009470/150‑1 (MQ=255)
                                                                                                                  cTGCGACTTTATCGACTACCTGGGGGGCAGCTAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAgtg                                     <  1:3049327/150‑1 (MQ=255)
                                                                                                                  cTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAgtg                                     <  1:5100288/150‑1 (MQ=255)
                                                                                                                  cTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAgtg                                     <  1:4508276/150‑1 (MQ=255)
                                                                                                                     cGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGttt                                  <  2:3751913/150‑1 (MQ=255)
                                                                                                                     cGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGttt                                  >  1:2351850/1‑150 (MQ=255)
                                                                                                                       aCTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTg                                <  2:3199004/150‑1 (MQ=255)
                                                                                                                        cTTTATCGACTACCTGGGGGGCAGCTAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTACTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGa                               >  2:1791928/1‑150 (MQ=255)
                                                                                                                         tttATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGAt                              <  2:4686954/150‑1 (MQ=255)
                                                                                                                         tttATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGAt                              <  2:735219/150‑1 (MQ=255)
                                                                                                                           tATCGACTACCTGGGGGGCAGCGAAATTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATAc                            <  2:566404/150‑1 (MQ=255)
                                                                                                                                   aCCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAAt                    <  1:2474495/150‑1 (MQ=255)
                                                                                                                                   aCCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAAt                    <  1:1216289/150‑1 (MQ=255)
                                                                                                                                    ccTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATg                   <  1:1021981/150‑1 (MQ=255)
                                                                                                                                     cTGGGGGGCAGCGTAACTGGTATTCGCGATCATGGCGATACGGAGACGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAAACGCACTGCGGATTACGTGGTAAATCAACAAATCTGAAAGTGTTTTGATACGGCAGAATGt                  <  2:4459866/150‑1 (MQ=255)
                                                                                                                                     cTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGt                  <  2:2867211/150‑1 (MQ=255)
                                                                                                                                     cTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGt                  <  2:587999/150‑1 (MQ=255)
                                                                                                                                      tGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGtt                 <  2:3322557/150‑1 (MQ=255)
                                                                                                                                       ggggggCAGCGAAATTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTc                >  2:1644435/1‑150 (MQ=255)
                                                                                                                                         ggggCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGa              <  2:2181290/150‑1 (MQ=255)
                                                                                                                                          gggCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGaa             <  2:2229943/150‑1 (MQ=255)
                                                                                                                                              aGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACgca         >  1:1033420/1‑150 (MQ=255)
                                                                                                                                              aGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACgca         >  2:131422/1‑150 (MQ=255)
                                                                                                                                               gCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACgcag        <  2:3406485/150‑1 (MQ=255)
                                                                                                                                                 gAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTAGAACgcagca      >  2:2611942/1‑150 (MQ=255)
                                                                                                                                                  aaaCTGGTATTCGCGATCATGGCGAGACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACGCAGCAt     <  2:2957864/150‑1 (MQ=255)
                                                                                                                                                   aaCTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACGCAGCAtc    >  1:2145726/1‑150 (MQ=255)
                                                                                                                                                   aaCTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACGCAGCAtc    <  2:4908185/150‑1 (MQ=255)
                                                                                                                                                     cTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACGCAGCAtctc  >  1:1352459/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCTTGCCTGAGCAGCGTTCAGTTTTACCGACTGCATTAATCCATCAACGTCGCTTTCCAGCACGACGCCGTTATTGACGACGGCTGCGACTTTATCGACTACCTGGGGGGCAGCGAAACTGGTATTCGCGATCATGGCGATACCGAGAAGCAGCGTTTTCCAGTTCTTCATACTTTTTCCATTTCAATTAACCGCACTGCGGATTACGTGGTAAATCAACAAATCACAAAGTGTTTTGATACGGCAGAATGTTCGAACGCAGCATCTC  >  CP009273/54501‑54799

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: