Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I1 R1 150 26.0 2053212 90.0% 1847890 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 787,882 +T intergenic (‑662/+630) h16_A0728 ← / ← h16_A0729 response regulator containing CheY‑like receiver, AAA‑type ATPase, and DNA‑binding domains/transcriptional regulator, LysR‑family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM260479787,8791.T96.8% 114.6 / ‑3.2 31intergenic (‑659/+633)h16_A0728/h16_A0729response regulator containing CheY‑like receiver, AAA‑type ATPase, and DNA‑binding domains/transcriptional regulator, LysR‑family
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base T (16/14);  total (16/15)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TCGGCGAACCCGTGCCAGTGTGTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGG  >  AM260479/787747‑788007
                                                                                                                                      |                                                                                                                                
tCGGCGAACCCGTGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAaggaaaaggaaa                                                                                                                                   >  2:359730/1‑134 (MQ=255)
tCGGCGAACCCGTGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAaggaaaaggaaa                                                                                                                                   >  2:325765/1‑134 (MQ=255)
            tGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGc                                                                                                                       >  1:118936/1‑134 (MQ=255)
                        gCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAAtt                                                                                                            <  1:272220/133‑1 (MQ=255)
                          cGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCt                                                                                                                          <  1:235391/117‑1 (MQ=255)
                          cGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCt                                                                                                                          >  2:235391/1‑117 (MQ=255)
                                                cTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCtt                                                                                   >  2:804384/1‑134 (MQ=255)
                                                         gCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACAt                                                                          >  2:290694/1‑134 (MQ=255)
                                                            tCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         <  2:836342/132‑1 (MQ=255)
                                                            tCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         >  1:836342/1‑132 (MQ=255)
                                                            tCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTg                                                                       >  2:515433/1‑134 (MQ=255)
                                                                ttCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  2:952113/133‑1 (MQ=255)
                                                                ttCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  1:804384/133‑1 (MQ=255)
                                                                ttCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  1:359730/133‑1 (MQ=255)
                                                                        gACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCtat                                                          <  1:887646/134‑1 (MQ=255)
                                                                            tatGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATAtt                                                       >  1:510226/1‑134 (MQ=255)
                                                                                         cGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCAtt                                          >  1:958097/1‑134 (MQ=255)
                                                                                           tGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCg                                        >  2:630184/1‑134 (MQ=255)
                                                                                              cAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCg                                     <  1:325765/134‑1 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCt                                                                                                         <  1:783544/54‑1 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCt                                                                                                         >  2:783544/1‑54 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAAtt                          <  1:515433/133‑1 (MQ=255)
                                                                                                          aTCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAAtt                          <  2:510226/133‑1 (MQ=255)
                                                                                                           tCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAg                        <  1:290694/134‑1 (MQ=255)
                                                                                                            ccTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    <  2:578930/89‑1 (MQ=255)
                                                                                                            ccTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCg                                                                    >  1:578930/1‑89 (MQ=255)
                                                                                                              tttGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGgcg                     >  2:269563/1‑134 (MQ=255)
                                                                                                                tGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCa                   >  1:748753/1‑134 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTg                                                                       <  2:614077/78‑1 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTg                                                                       >  1:614077/1‑78 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAg                        <  2:335285/125‑1 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAg                        >  1:335285/1‑125 (MQ=255)
                                                                                                                    aaCAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGa               >  2:179145/1‑134 (MQ=255)
                                                                                                                          aggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATgg                <  2:492812/127‑1 (MQ=255)
                                                                                                                          aggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATgg                >  1:492812/1‑127 (MQ=255)
                                                                                                                          aggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAg         >  2:404748/1‑134 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         <  2:661129/69‑1 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         <  2:231895/69‑1 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         >  1:661129/1‑69 (MQ=255)
                                                                                                                           ggaaaaggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATc                                                                         >  1:231895/1‑69 (MQ=255)
                                                                                                                                 ggaaatTTTGGGCTGGGCTTGGGTAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTgg  <  2:748753/134‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                
TCGGCGAACCCGTGCCAGTGTGTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGG  >  AM260479/787747‑788007

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TCTTCGGCGAACCCGTGCCAGTGTGTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGGCAAGCGAGTGT  >  AM260479/787744‑788018
                                                                                                                                            |                                                                                                                                        
TCTTCGGCGAACCCGTGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCC                                                                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1108:11894:90197/1‑150 (MQ=60)
           cagGTGCCAGTGTGTTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTG                                                                                                                      <  GWNJ‑0478:712:GW2002102894th:2:1114:8963:87385/147‑1 (MQ=60)
                          TGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAA                                                                                                        <  GWNJ‑0478:712:GW2002102894th:2:1203:14743:42398/149‑1 (MQ=60)
                            caagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacaGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGT                                                                                                     <  GWNJ‑0478:712:GW2002102894th:2:2203:13131:20801/85‑1 (MQ=60)
                                 agtagctattgtatttcttttctttttttttttcaagcagaagacggcatacgagatgtagctccgtctcgtgggctcggagatgtgtataagagacagggaaatTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGC                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1215:2486:30476/45‑1 (MQ=60)
                                                CCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGG                                                                                 >  GWNJ‑0478:712:GW2002102894th:2:2205:2000:69490/1‑150 (MQ=60)
                                                                CTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTG                                                                 >  GWNJ‑0478:712:GW2002102894th:2:1216:2014:75569/1‑150 (MQ=60)
                                                                  GTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGC                                                                <  GWNJ‑0478:712:GW2002102894th:2:2202:6669:59523/149‑1 (MQ=60)
                                                                  GTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGC                                                                <  GWNJ‑0478:712:GW2002102894th:2:1208:14139:10189/149‑1 (MQ=60)
                                                                          GGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCAT                                                      <  GWNJ‑0478:712:GW2002102894th:2:2207:2431:89712/150‑1 (MQ=60)
                                                                             CATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTC                                                    >  GWNJ‑0478:712:GW2002102894th:2:2213:2886:11408/1‑150 (MQ=60)
                                                                                             gggctcggagatgtgtataagagacagACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATT                                     <  GWNJ‑0478:712:GW2002102894th:2:2114:18241:27582/122‑1 (MQ=60)
                                                                                                CCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCctgtctcttatacacatctccgagcccacg                                  >  GWNJ‑0478:712:GW2002102894th:2:2104:10827:52120/1‑119 (MQ=60)
                                                                                                CCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGC                                 <  GWNJ‑0478:712:GW2002102894th:2:1204:7393:93372/150‑1 (MQ=60)
                                                                                                    ATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCAT                             >  GWNJ‑0478:712:GW2002102894th:2:2114:9275:91560/1‑150 (MQ=60)
                                                                                                        GGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATctgtctcttatacacatctccgagcccacgagacggagctac                         >  GWNJ‑0478:712:GW2002102894th:2:2106:3080:81837/1‑108 (MQ=60)
                                                                                                        GGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGAT                         >  GWNJ‑0478:712:GW2002102894th:2:1207:13783:55806/1‑150 (MQ=60)
                                                                                                            AATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCAC                      <  GWNJ‑0478:712:GW2002102894th:2:2101:13898:12780/149‑1 (MQ=60)
                                                                                                             ATCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAG                    <  GWNJ‑0478:712:GW2002102894th:2:1202:11480:65555/150‑1 (MQ=60)
                                                                                                              TCCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGT                   >  GWNJ‑0478:712:GW2002102894th:2:1214:8132:73107/1‑150 (MQ=60)
                                                                                                               CCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATctgtctcttatacacatctccgagcccacgagacggagctacatctcgtat                  >  GWNJ‑0478:712:GW2002102894th:2:2111:11338:72556/1‑99 (MQ=60)
                                                                                                               CCTTTGATAACAAAAGGAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATctgtctcttatacacatctccgagcccacgagacggagctacatctcgtat                  >  GWNJ‑0478:712:GW2002102894th:2:1114:7435:66043/1‑99 (MQ=60)
                                                                                                                               GAAAAGGAAATTTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGGCAAGCGAGTGT  >  GWNJ‑0478:712:GW2002102894th:2:2105:7644:70994/1‑150 (MQ=60)
                                                                                                                                            |                                                                                                                                        
TCTTCGGCGAACCCGTGCCAGTGTGTGCCGCTTCGGTTTCAGACGTCCCCCTGTTCCGGGCTTCTGTTCGCGCGGACATATGTGCGACGAACGTGCCAGATATGGCTAATCCTTTGATAACAAAAGGAAAAGGAAATTTCGGCTGGCCTTGGGCAATTGCTGGCGATCAATAACGTTCCGCTTCTGATACATCTGGCGGGATTTCTATATTGCAGCGCAGCATTCGGCGAGCGGGTAATTAGGCGCATGGATCACAGTGCGTGGCAAGCGAGTGT  >  AM260479/787744‑788018

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