Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 2,087,635 T→C 11.5% intergenic (+448/+470) PP_1865 → / ← aroH transposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154512,087,6350TC11.5% 179.8 / 13.6 61intergenic (+448/+470)PP_1865/aroHtransposase/Phospho‑2‑dehydro‑3‑deoxyheptonate aldolase
Reads supporting (aligned to +/- strand):  ref base T (35/19);  new base C (3/4);  total (38/23)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.83e-01

CGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTT  >  AE015451/2087438‑2087823
                                                                                                                                                                                                     |                                                                                                                                                                                              
cgcTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCa                                                                                                                                                                                     >  2:1179183/1‑209 (MQ=255)
 gcTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAgc                                                                                                                                                             >  2:952887/1‑232 (MQ=255)
    tAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATACGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTTGCTGTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCaggtgg                                                                                                                                                                   >  1:264661/1‑223 (MQ=255)
                cAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGc                                                                                                                                                                                      >  1:528526/1‑192 (MQ=255)
                      cGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCagga                                                                                                                                                                     >  1:801259/1‑203 (MQ=255)
                                        aTCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGGTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGTGGAGGGAGGTGACGGGcct                                                                                                                               <  1:643509/223‑3 (MQ=255)
                                          cggccggaACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTg                                                                                                                                                                                <  1:202238/172‑1 (MQ=255)
                                            gccggcACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTACTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                               <  1:562512/235‑1 (MQ=255)
                                                 cACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCggagggag                                                                                                                                          <  1:292286/203‑1 (MQ=255)
                                                              gtCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAg                                                                                                                                                              <  1:768629/169‑1 (MQ=255)
                                                                       gTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGAt                                                                                                                                                                                  <  2:1147618/141‑1 (MQ=255)
                                                                         tCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGAt                                                                                                                           >  2:303665/1‑194 (MQ=255)
                                                                            gCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTggg                                                                                                         >  2:1207271/1‑208 (MQ=255)
                                                                                ggTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAg                                                                                                                                                              >  1:152081/1‑152 (MQ=255)
                                                                                   cATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCgg                                                                      <  1:1179186/235‑1 (MQ=255)
                                                                                     tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                             >  1:684886/1‑180 (MQ=255)
                                                                                     tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                             <  2:684883/180‑1 (MQ=255)
                                                                                     tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATg                                                                                                                              >  1:1053373/1‑179 (MQ=255)
                                                                                     tCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGGCTGCAgcg                                                                                  >  1:849487/1‑221 (MQ=255)
                                                                                        ggATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                               <  1:28787/191‑1 (MQ=255)
                                                                                               gATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                         >  1:541889/1‑236 (MQ=255)
                                                                                                  atCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAg                                                        <  2:541886/234‑1 (MQ=255)
                                                                                                    cGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGTGCTTGCGGAGGGAGGTGa                                                                                                                                      >  1:132985/1‑156 (MQ=255)
                                                                                                      ccGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGGCGTTGCTCTAGCTAGTGCCTTTGCTTAAGCTCTTGCTTCTTAGCGGCCGGTAGTTCTGgcggcgttgc                                                                                                                                                                                   >  1:980700/1‑105 (MQ=16)
                                                                                                       cGGCAGGGCTGCCCTCGGCCTGGTGTACCGGGGAGAACGCCGTTGCTCTGGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTg                                                                                                                                                   <  2:1147212/140‑1 (MQ=255)
                                                                                                            tGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCt                                                                                        >  2:184781/1‑192 (MQ=255)
                                                                                                             ggCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAg                                                                                                                                                                                    <  2:1130668/101‑1 (MQ=255)
                                                                                                                tGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGag                                        <  1:184781/236‑1 (MQ=255)
                                                                                                                  cccTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcgtgctccc                                                                         <  2:821251/203‑8 (MQ=14)
                                                                                                                  cccTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCATGGATGGGGCCTGCAGcgcgtgctccc                                                                         >  1:821254/1‑196 (MQ=14)
                                                                                                                    atCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTATTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                             <  1:344861/148‑1 (MQ=255)
                                                                                                                     tCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                      >  1:1044378/1‑170 (MQ=255)
                                                                                                                       ggCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGa                                                         >  1:1188950/1‑212 (MQ=255)
                                                                                                                          cTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGGGTGGAGGCATTGCGGAGGGAGATGACGGGCATGGATGccc                                                                                                                       >  1:989333/1‑149 (MQ=17)
                                                                                                                               gtACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAgg                                                                                                                                                                                            >  1:865053/1‑75 (MQ=255)
                                                                                                                                tACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGcgccgc                                                                                                                                                                                      >  2:242638/1‑80 (MQ=255)
                                                                                                                                tACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTg                                                                                                            >  2:927381/1‑154 (MQ=255)
                                                                                                                                     ggAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGGGAGGGGACCCCGGAGCGTAGCGTAggg                            >  2:888717/1‑227 (MQ=255)
                                                                                                                                        agaaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGccc                                                                                                      >  2:408549/1‑151 (MQ=255)
                                                                                                                                          aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAgcg                                                                                                                                                            >  1:1022812/1‑96 (MQ=255)
                                                                                                                                          aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATgg                                                                                                                             >  1:132198/1‑127 (MQ=255)
                                                                                                                                          aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                >  2:1018300/1‑140 (MQ=255)
                                                                                                                                          aaCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAgc                                                                                                                <  1:1018303/140‑1 (MQ=255)
                                                                                                                                           aCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTcg                                                                                                                                                                        >  1:527441/1‑82 (MQ=255)
                                                                                                                                           aCGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGcccaa                                                                                                    >  2:942536/1‑148 (MQ=255)
                                                                                                                                              ccGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGTGGAGGCATTGCgg                                                                                                                                                <  1:788218/104‑1 (MQ=17)
                                                                                                                                               cGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgt                                     >  2:700757/1‑208 (MQ=40)
                                                                                                                                               cGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTCGCGTAGGGGCCGGATGATGGGAGcc           >  2:974440/1‑234 (MQ=255)
                                                                                                                                               cGTTGCTCTTGCTCTTGCGCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGcgcg                                                                                <  1:976436/165‑1 (MQ=255)
                                                                                                                                                gTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGAc        <  1:555439/236‑1 (MQ=40)
                                                                                                                                                 ttgctcttgctcttgctcTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAg                                                                                                                 <  1:632401/132‑1 (MQ=255)
                                                                                                                                                  tgctcttgctcttgctctTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGTCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGACTGCAGCGCGGTCCTccc                                                                        >  1:97878/1‑170 (MQ=255)
                                                                                                                                                   gctcttgctcttgctcttGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATgggg              >  2:555436/1‑226 (MQ=39)
                                                                                                                                                      ccgaattcgcttctcgctTTTGATCTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGtttt  >  1:882625/12‑236 (MQ=11)
                                                                                                                                                       cgaattcgcttctcgctTTTGATCTTGCTTCTAAGCGCGCGGTAGTCCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGa              <  2:882622/213‑1 (MQ=11)
                                                                                                                                                                  cttgctTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGTGCGTGCggagc                                                                                                                                             >  2:800481/1‑86 (MQ=255)
                                                                                                                                                                            cTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                  >  1:1277325/1‑118 (MQ=32)
                                                                                                                                                                            cTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAg                                                                                                  <  2:1277322/118‑1 (MQ=32)
                                                                                                                                                                                  cTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAggg                                                  >  1:614367/1‑160 (MQ=31)
                                                                                                                                                                                    tCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGagcgtagc                                  <  1:303666/174‑1 (MQ=31)
                                                                                                                                                                                        aGCGCGCGGTAGTTCAGGCGACGTAGATTGCGAGTGCAGGAGGc                                                                                                                                                                  >  2:461390/1‑44 (MQ=11)
                                                                                                                                                                                                     |                                                                                                                                                                                              
CGCTTAGAAGCAAGATCAAAAGCGAGAAGCGAATTCGGTGATCGGCCGGCACCGCCGATGCATTCGGCGAGGTTCTGCCCGGTCATCTGGATAAGGATATCGCCGGCATGGCTGCCCTCGGCCTGGTGTACCTGGAAGAACGCCGTTGCTCTTGCTCTTGCTCTTGCTTTGGCTTTTGCTTCTAAGCGCGCGGTAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGAGCGCAGGGCTTGCGGAGGGAGGTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGTGAGGGGACCCCGGAGCGTAGCGTAGGGGCCGGATGATGGGAGCCGACGGTTTT  >  AE015451/2087438‑2087823

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