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
A0 F0 I1 R1 146 27.9 2327082 89.7% 2087392 107.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA AM260479 3,529,526 +G intergenic (‑40/‑264) secA ← / → h16_A3265 preprotein translocase subunit SecA/Hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AM2604793,529,5241.G100.0% 143.9 / NA 38intergenic (‑38/‑266)secA/h16_A3265preprotein translocase subunit SecA/Hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (18/20);  total (18/20)

GCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTG  >  AM260479/3529398‑3529648
                                                                                                                               |                                                                                                                             
gcATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAgggatgc                                                                                                                        >  1:375522/1‑129 (MQ=255)
gcATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAgggatgc                                                                                                                        >  1:754138/1‑129 (MQ=255)
            gcCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCcgg                                                                                                             <  1:955528/134‑1 (MQ=255)
               aCCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCcggcgg                                                                                                          >  1:1042245/1‑134 (MQ=255)
                    gcggcgATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGCCTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCATCCCAGGGATGCGGGCTGCTCCGGCGGCCgg                                                                                                      <  1:838094/133‑1 (MQ=255)
                    gcggcgATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCggg                                                                                                     >  2:804055/1‑134 (MQ=255)
                    gcggcgATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCgg                                                                                                      >  2:838094/1‑133 (MQ=255)
                      ggcgATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCg                                                                                                   <  1:791362/134‑1 (MQ=255)
                      ggcgATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCg                                                                                                   <  2:754138/134‑1 (MQ=255)
                                 ttGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGa                                                                                          >  2:331320/1‑132 (MQ=255)
                                            cTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             <  1:662531/134‑1 (MQ=255)
                                            cTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             <  1:582371/134‑1 (MQ=255)
                                               aTTGCGGCTTCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaa                                                                            <  2:1042245/132‑1 (MQ=255)
                                               aTTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaaa                                                                           <  2:375522/133‑1 (MQ=255)
                                                   cGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaaa                                                                           >  2:770276/1‑129 (MQ=255)
                                                   cGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaaa                                                                           <  1:770276/129‑1 (MQ=255)
                                                             aGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGcc                                                                                                        >  1:637380/1‑90 (MQ=255)
                                                             aGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGcc                                                                                                        <  2:637380/90‑1 (MQ=255)
                                                                       gagaAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGt                                                  >  2:211360/1‑134 (MQ=255)
                                                                                 gTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCgg                                                     <  2:704089/121‑1 (MQ=255)
                                                                                 gTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCgg                                                     >  1:704089/1‑121 (MQ=255)
                                                                                     tcatcGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCcgg                                                                                                             <  2:379242/61‑1 (MQ=255)
                                                                                     tcatcGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCcgg                                                                                                             >  1:379242/1‑61 (MQ=255)
                                                                                          gATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCt                                <  1:331320/133‑1 (MQ=255)
                                                                                                     ccTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             >  1:191192/1‑77 (MQ=255)
                                                                                                     ccTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             <  2:191192/77‑1 (MQ=255)
                                                                                                      cTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                     <  2:866187/116‑1 (MQ=255)
                                                                                                      cTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                     >  1:866187/1‑116 (MQ=255)
                                                                                                         cgcgGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAg                >  1:1097360/1‑134 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             <  2:286746/69‑1 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             <  2:549891/69‑1 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             >  1:549891/1‑69 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaa                                                                             >  1:286746/1‑69 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                     <  2:293347/109‑1 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAAcg                                     >  1:293347/1‑109 (MQ=255)
                                                                                                             gAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAggg              >  1:870022/1‑132 (MQ=255)
                                                                                                                 cGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGggcaggca        <  2:1069804/134‑1 (MQ=255)
                                                                                                                       gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaa                                                                            >  2:681712/1‑60 (MQ=255)
                                                                                                                       gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCaaa                                                                            <  1:681712/60‑1 (MQ=255)
                                                                                                                       gcACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTg  >  2:1080392/1‑134 (MQ=255)
                                                                                                                               |                                                                                                                             
GCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTG  >  AM260479/3529398‑3529648

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGAACTTCGGCTCCAGCGCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGT  >  AM260479/3529381‑3529669
                                                                                                                                                  |                                                                                                                                                
CGAACTTCGGCTCCAGCGCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGga                                                                                                                                                <  GWNJ‑0478:712:GW2002102894th:2:1215:15494:46041/148‑3 (MQ=60)
  AACTTCGGCTCCAGCGCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCG                                                                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:2110:7949:14158/1‑150 (MQ=60)
  AACTTCGGCTCCAGCGCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCG                                                                                                                                             >  GWNJ‑0478:712:GW2002102894th:2:1207:2116:41223/1‑150 (MQ=60)
                 GCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGC                                                                                                                              >  GWNJ‑0478:712:GW2002102894th:2:2210:17490:6221/1‑150 (MQ=60)
                            CGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAA                                                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2209:7155:25859/150‑1 (MQ=60)
                                   GGGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGCCTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCATCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAG                                                                                                            <  GWNJ‑0478:712:GW2002102894th:2:2114:1980:41811/150‑1 (MQ=60)
                                      CGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATT                                                                                                         <  GWNJ‑0478:712:GW2002102894th:2:2112:18408:7278/150‑1 (MQ=60)
                                                            GCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATC                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:1216:13615:97001/150‑1 (MQ=60)
                                                            GCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATC                                                                                   <  GWNJ‑0478:712:GW2002102894th:2:2104:4382:98348/150‑1 (MQ=60)
                                                             gaagacggcatacgagattagcgagtgtctcgtgggctcggagatgtgtataagagacagGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCA                                                                                  <  GWNJ‑0478:712:GW2002102894th:2:2106:20110:10267/90‑1 (MQ=60)
                                                               CATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATctgtctcttatacacatctccgagcccacg                                                                                >  GWNJ‑0478:712:GW2002102894th:2:2105:6997:46042/1‑120 (MQ=60)
                                                               CATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCAT                                                                                 <  GWNJ‑0478:712:GW2002102894th:2:2110:15217:99548/149‑1 (MQ=60)
                                                                                   TTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGA                                                             >  GWNJ‑0478:712:GW2002102894th:2:2108:4606:42837/1‑149 (MQ=60)
                                                                                       TGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAActgtctcttatacacatctccgagcccacgagacactcgctaatctcgtatgccgtctt                                                        >  GWNJ‑0478:712:GW2002102894th:2:1207:15077:61598/1‑91 (MQ=60)
                                                                                                       CATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCctgtctcttatacacatctccgagcccacgagacactcgcta                                         >  GWNJ‑0478:712:GW2002102894th:2:1113:10036:13509/1‑107 (MQ=60)
                                                                                                        ATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCtgt                                       >  GWNJ‑0478:712:GW2002102894th:2:2116:4815:80763/1‑147 (MQ=60)
                                                                                                          CGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAA                                      <  GWNJ‑0478:712:GW2002102894th:2:1202:11854:78814/149‑1 (MQ=60)
                                                                                                           GATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGG                                    >  GWNJ‑0478:712:GW2002102894th:2:2213:16402:87131/1‑150 (MQ=60)
                                                                                                               GCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtat                                >  GWNJ‑0478:712:GW2002102894th:2:1214:7804:40091/1‑99 (MQ=60)
                                                                                                               GCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCctgtctcttatacacatctccgagcccacgagacactcgctaatctcgtat                                >  GWNJ‑0478:712:GW2002102894th:2:1116:11292:99971/1‑99 (MQ=60)
                                                                                                               GCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCtgtctcttat                                >  GWNJ‑0478:712:GW2002102894th:2:1201:9992:43146/1‑140 (MQ=60)
                                                                                                               GCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGC                                  >  GWNJ‑0478:712:GW2002102894th:2:2201:21148:2067/1‑148 (MQ=60)
                                                                                                                                             CCAGGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGT  >  GWNJ‑0478:712:GW2002102894th:2:1115:20099:31080/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                
CGAACTTCGGCTCCAGCGCATTGATCTGCGCCACCGTGCGGCGATATTGTTTGATCAGCCGCTCATTGCGGCTGCCGAAGACTTTCTTGAGAAGGCCGGTGATCATCGATTGCTTTCACCTGCGCGGAACCGGGCGGCACCCCAGGATGCGGGCTGCTCCGGCGGCCGGGCGAGGAATGGATAGATTGGCGCAAAACGCCGAGTTTATCATGTGCCGGGGTCACGGCGGAACGCGCCTGGCCCCGTCCCGCAAGGGCAGGCAAGCCTGCCGGATGTCGGCCACGCGTGT  >  AM260479/3529381‑3529669

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