Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 5,674,750:1 +G 100% intergenic (+112/+28) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding protein, acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154515,674,7501.G97.1% 245.9 / 0.4 69intergenic (+112/+28)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding protein, acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (31/36);  total (31/38)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGT  >  AE015451/5674523‑5674975
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
cAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggccccct                                                                                                                                                                                                                           >  1:713149/1‑229 (MQ=255)
            gCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCGGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                                                                                                          >  1:265707/1‑226 (MQ=255)
                       aaGCTGGACCGTTGCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTa                                                                                                                                                                                                     <  1:143242/236‑1 (MQ=255)
                        aGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTa                                                                                                                                                                                                     <  1:672225/235‑1 (MQ=255)
                                 gTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGc                                                                                                                                                                                              >  2:399668/1‑233 (MQ=255)
                                      ccGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGc                                                                                                                                                                                                                 <  1:340470/209‑1 (MQ=255)
                                        gACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                                                                                                    <  1:837087/236‑1 (MQ=255)
                                               gAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGaca                                                                                                                                                                               <  2:736690/234‑1 (MQ=255)
                                                   aaCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGata                                                                                                                                                                         >  1:423909/1‑236 (MQ=255)
                                                    aCTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                                                                                                          <  2:806104/234‑1 (MQ=255)
                                                             ccTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                                                                                                                 >  1:910349/1‑218 (MQ=255)
                                                                  cAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtc                                                                                                                                                                       >  2:492209/1‑223 (MQ=255)
                                                                      ccgcaCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGAGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtc                                                                                                                                                                 <  1:225502/225‑1 (MQ=255)
                                                                           cACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt                                                                                                                                                                             >  1:642819/1‑208 (MQ=255)
                                                                                         ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCttt                                                                                                                                   >  1:188858/1‑236 (MQ=255)
                                                                                          cAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGcc                                                                                                                                      <  1:367344/232‑1 (MQ=255)
                                                                                               ctgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGACGTAGCACCCGGGCCTTTGATg                                                                                                                               >  2:867913/1‑234 (MQ=255)
                                                                                                     aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                                                                                                          <  2:712976/185‑1 (MQ=255)
                                                                                                               gCTGATAGAGCTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGTTAGGGCCGCAAAGCGGCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCCTCACACCtt                                                                                                                                                                                    <  2:695158/164‑1 (MQ=38)
                                                                                                                cTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                                                                                                            <  1:588057/236‑1 (MQ=255)
                                                                                                                        aGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                                                                                                    <  1:701598/236‑1 (MQ=255)
                                                                                                                          tACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                                                                                                    <  1:857658/154‑1 (MQ=255)
                                                                                                                                  gggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt                                                                                            >  2:1066743/1‑234 (MQ=255)
                                                                                                                                    gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg                                                                                           >  2:988081/1‑233 (MQ=255)
                                                                                                                                       gctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTTGCGTTGAAGc                                                                                       >  2:962364/1‑234 (MQ=255)
                                                                                                                                         ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                                                                                                         >  1:1010990/1‑182 (MQ=255)
                                                                                                                                          tCGCAGCCCTATCACGACACAAGGCCGATCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATGGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTGTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGAGACGTTGAAGCCa                                                                                     >  1:1068473/1‑233 (MQ=255)
                                                                                                                                            gCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAgc                                                                                   <  2:910348/233‑1 (MQ=255)
                                                                                                                                               gCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                                                                                                          <  2:536392/143‑1 (MQ=255)
                                                                                                                                                    aTCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCt                                                                              <  1:130604/230‑1 (MQ=255)
                                                                                                                                                             caAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                                                                                                    >  2:883847/1‑119 (MQ=255)
                                                                                                                                                             caAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                                                                                                    <  1:883848/119‑1 (MQ=255)
                                                                                                                                                             caAGGCCGCTCCTGCAGAGAATTTATTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCACAGCGGCCCCCTCGGCTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt                                                                                                                                                                             >  2:1096111/1‑126 (MQ=38)
                                                                                                                                                              aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAgc                                                                                   <  2:773683/215‑1 (MQ=255)
                                                                                                                                                                 gCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGt                                                                                                                                                           <  2:128289/140‑1 (MQ=255)
                                                                                                                                                                    gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCtt                                                                                                                                    <  2:741119/160‑1 (MQ=255)
                                                                                                                                                                     cTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCt                                                       <  1:245433/236‑1 (MQ=255)
                                                                                                                                                                      tCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCt                                                       <  1:1080140/235‑1 (MQ=255)
                                                                                                                                                                        cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCa                                                    >  1:870520/1‑236 (MQ=255)
                                                                                                                                                                        cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCa                                                    >  1:852398/1‑236 (MQ=255)
                                                                                                                                                                           cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCt                                                                                                                                                                                             >  2:252273/1‑96 (MQ=255)
                                                                                                                                                                               gaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCAcc                                                                       <  1:364586/210‑1 (MQ=255)
                                                                                                                                                                                  tttCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCt                                                                              <  2:1047246/200‑1 (MQ=255)
                                                                                                                                                                                            tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTa                                        <  1:260123/228‑1 (MQ=255)
                                                                                                                                                                                                   gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGTGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCTAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCAc                          <  1:634203/235‑1 (MQ=255)
                                                                                                                                                                                                   gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCAc                          <  1:407099/235‑1 (MQ=255)
                                                                                                                                                                                                      cGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTGTCAGCTGGCAGAAATTACTGCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTTGTTCAGGTTCAGCgtggtg                                                                                                    <  2:169601/158‑1 (MQ=255)
                                                                                                                                                                                                       ggCCTTGTGTCGTGATAGCGCCGCAAAGCGGCACCCCCGGTTTTCAGCTGGCAGAAATTACTACTGc                                                                                                                                                                                              <  2:1063168/67‑1 (MQ=255)
                                                                                                                                                                                                         ccTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCAACGTTGAAGGTAGccg                                    >  1:640398/1‑218 (MQ=255)
                                                                                                                                                                                                         ccTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGACAGCGCGCGCGCCACCCCATGGGTCCTTGTCTTccc                                                   >  2:330991/1‑202 (MQ=255)
                                                                                                                                                                                                           ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCtt                    >  2:724288/1‑233 (MQ=255)
                                                                                                                                                                                                           ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTgg                  >  1:613277/1‑235 (MQ=255)
                                                                                                                                                                                                           ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTgg                  >  2:790415/1‑235 (MQ=255)
                                                                                                                                                                                                           ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGt                 >  1:92869/1‑236 (MQ=255)
                                                                                                                                                                                                           ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCACTGGGTCCt                                                            >  1:88235/1‑193 (MQ=255)
                                                                                                                                                                                                           ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCGTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGGAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCctggtggtg                                                                                                 >  2:1085133/1‑156 (MQ=255)
                                                                                                                                                                                                              tgtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTCCTTCTGCTTCACAACTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCGCCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTg                <  1:284164/234‑1 (MQ=255)
                                                                                                                                                                                                                tCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGa                                             <  2:852397/203‑1 (MQ=255)
                                                                                                                                                                                                                tCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGacac            <  1:243140/236‑1 (MQ=255)
                                                                                                                                                                                                                   tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc                                                                                                                             <  2:902859/120‑1 (MQ=255)
                                                                                                                                                                                                                      tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCa                                                    <  1:354584/190‑1 (MQ=255)
                                                                                                                                                                                                                      tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCa                                                    >  2:354583/1‑190 (MQ=255)
                                                                                                                                                                                                                      tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCa      >  1:606251/1‑236 (MQ=255)
                                                                                                                                                                                                                      tAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGGGGTGGTGGCGTTTAAGCCAgcg                                                                                  >  2:257344/1‑160 (MQ=255)
                                                                                                                                                                                                                        gggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTc                            <  1:985435/212‑1 (MQ=255)
                                                                                                                                                                                                                        ggaccGCAAAGAGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCg                                                                            <  2:474970/161‑1 (MQ=255)
                                                                                                                                                                                                                         ggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTc                            <  1:244133/211‑1 (MQ=255)
                                                                                                                                                                                                                          gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTAGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGa               <  1:887704/223‑1 (MQ=255)
                                                                                                                                                                                                                           ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGTAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGTAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGCAGCCAGTGCGCTCGCCACCCCATGCGTCCTTGCCTTCAcc                                                  >  2:187191/1‑187 (MQ=255)
                                                                                                                                                                                                                           ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGt  >  1:843385/1‑235 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                 
CAGTGATCATTCGCAAGGAGTCCAAGCTGGACCGTTTCCCGACGCCTGAGCAACTGGATGACCTGCCAAGCCGCCCACGCGGTGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCACCGTTGAAGGTAGCCTTGAACGTCACCGGCTTGGTGACACCGTGCATGGT  >  AE015451/5674523‑5674975

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