Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,530,173 A→C 88.1% *399S (TAA→TCA) 
noncoding (62/244 nt)
ftsZ →
PP_mr19 →
Cell division protein FtsZ
Xoo8 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,530,1730AC88.1% 217.7 / 18.9 67*399S (TAA→TCA) 
noncoding (62/244 nt)
ftsZ
PP_mr19
Cell division protein FtsZ
Xoo8 ncRNA
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base C (23/36);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

AAGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACACCGGCATC  >  AE015451/1529950‑1530397
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                
aaGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       <  1:663389/235‑1 (MQ=255)
aaGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       <  1:318799/235‑1 (MQ=255)
         gACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGtata                                                                                                                                                                                                             >  1:626015/1‑236 (MQ=255)
              tACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGt                                                                                                                                                                                                        <  1:276668/236‑1 (MQ=255)
               ccccTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       <  1:228794/219‑1 (MQ=255)
                   ggCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTg                                                                                                                                                                                                   <  1:231429/235‑1 (MQ=255)
                        aCTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGt                                                                                                                                                                                                                >  1:320370/1‑218 (MQ=255)
                        aCTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTa                                                                                                                                                                                                >  2:425504/1‑234 (MQ=255)
                              cagcagGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGGAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGGGATTGGTATTCATCaa                                                                                                                                                                                        >  1:179424/1‑236 (MQ=255)
                               agcagGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCaaa                                                                                                                                                                                       >  1:86577/1‑236 (MQ=255)
                                  agGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       <  1:195523/201‑1 (MQ=255)
                                         cGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAAtt                                                                                                                                                                                                                          >  2:540896/1‑191 (MQ=255)
                                         cGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAAtt                                                                                                                                                                                                                          <  1:540898/191‑1 (MQ=255)
                                            aGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTg                                                                                                                                                                            >  2:679475/1‑234 (MQ=255)
                                                 ccAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTatcat                                                                                                                                                                     >  1:232095/1‑236 (MQ=255)
                                                       cTGCGCCGGTGCGCCAGGCGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTg                                                                                                                                                                                                   <  2:553644/200‑1 (MQ=255)
                                                           gcCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCtt                                                                                                                                                           >  1:509970/1‑236 (MQ=255)
                                                                  cgcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGAt                                                                                                                                                    >  1:195879/1‑236 (MQ=255)
                                                                   gcCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAgg                                                                                                                                                                                     <  2:692529/202‑1 (MQ=255)
                                                                    ccAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       <  2:98640/167‑1 (MQ=255)
                                                                             gccagcGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGc                                                                                                                                          <  1:574119/235‑1 (MQ=255)
                                                                                  cGGTCAACTACCGTGACCTGGAGCGCACGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAAttt                                                                                                                                     <  1:772611/235‑1 (MQ=255)
                                                                                         cTACCGGGCCCGGGCGCGCCCCACCGTCATGCGCCACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGCCTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGcc                                                                                                                                                             <  2:92944/204‑1 (MQ=255)
                                                                                                 aCCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                                     <  1:181575/236‑1 (MQ=255)
                                                                                                       aGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATcc                                                                                                                                                                   >  1:205953/1‑184 (MQ=255)
                                                                                                          gcCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAAttt                                                                                                                                     <  1:304930/211‑1 (MQ=255)
                                                                                                                ccGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaaca                                                                                                      <  1:679477/236‑1 (MQ=255)
                                                                                                                    cATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGc                                                                                                  <  1:69649/236‑1 (MQ=255)
                                                                                                                    cATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGATAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGc                                                                                                  <  1:262974/236‑1 (MQ=255)
                                                                                                                        cgcAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaacaac                                                                                                    <  1:138900/230‑1 (MQ=255)
                                                                                                                        cgcAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCAccc                                                                                              <  1:766829/236‑1 (MQ=255)
                                                                                                                           aaCCAGGCGCATGCAGGTGCTGCCGCAGACGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGc                                                                                                                                          <  1:715311/189‑1 (MQ=255)
                                                                                                                                ggCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTgg                                                                                                                                                                                 >  1:438117/1‑145 (MQ=255)
                                                                                                                                ggCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCaa                                                                                                                                        >  2:789923/1‑186 (MQ=255)
                                                                                                                                      tGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAg                                                                                                                                                               <  2:721518/157‑1 (MQ=255)
                                                                                                                                      tGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCGGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaaca                                                                                                      <  2:5312/214‑1 (MQ=255)
                                                                                                                                       gCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCAccc                                                                                              <  1:27131/221‑1 (MQ=255)
                                                                                                                                           gTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATg                                                                                                              <  1:329843/201‑1 (MQ=255)
                                                                                                                                                         gCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCa                                                                                                                      >  2:26315/1‑179 (MQ=255)
                                                                                                                                                           tAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCa                                                                                                                <  2:342763/183‑1 (MQ=255)
                                                                                                                                                                   aCCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCGGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCAccc                                                                                              <  2:199153/193‑1 (MQ=255)
                                                                                                                                                                      cTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGga                                                >  1:680671/1‑236 (MQ=255)
                                                                                                                                                                              gACCTGGACTACCTCTATATCCCGGCATTCCTGCGTCGTCACGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCaa                                                                                                                     <  1:784799/159‑1 (MQ=255)
                                                                                                                                                                                  tGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTAc                                      >  2:415601/1‑234 (MQ=255)
                                                                                                                                                                                   ggACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaaca                                                                                                      <  2:345444/169‑1 (MQ=255)
                                                                                                                                                                                   ggACTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAaaca                                                                                                      >  1:345446/1‑169 (MQ=255)
                                                                                                                                                                                    gACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGATCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCcac                                                                        >  2:112108/1‑198 (MQ=255)
                                                                                                                                                                                     aCTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGAc                                 <  1:208717/236‑1 (MQ=255)
                                                                                                                                                                                      cTACCTCGATATCCCGGCATTCCTGCGTTGTCAGGCTGATTAATGGAATTTAt                                                                                                                                                                                                                       <  2:254631/53‑1 (MQ=255)
                                                                                                                                                                                      cTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGc                                                                           <  2:86576/193‑1 (MQ=255)
                                                                                                                                                                                        aCCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTAc                                      >  1:635409/1‑228 (MQ=255)
                                                                                                                                                                                        aCCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTAc                                      <  2:635407/228‑1 (MQ=255)
                                                                                                                                                                                        aCCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGAccc                               >  1:751729/1‑235 (MQ=255)
                                                                                                                                                                                            cGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCAGCCGTg                                                                            >  2:120902/1‑186 (MQ=255)
                                                                                                                                                                                              atatCCCGGCATTCCTGCGTCGTCAGGCTGCTTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGc                         <  1:147307/235‑1 (MQ=255)
                                                                                                                                                                                                 tCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       <  2:747666/42‑1 (MQ=255)
                                                                                                                                                                                                 tCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTAt                                                                                                                                                                                                                       >  1:747668/1‑42 (MQ=255)
                                                                                                                                                                                                 tCCCGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTg                      >  1:529917/1‑235 (MQ=255)
                                                                                                                                                                                                 tCCCGGCATTACTGCTTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGAATCCATCAAGGGTTGTGTCTGTTATCATCCCCAGCCTTTCTTGATACCTGTTCGCAATTTGCGCTGagg                                                                                                                            >  2:782903/1‑131 (MQ=255)
                                                                                                                                                                                                   ccGGCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAAttt                                                                                                                                     <  2:666122/122‑1 (MQ=255)
                                                                                                                                                                                                      gCATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCAcc                  >  2:414067/1‑234 (MQ=255)
                                                                                                                                                                                                       cATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGt               <  1:415603/236‑1 (MQ=255)
                                                                                                                                                                                                       cATTCCTGCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGt               <  1:425506/236‑1 (MQ=255)
                                                                                                                                                                                                              gCGTCGTCAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGac           >  2:692907/1‑233 (MQ=255)
                                                                                                                                                                                                                    tcAGGCTGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACACCGGCATc  <  1:480184/236‑1 (MQ=255)
                                                                                                                                                                                                                       ggCTGATTCATGGAATTTATCAAGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCg                                                                              >  2:577422/1‑157 (MQ=255)
                                                                                                                                                                                                                          tGATTCATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTAGCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGGCGGTCt                                                            <  2:239877/172‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                
AAGGTGGTCGACAATACCCTGCAGACTGCCCAGCAGGCGTACGAAGCGTCCAATCCTGCGCCGGTGCGCCAGGAGCAGCCAGCGGTCAACTACCGTGACCTGGAGCGCCCGACCGTCATGCGCAACCAGGCGCATGCAGGTGCTGCCGCAGCCGCTAAACTGAACCCTCAGGATGACCTGGACTACCTCGATATCCCGGCATTCCTGCGTCGTCAGGCTGATTAATGGAATTTATCAGGGGTATAAAGGTGATTGGTATTCATCAAAGGTTGGGTCTGTTATCATCCCCAGCCTTTGTTGATACCTGTTCGCAATTTGCGCTGAAGCGGCCAATGCCATGATTAAACAACGCACCCTGAAGAATACCATCCGTGCCACAGGTGTCGGTCTGCACTCCGGGGAGAAGGTTTACCTGACCCTCAAGCCTGCACCTGTTGACACCGGCATC  >  AE015451/1529950‑1530397

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: