Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 1,831,780 T→C 6.7% L220P (CTG→CCG)  yfiP → conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,831,7800TC 6.7% 200.9 / 4.5 60L220P (CTG→CCG) yfiPconserved protein
Reads supporting (aligned to +/- strand):  ref base T (26/30);  new base C (2/2);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.35e-01

GCGCAAGATGTTCCGCAAAAGCCCGTACCTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATCCTTAGATC  >  AE015451/1831550‑1832010
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      
gcgcAAGATGTTCCGCAAAAGCCCGTACCTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCg                                                                                                                                                                                                                                    >  1:958622/1‑235 (MQ=255)
     ggATGTTCCGCAAAAGCCCGTACCTGGACCGCTGCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                                                                                                              <  1:507296/235‑1 (MQ=255)
           tCCGCAAAAGCCCGTACCTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCa                                                                                                                                                                                                                         <  1:987092/235‑1 (MQ=255)
             cGCAAAAGCCCGTACCTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAgg                                                                                                                                                                                                                       <  1:935084/235‑1 (MQ=255)
                          aCCTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAg                                                                                                                                                                                                         >  1:289620/1‑236 (MQ=255)
                            cTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGc                                                                                                                                                                                                        >  1:694777/1‑235 (MQ=255)
                                  cGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                                                                                                              >  2:1223976/1‑207 (MQ=255)
                                  cGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                                                                                                              <  1:1223979/207‑1 (MQ=255)
                                       cccGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGCGCGGCAAGCGGCGGCTGCCGCTTGACGAGc                                                                                                                                                                                                                          >  1:890619/1‑206 (MQ=255)
                                       cccGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGat                                                                                                                                                                                                                              >  1:326776/1‑201 (MQ=255)
                                       cccGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAgccgcc                                                                                                                                                                                                                                     >  2:1141535/1‑195 (MQ=255)
                                       cccGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTATACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATTTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCACCCAGCGTTTGCATAc                                                                                                                                                                                              >  1:648769/1‑234 (MQ=255)
                                                gAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCa                                                                                                                                                                                                 <  2:669040/222‑1 (MQ=255)
                                                 aGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGt                                                                                                                                                                                   >  2:540281/1‑235 (MQ=255)
                                                                   aTGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTg                                                                                                                                                                <  1:1135203/236‑1 (MQ=255)
                                                                     gTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGttt                                                                                                                                                                                                    <  1:941174/198‑1 (MQ=255)
                                                                        gcgcTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATa                                                                                                                                                                                               <  1:165499/200‑1 (MQ=255)
                                                                         cgcTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAg                                                                                                                                                                                                         <  1:200833/189‑1 (MQ=255)
                                                                         cgcTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCtt                                                                                                                                                             <  2:512057/233‑1 (MQ=255)
                                                                         cgcTACCGCCGGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTg                                                                                                                                                           <  1:832601/235‑1 (MQ=255)
                                                                            tACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCtt                                                                                                                                                                                           <  1:782977/200‑1 (MQ=255)
                                                                               cGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATa                                                                                                                                                                                               <  2:344538/193‑1 (MQ=255)
                                                                                gCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATa                                                                                                                                                                                               <  2:773228/192‑1 (MQ=255)
                                                                                   tGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATa                                                                                                                                                    <  2:615670/232‑1 (MQ=255)
                                                                                    gcggcgCTCGAAGCGCGATGATCATTTCGGTCCGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGAc                                                                                                                                                                                                                              <  2:314361/157‑1 (MQ=255)
                                                                                    gcggcgCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACGTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTAGCATACCTTACTATAGTCCAAGGGGCACAACGCCt                                                                                                                                                                 >  2:1208855/1‑218 (MQ=255)
                                                                                      ggcgCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCAt                                                                                                                                                                                                <  2:597101/185‑1 (MQ=255)
                                                                                           tCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccac                                                                                                                                                                                                            <  2:74576/168‑1 (MQ=255)
                                                                                            ccAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGAcacccacc                                                                                                                                                                                                           <  2:598035/168‑1 (MQ=255)
                                                                                                              tCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCgag                                                                                                                      >  1:1172747/1‑235 (MQ=255)
                                                                                                                     ggCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGGGGGCGAGAc                                                                                                                    >  1:290932/1‑230 (MQ=255)
                                                                                                                          aaGTTGCGGCGATGTGCCTGGACCTTGCTGTCGATGGCCAGTCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGa                                                                                                                                                                                                                               >  2:368205/1‑118 (MQ=255)
                                                                                                                             ttGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCCGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGAtt                                                                                                        <  1:213459/234‑1 (MQ=255)
                                                                                                                                cggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCtt                                                                                                   <  1:1259983/236‑1 (MQ=255)
                                                                                                                                cggcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCt                                                                                                    <  1:116024/235‑1 (MQ=255)
                                                                                                                                  gcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTc                                                                                                                                                                                          <  2:865318/147‑1 (MQ=255)
                                                                                                                                  gcgATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTg                                                                                                  <  1:540284/235‑1 (MQ=255)
                                                                                                                                                      tGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTc                                                                               >  1:719566/1‑234 (MQ=255)
                                                                                                                                                                cAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCGCCCTTGACCGTGCt                                                                                >  1:973831/1‑223 (MQ=255)
                                                                                                                                                                          aGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCGCGCGGCGGCTGCCGGTTGACGAGCAGGACGACCCCCACCAgg                                                                                                                                                                                                        >  2:189367/1‑92 (MQ=255)
                                                                                                                                                                           ggCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGt                                                                                   >  2:841094/1‑209 (MQ=255)
                                                                                                                                                                                tGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGACGAGCa                                                                                                                                                                                                                         >  2:1282070/1‑70 (MQ=255)
                                                                                                                                                                                tGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCAGCCGCCGCTTGACGAGCa                                                                                                                                                                                                                         <  1:1282073/70‑1 (MQ=255)
                                                                                                                                                                                  gATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCa                                                                                                                >  2:575109/1‑173 (MQ=255)
                                                                                                                                                                                   aTGCCGATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAgggg                                                                                                                                                                           <  2:175390/113‑1 (MQ=255)
                                                                                                                                                                                                    ttCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCCGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGCtt                                                            >  1:249032/1‑207 (MQ=255)
                                                                                                                                                                                                      ccgCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGt                              <  1:407145/233‑1 (MQ=255)
                                                                                                                                                                                                      cAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGt                              <  1:1212807/235‑1 (MQ=255)
                                                                                                                                                                                                       aGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGt                              <  2:652133/234‑1 (MQ=255)
                                                                                                                                                                                                            gCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGGGGGCGAGACGGCa                                                                                                                >  2:289546/1‑147 (MQ=255)
                                                                                                                                                                                                               ttACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCt                                                                                          >  2:1129692/1‑166 (MQ=255)
                                                                                                                                                                                                                   cTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCACAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGATCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAAc                <  1:1208858/236‑1 (MQ=255)
                                                                                                                                                                                                                         ggCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTcccgg                                                       >  2:650196/1‑189 (MQ=255)
                                                                                                                                                                                                                          gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCggg                                                  >  2:63957/1‑195 (MQ=255)
                                                                                                                                                                                                                          gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATcc         >  1:477476/1‑236 (MQ=255)
                                                                                                                                                                                                                          gCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATc          >  1:947034/1‑235 (MQ=255)
                                                                                                                                                                                                                              gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATCATGACGGAGCGGGCAGCCAGGGTGGGCGAGACgg                                                                                                                  >  2:1013941/1‑127 (MQ=255)
                                                                                                                                                                                                                              gcggcgGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATCCTTAg     >  1:1115559/1‑236 (MQ=255)
                                                                                                                                                                                                                                ggcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATCCTTAGa    >  1:1011148/1‑235 (MQ=255)
                                                                                                                                                                                                                                 gcggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATCCTTAGATc  <  1:63957/236‑1 (MQ=255)
                                                                                                                                                                                                                                  cggcTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGc                       <  2:578519/214‑1 (MQ=255)
                                                                                                                                                                                                                                      tGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGcttcc                                           >  1:1175314/1‑190 (MQ=255)
                                                                                                                                                                                                                                      tGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAAc                <  2:694774/217‑1 (MQ=255)
                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                      
GCGCAAGATGTTCCGCAAAAGCCCGTACCTGGACCGCTTCCCGGTGTTGAGCCTGCAAGCCGAGCAGATGTCGCGCTACCGCCTGCGGCGCTCCAAGCGCGATGATCATTTCTGTACGGCAGAAGTTGCGGCGATGTGCCTGGACCTTGCTGGCGATGGCCAGGCTTCGCAGGCGCTGGATGCCTATCTGGATGTGTTCAGCCTGCATTACCTGAGCGGCAAGCGGCGGCTGCCGCTTGACGAGCAGGACGACACCCACCAGCGTTTGCATACCTTCCTATAGTCCAAGGGGCACAACCCCTGCTTTGGTTCATAATGACGGAGCTGGCAGCCAGGGTGGGCGAGACGGCAGGGCGATTGGCTTGATCTCCCTTGACCGTGCTCGGCATCCTTGGCGCCGATTCCCCGCCGGGCGCTTCCCTTCCCCTCGAGTGCCTGGCACAGAACAGGATCCTTAGATC  >  AE015451/1831550‑1832010

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