Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 1,419,549 +C intergenic (‑38/+98) PP_1242 ← / ← PP_1244 conserved protein of unknown function/conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154511,419,5481.C96.8% 110.8 / ‑2.9 31intergenic (‑37/+99)PP_1242/PP_1244conserved protein of unknown function/conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/17);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CATCCCCAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCAGAGTACGTACT  >  AE015451/1419251‑1419831
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                           
catcCCCAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggcct                                                                                                                                                                                                                                                                                          <  4:179440/301‑3 (MQ=255)
catcCCCAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTGGTGTTTTCGGTGGGCATTTTTTTTTTTCCCTGGGTTCTGGGTTTCCTTggcgggccc                                                                                                                                                                                                                                                                                          >  2:2746/1‑299 (MQ=255)
  cccccAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggcctct                                                                                                                                                                                                                                                                                        <  1:605386/300‑5 (MQ=255)
   ccccAACTCGCCACAATCCTGCAGCGCCTTGGCAATCGAGGTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGggctggcctctt                                                                                                                                                                                                                                                                                       <  2:198835/301‑6 (MQ=255)
                      agcagcGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTccc                                                                                                                                                                                                                                                                     <  4:17301/301‑1 (MQ=255)
                                  cccATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACtgtg                                                                                                                                                                                                                                                         <  1:526205/300‑1 (MQ=255)
                                     tcttccgagctgTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAg                                                                                                                                                                                                                                                      <  2:147738/290‑1 (MQ=255)
                                        gcgcTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGc                                                                                                                                                                                                                                                   <  4:2746/301‑1 (MQ=255)
                                                gTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAgagatcggaaga                                                                                                                                                                                                                                           >  4:147768/1‑290 (MQ=255)
                                                          cccccgATCAGTTTCTGGACGCGGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGAGGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGAcc                                                                                                                                                                                                                                 <  2:709190/296‑1 (MQ=255)
                                                             cAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCtg                                                                                                                                                                                                                              >  4:359193/1‑301 (MQ=255)
                                                                aTCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTgg                                                                                                                                                                                                                           <  4:410139/301‑1 (MQ=255)
                                                                 ccggTTTGTGGTCGCCGTCCACCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTggg                                                                                                                                                                                                                          <  2:696534/298‑1 (MQ=255)
                                                                       tCTGTTCGCCGTCCATCAGGCCGGGGAGGATCAGGTTGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCggg                                                                                                                                                                                                                    <  2:359120/301‑1 (MQ=255)
                                                                                          gCCGGGGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCGTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCt                                                                                                                                                                                                 <  3:487043/301‑1 (MQ=255)
                                                                                                   gtcgtgtgctcttccgatctccATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTa                                                                                                                                                                                        <  4:681598/281‑1 (MQ=255)
                                                                                                      cAGGTCGCGCACGTAACCCATCAGTGCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGGTGCGGAACAACGGGTGAGTGTTGTTTTCGCATTGGTAGGAGCAGGTCTTGCCGGCTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATg                                                                                                                                                                                     <  4:390121/301‑1 (MQ=255)
                                                                                                             tgtgctcttccgatctgTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGgtg                                                                                                                                                                              <  3:52370/285‑1 (MQ=255)
                                                                                                                      cccATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGGTGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGctct                                                                                                                                                                     <  1:427919/301‑1 (MQ=255)
                                                                                                                       ccATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGGGCGTGACGGCGGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCttagatccgaggtgcgccgggg                                                                                                                                                                    >  2:681422/1‑280 (MQ=255)
                                                                                                                        cgtgaggtCGGAAGCCTGGTCGCGGGGGTCCGGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGCGTGGCTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTTTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTATGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCCTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCtgt                                                                                                                                                                   <  1:469366/294‑1 (MQ=255)
                                                                                                                             gTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGgtgagatcggaagagcgtc                                                                                                                                                              >  1:52369/1‑285 (MQ=255)
                                                                                                                                   ttcctgcagtgtgctcttccgttctggggATGCTGGGTTCGATGCGTAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTTGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCg                                                                                                                                                        <  4:219168/276‑1 (MQ=255)
                                                                                                                                                                     gggTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGTGAGCGGGAGTGCCCGCGACGAGGCCTGACCTGCTAATGCACGGTGTGAGGCTATGTGGGTAAATGCGACCGCATCTTCAACTCCATCACCGCCGACTCGAc                                                                                                                      >  2:688658/1‑301 (MQ=255)
                                                                                                                                                                           gATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGccc                                                                                                                >  3:186506/1‑301 (MQ=255)
                                                                                                                                                                                  aaCAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTGTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGAGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTAcc                                                                                                         <  4:688840/301‑1 (MQ=255)
                                                                                                                                                                                                gttgttTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCCACTGCATGAACGCCGACTCGACCAGCACCTGTACCGCTTTCATTTCATg                                                                                           >  1:351468/1‑301 (MQ=255)
                                                                                                                                                                                                          ccTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGcc                                                                                 >  3:231933/1‑301 (MQ=255)
                                                                                                                                                                                                           cTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGGTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCt                                                                                <  4:376816/301‑1 (MQ=255)
                                                                                                                                                                                                                                          gggTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGttt                                                 >  3:478874/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                           gtttCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGTCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAgggg                >  1:548830/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                              tCCCTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGgcgcg             >  4:343412/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                ccTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCAACGGACAGGGCGCagg           >  1:221364/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                 cTGGATTCTGGGTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCagag          >  3:79010/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                         ggggTTGCCTGGGCTGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGTCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCAGAGTACGTAct  <  3:548844/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                           
CATCCCCAACTCGGCATCATCCAGCAGCGCCTTGGCCATCGCGCTCAGGTAATGCGAGGCCCAGATCAGTTTCTGGTCGCCGTCCATCAGGCCGGTGATGATCAGGTCGCGCACGTAACCCATCAGTTCGGAAGCCTGTTCGCGGGCGTCCTGGCAGGGGATGCCGGGTTCGATGCGGAACAGCGGGTGAGTGTTGTTTTCGCCTTGGTAGAAGCAGGTCTTGCCGACTGTGGTGGGTTCTGTTGTGTTGTCGGTGGGCATTGTTTTGTTTCCCTGGATTCTGGGTTGCCTGGGCTGGCTCTTCGCGGGCACGCCCGCTCCCACAGGTACTGTGCAGTGCGTGACGGCTGTGGTGACCCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCTGACCTGCTAATGCAAGGTGTGAGGCTCTGTGGGTAAATGCGACCGCATCTGCAACTGCATCAACGCCGACTCGACCAGCCCCTGTACCGCTTCCATTTCATGGATCACGGCCTGCATCACGATCGAGCTGTGGGCGGTGGGTTTGAGTTGCACGGATTGGCGGGCCACGGACAGGGCGCAGAGTACGTACT  >  AE015451/1419251‑1419831

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