Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,481,628 4 bp→TACC intergenic (+279/‑727) unknown → / → unknown Zona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,481,6280CT97.4% 119.8 / ‑5.8 38intergenic (+279/‑730)unknown/unknownZona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (15/22);  total (15/23)
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.00e+00
*11,481,6290TA94.7% 114.3 / ‑5.1 38intergenic (+280/‑729)unknown/unknownZona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (15/21);  minor base C (0/2);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*11,481,6300GC100.0% 124.2 / NA 38intergenic (+281/‑728)unknown/unknownZona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (15/23);  total (15/23)
*11,481,6310AC94.7% 108.8 / ‑5.7 38intergenic (+282/‑727)unknown/unknownZona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (15/21);  total (15/23)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GG‑AAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAATAATGCTCGATGTTATTCGGGCATTTTTATTTGTTGCGCG‑CTGATGATTTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACT  >  1/1481407‑1481815
                                                                                                                                                                                                                                 ||||                                                                                                                                                                                        
gtcTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAAt                                                                                                                                                                  >  1:850041/3‑251 (MQ=255)
          aaataaaaTAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGccc                                                                                                                                                         >  1:712035/8‑250 (MQ=255)
               aaaTAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCaa                                                                                                                                                   <  1:689263/249‑1 (MQ=255)
                            tgAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTAA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGc                                                                                                                                      <  1:484113/250‑1 (MQ=255)
                            tgAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTGGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGc                                                                                                                                      <  1:575145/250‑1 (MQ=255)
                                   ggTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGtaa                                                                                                                               >  1:875811/1‑251 (MQ=255)
                                       tttCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGtaataa                                                                                                                            >  1:521567/1‑250 (MQ=255)
                                       tttCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAAc                                                                                                                           <  1:41444/251‑1 (MQ=255)
                                          cATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAActgc                                                                                                                        <  1:498703/251‑1 (MQ=255)
                                                   cTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGAt                                                                                                               <  1:790901/251‑1 (MQ=255)
                                                           gtttgatagttttttaacgaaattGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGAt                                                                                                       >  1:297496/24‑251 (MQ=255)
                                                                      aaaaTCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGaa                                                                                            <  1:341661/251‑1 (MQ=255)
                                                                                  tGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAAc                                                                                                                           >  1:400230/1‑208 (MQ=255)
                                                                                     tACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGc                                                                                                          <  1:531319/222‑1 (MQ=255)
                                                                                     tACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATAc                                                                                                     >  1:303813/1‑227 (MQ=255)
                                                                                     tACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACg                                                                                                    <  1:613501/228‑1 (MQ=255)
                                                                                      aCTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACAtttt                                                                            >  1:739068/1‑251 (MQ=255)
                                                                                       cTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACAttttt                                                                           >  1:756241/1‑251 (MQ=255)
                                                                                           gACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTaa                                                                                                                                         <  1:723725/185‑1 (MQ=255)
                                                                                           gACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTa                                                                                        >  1:791023/1‑234 (MQ=255)
                                                                                             cAGATATAGTGATTTTCCTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCAGGCAGAATGTTATTCGGGCATTTTTATTTTTTCCGCG‑TCCATAAATTAAATCGAGAATTGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATAc                                                                                                     <  1:226758/219‑1 (MQ=255)
                                                                                                atatAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGAt                                                                                                               <  1:171824/206‑1 (MQ=255)
                                                                                                      tGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGGTTGCAGCCCGATAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGTCCCA‑AAATTTAATCAAGAAATAAAAATTTACCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATa                                                            <  1:711642/251‑1 (MQ=255)
                                                                                                        aTTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAAc                                                                                                                           >  1:799159/1‑186 (MQ=255)
                                                                                                        aTTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTc                                                                  <  1:426022/243‑1 (MQ=255)
                                                                                                        aTTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATaaa                                                          <  1:268366/251‑1 (MQ=255)
                                                                                                         ttttACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAAt                                                         >  1:784907/1‑251 (MQ=255)
                                                                                                                  aaGTGAAACTTGTCTAACTATTTTTGAGCAAGC‑TCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAAtattta                                                >  1:844796/1‑251 (MQ=255)
                                                                                                                                  aactatttttgggcaagctcaaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAgg                                  >  1:346280/22‑250 (MQ=255)
                                                                                                                                  aactatttttgagcaagctcaaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGc                                 <  1:236940/230‑1 (MQ=255)
                                                                                                                                                  gctcaaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATa                                                                                                                                                <  1:720040/119‑1 (MQ=255)
                                                                                                                                                     caaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAg                                                                                                                                       <  1:233728/128‑1 (MQ=255)
                                                                                                                                                                  cGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAAcc  <  1:367120/250‑2 (MQ=255)
                                                                                                                                                                              gCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACaa                                    <  1:828074/204‑1 (MQ=255)
                                                                                                                                                                               cccGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAatat                                                   <  1:530159/188‑1 (MQ=255)
                                                                                                                                                                                                 aTGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGc                                                                                                                                      >  1:791278/1‑87 (MQ=255)
                                                                                                                                                                                                 aTGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGc                                                                                     <  1:51778/136‑1 (MQ=255)
                                                                                                                                                                                                 aTGTTATTCGGGCATTTTTATTTTTTGCGAA‑TACCTAAATTTAACCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGc                                                                                     <  1:785375/136‑1 (MQ=255)
                                                                                                                                                                                                                                 ||||                                                                                                                                                                                        
GACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GG‑AAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAATAATGCTCGATGTTATTCGGGCATTTTTATTTGTTGCGCG‑CTGATGATTTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACT  >  1/1481407‑1481815

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