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
A2 F5 I184 R1 278 29.6 1166432 97.3% 1134938 124.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 737,899 A→G N38N (AAT→AAC dtpD ← dipeptide permease DtpD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273737,8990AG100.0% 134.5 / NA 37N38N (AAT→AACdtpDdipeptide permease DtpD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/20);  total (17/20)

CCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAAT  >  NZ_CP009273/737766‑738015
                                                                                                                                     |                                                                                                                    
cccccAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTtggttggt                                                                                                                 <  2:172621/139‑1 (MQ=255)
cccccAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTtggttggt                                                                                                                 <  1:572765/139‑1 (MQ=255)
   ccAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTtggttggt                                                                                                                 >  1:336104/1‑136 (MQ=255)
            ccGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAAt                                                                                                     <  2:379421/139‑1 (MQ=255)
               ccATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcag                                                                                                  <  1:156955/139‑1 (MQ=255)
               ccATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcag                                                                                                  >  1:507650/1‑139 (MQ=255)
                 aTGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAACTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATcag                                                                                                  >  2:371717/1‑137 (MQ=255)
                          ccGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCa                                                                                       <  2:387704/139‑1 (MQ=255)
                           cGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCAt                                                                                      <  2:580127/139‑1 (MQ=255)
                            gagaACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATg                                                                                     >  2:448866/1‑139 (MQ=255)
                                    tATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTaa                                                                               >  2:562740/1‑137 (MQ=255)
                                                          gATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGc                                                       <  1:178105/139‑1 (MQ=255)
                                                                 gTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCaa                                                <  2:50842/139‑1 (MQ=255)
                                                                    aCATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                              <  2:507650/138‑1 (MQ=255)
                                                                     cATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc                                            <  2:336104/139‑1 (MQ=255)
                                                                              gCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                    >  2:184475/1‑138 (MQ=255)
                                                                              gCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                    >  2:469007/1‑138 (MQ=255)
                                                                               cGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTaaa                                    >  2:377394/1‑137 (MQ=255)
                                                                                gAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAg                                 >  2:15163/1‑139 (MQ=255)
                                                                                 aaCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                                   >  1:433924/1‑120 (MQ=255)
                                                                                 aaCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAgcg                                                   <  2:433924/120‑1 (MQ=255)
                                                                                 aaCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc                                            <  1:359302/127‑1 (MQ=255)
                                                                                 aaCAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAc                                            >  2:359302/1‑127 (MQ=255)
                                                                                    aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                             <  1:184475/139‑1 (MQ=255)
                                                                                    aGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAgcgcg                             <  1:562740/139‑1 (MQ=255)
                                                                                      tAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGcgcgcg                           >  1:146346/1‑139 (MQ=255)
                                                                                               aaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATg                  <  1:552895/139‑1 (MQ=255)
                                                                                               aaTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATg                  <  1:273288/139‑1 (MQ=255)
                                                                                                aTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                     >  1:290059/1‑103 (MQ=255)
                                                                                                aTAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAg                                                     <  2:290059/103‑1 (MQ=255)
                                                                                                    cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt             >  1:75401/1‑139 (MQ=255)
                                                                                                    cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCgtgt             <  2:146346/139‑1 (MQ=255)
                                                                                                    cTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCAGCTGTGATGCgtgt             <  2:276524/139‑1 (MQ=255)
                                                                                                         aGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                              <  2:309851/101‑1 (MQ=255)
                                                                                                         aGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                              >  1:309851/1‑101 (MQ=255)
                                                                                                         aGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt        >  2:355026/1‑139 (MQ=255)
                                                                                                         aGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTAtt        >  1:507800/1‑139 (MQ=255)
                                                                                                                aGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCataat  <  1:377394/138‑1 (MQ=255)
                                                                                                                          tATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                              >  2:198064/1‑84 (MQ=255)
                                                                                                                          tATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACg                                              <  1:198064/84‑1 (MQ=255)
                                                                                                                                     |                                                                                                                    
CCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAAT  >  NZ_CP009273/737766‑738015

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAACAACGCCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAG  >  NZ_CP009273/737758‑738048
                                                                                                                                             |                                                                                                                                                     
CAACAACGCCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAG                                                                                                                                               <  SRR3722139.202430/150‑1 (MQ=60)
CAACAACGCCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTG                                                                                                                                                 >  SRR3722139.340209/1‑148 (MQ=60)
       GCCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAA                                                                                                                                       <  SRR3722139.579313/151‑1 (MQ=60)
            CAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGC                                                                                                                                  >  SRR3722139.513498/1‑151 (MQ=60)
                      GCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCA                                                                                                                        <  SRR3722139.158873/151‑1 (MQ=60)
                                                                 GGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGAC                                                                             <  SRR3722139.180311/151‑1 (MQ=60)
                                                                            AGATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAG                                                                  <  SRR3722139.363631/151‑1 (MQ=60)
                                                                              ATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGctgtctctt                                                                >  SRR3722139.439024/1‑142 (MQ=60)
                                                                               agtacggtctcgtgggctcggagatgtgtataagagacagGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAAT                                                                    <  SRR3722139.200505/106‑1 (MQ=60)
                                                                                   CCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGG                                                           >  SRR3722139.148129/1‑151 (MQ=60)
                                                                                           CAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATG                                                   <  SRR3722139.186751/151‑1 (MQ=60)
                                                                                           CAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATG                                                   <  SRR3722139.569195/151‑1 (MQ=60)
                                                                                             GTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATctgtctcttatacacatctccgagcc                                                 >  SRR3722139.293710/1‑125 (MQ=60)
                                                                                                 GCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTT                                             >  SRR3722139.76272/1‑151 (MQ=60)
                                                                                                      CAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTC                                        >  SRR3722139.513652/1‑151 (MQ=60)
                                                                                                      AAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATctgtctcttatacacatctccgagccc                                         >  SRR3722139.313686/1‑123 (MQ=60)
                                                                                                      AAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTC                                        <  SRR3722139.276758/151‑1 (MQ=60)
                                                                                                      AAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGGGTTTATTC                                        <  SRR3722139.559249/151‑1 (MQ=60)
                                                                                                                       GAGTATCGTTGTATTTTAGTTGGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAAT                        <  SRR3722139.381905/150‑1 (MQ=60)
                                                                                                                                            GGTTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAG  <  SRR3722139.287598/151‑1 (MQ=60)
                                                                                                                                             |                                                                                                                                                     
CAACAACGCCCCCAGCATCACCGCCATGCGATTGCCGAGAACTTTATCCGCCAAAAAGCCACCGAGGATTGGCGTGACATACACCAGCGAACAGTAGGCGCTAAATAACTCGTAGGCGTGAGTATCGTTGTATTTTAGTTGATTGGTGAGATAGAGAATCAGCAGGGCACGCATGCCGTAAAAGCTAAAATATTCCCAGATTTGCAGCGCAACGACATAGTAAATAGCGCGCGGCTGTGATGCGTGTTTATTCATAATAATTCCACAATTGACGGCAACGTCACGTAATAG  >  NZ_CP009273/737758‑738048

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