Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,962,159:1 +C 7.4% coding (422/1833 nt) PP_1760 → ABC transporter, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,962,1591.C 7.4% 218.8 / 11.3 54coding (422/1833 nt)PP_1760ABC transporter, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (26/24);  new base C (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

TCTGCTCGGCCATGGCGATGCTGCGCTGAAGTGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCTACTTCACCGTGCAGTTGGGCAA  >  AE015451/1961935‑1962384
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
tCTGCTCGGCCATGGCGATGCTGCGCTGAAGTGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGt                                                                                                                                                                                                                        <  1:818470/236‑1 (MQ=255)
tCTGCTCGGCCATGGCGATGCTGCGCGGAAGTGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGt                                                                                                                                                                                                                        <  1:373082/236‑1 (MQ=255)
     tCGGCCATGGCGATGCTGCGCTGAAGTGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGt                                                                                                                                                                                                                        <  1:95001/231‑1 (MQ=255)
          cATGGCGATGCTGCGCTGAAGTGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCt                                                                                                                                                                                                              >  1:300491/1‑236 (MQ=255)
                               tGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGc                                                                                                                                                                                                                                 <  2:829354/196‑1 (MQ=255)
                                         aTCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGTCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGACTCAAGCGTATTTCGCTGCAagag                                                                                                                                                                                                                         >  2:187320/1‑194 (MQ=255)
                                              gtgccGTCCAGCTGGCAGGGGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtgg                                                                                                                                                                                 >  2:81664‑M2/2‑226 (MQ=255)
                                                          tGGCGGGGGGCGGCCGGCTACATCGGCCTGATGCTGGCGCTCCCCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGt                                                                                                                                                                                                                        <  2:346014/178‑1 (MQ=255)
                                                                 tGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGCCACCGACCTGGACACCCCCGCCaa                                                                                                                                                        >  1:256739/1‑235 (MQ=255)
                                                                     aGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCt                                                                                                                                                              >  2:656231/1‑225 (MQ=255)
                                                                                           cTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtggtcaccgtgccgctgcccaggctttcgtactctttcagcgaaatacgcctgt                                                                                                                              <  1:156711‑M2/236‑55 (MQ=255)
                                                                                           cTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtggtcaccgtgccgctgcccaggctttcgtactctttcagc                                                                                                                                           <  1:156709‑M2/223‑42 (MQ=255)
                                                                                           cTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCtt                                                                                                                             >  1:558381/1‑236 (MQ=255)
                                                                                                      ccTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGcc                                                                                                                   <  1:438706/235‑1 (MQ=255)
                                                                                                      ccTGGGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGGTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtgg                                                                                                                                                                                 >  2:722554‑M2/1‑170 (MQ=255)
                                                                                                                   cTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACcacccaccg                                                                                                                                                                                   >  2:432977/1‑157 (MQ=255)
                                                                                                                   cTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGCCAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACt                                                                                                        >  2:778622/1‑233 (MQ=255)
                                                                                                                   cTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGAc                                                                                                     >  1:940555/1‑236 (MQ=255)
                                                                                                                     ggCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGAcc                                                                                                    >  1:1030842/1‑235 (MQ=255)
                                                                                                                       cTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCg                                                                                                   >  2:984962/1‑234 (MQ=255)
                                                                                                                        tGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGc                                                                                                                                                                                                                                 >  2:296041/1‑107 (MQ=255)
                                                                                                                        tGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGc                                                                                                                                                                                                                                 <  1:296044/107‑1 (MQ=255)
                                                                                                                         gCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCAc                                                                                               >  1:1008529/1‑236 (MQ=255)
                                                                                                                                              cAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCtggt                                                                                                                       <  1:213028/191‑1 (MQ=255)
                                                                                                                                                    aaGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtgg                                                                                                                                                                                 <  1:843728‑M2/127‑4 (MQ=255)
                                                                                                                                                    aaGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtgg                                                                                                                                                                                 >  2:843725‑M2/1‑124 (MQ=255)
                                                                                                                                                    aaGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtgg                                                                                                                                                                                 >  2:4044‑M2/1‑124 (MQ=255)
                                                                                                                                                        tGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAg                                                                >  1:975669/1‑236 (MQ=255)
                                                                                                                                                         gTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCa                                                                 >  2:400770/1‑234 (MQ=255)
                                                                                                                                                                         aaGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCAC                                                                                                                                                                                                  <  2:238984‑M1/89‑1 (MQ=255)
                                                                                                                                                                                  ctctccCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtggtcaccgtgccgctgcccaggctttcgtactctt                                                                                                                                                >  2:156707‑M2/2‑95 (MQ=255)
                                                                                                                                                                                    cTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGcc                                                                                                                   <  1:223825/157‑1 (MQ=255)
                                                                                                                                                                                    cTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAAcc                                    <  1:367651/236‑1 (MQ=255)
                                                                                                                                                                                     tACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCAggtgggtggtcaccgtgccgctgcccaggctttcgtactctttcagcgaaatacgcttgaggcgctcgatcaggcggatgcgcaggcgg                                                                                                 >  1:168650‑M1/1‑85 (MQ=255)
                                                                                                                                                                                       ccGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtggtcaccgtgccgctgcccaggctttcgtactctttcagcgaaatacgcttgaggcgctcgatcaggcggatgcgcaggcggta                                                                                               >  1:168650‑M2/1‑89 (MQ=255)
                                                                                                                                                                                        cGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCtggtcaccgtgccgctgcccaggctttcgtactctt                                                                                                                                                >  2:156709‑M2/1‑89 (MQ=255)
                                                                                                                                                                                        cGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCt                                >  1:303417/1‑236 (MQ=255)
                                                                                                                                                                                        cGCCAGCGCATCCGCCTGATCGAGCGCCTCAGGCGTATTTCGCTGAAAGAGTACGAAGGCCTGGGCAGCGGCACGGTGACCACCCGCCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCt                                                         <  2:350992/211‑1 (MQ=255)
                                                                                                                                                                                         tccTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCAcgg                                                                                             <  1:84517/173‑1 (MQ=255)
                                                                                                                                                                                         gCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTg                               <  1:834816/236‑1 (MQ=255)
                                                                                                                                                                                          ccTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAAccc                                   >  1:493692/1‑231 (MQ=255)
                                                                                                                                                                                             gcgcATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCATCCCGCTGGTGa                           >  1:885032/1‑236 (MQ=255)
                                                                                                                                                                                                cATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCCCCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCt                        <  1:326351/236‑1 (MQ=255)
                                                                                                                                                                                                   ccGCCTGATCGGGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCCCCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCAcggcg                                                                                           <  2:336803/166‑1 (MQ=255)
                                                                                                                                                                                                     gCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCCACCCg                                  >  1:939852/1‑221 (MQ=255)
                                                                                                                                                                                                         gATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCTACTTCa                  <  2:951074/233‑1 (MQ=255)
                                                                                                                                                                                                                 gcCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCt                                                                                                                                                              <  2:170295/85‑1 (MQ=255)
                                                                                                                                                                                                                 gcCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCt                                                                                                                                                              >  1:170297/1‑85 (MQ=255)
                                                                                                                                                                                                                   cTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCTACTTCACCGTGCAGTTgg      <  1:1023445/235‑1 (MQ=255)
                                                                                                                                                                                                                    tCAGGCGTATTTCGCGGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCTACTTCACCGTGCAGTTgg      <  2:818177/234‑1 (MQ=255)
                                                                                                                                                                                                                    acaggcgtaTTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCAggtgggtggtcaccgtgccgctgcccaggctttcgtactctttcagcggaaatacgcttgaggcgctcgatcaggcggatgcgcaggcggtagacgatgtccttggccagcccggcgaacagcttggactggataacgttgaacgccagtgcagccaggcgcaggcacagggtgagcgccag    <  1:156711‑M1/232‑183 (MQ=255)
                                                                                                                                                                                                                        gCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCTACTTCACCGTGCAGTTGGGCaa  >  2:356836/1‑234 (MQ=255)
                                                                                                                                                                                                                            agttCGCCGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCGCCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGc                                                                        <  2:142706/158‑1 (MQ=255)
                                                                                                                                                                                                                            aTTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAAcc                                    <  2:1088809/196‑1 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                 
TCTGCTCGGCCATGGCGATGCTGCGCTGAAGTGGATGAACCATCTGCTGCCGTCCAGCTGGCAGGTGGCAGCCGGCTACATCGGCCTGATGCTGGCGCTCACCCTGTGCCTGCGCCTGGCTGCACTGGCGTTCAACGTTATCCAGTCCAAGCTGTTCGCCGGGCTGGCCAAGGACATCGTCTACCGCCTGCGCATCCGCCTGATCGAGCGCCTCAAGCGTATTTCGCTGAAAGAGTACGAAAGCCTGGGCAGCGGCACGGTGACCACCCACCTGGTCACCGACCTGGACACCCTCGACAAGTTCGTTGGCGAAACCCTCAGCCGCTTCCTGGTGGCCATGCTGACACTGACCGGCACGGCGGCCATCCTGATCTGGATGCACTGGCAGCTGGCCTTGCTGATCTTGCTGTTCAACCCGCTGGTGATCTACTTCACCGTGCAGTTGGGCAA  >  AE015451/1961935‑1962384

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