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
A4 F4 I1 R1 17 37.4 1246952 98.2% 1224506 142.8

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,992,784 C→A Y544* (TAC→TAA cyaA → adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,992,7840CA100.0% 120.3 / NA 35Y544* (TAC→TAAcyaAadenylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (23/12);  total (23/12)

CTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAG  >  NC_000913/3992658‑3992926
                                                                                                                              |                                                                                                                                              
cTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATc                                                                                                                          >  4:23513/1‑149 (MQ=255)
 tCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCt                                                                                                                         >  6:130530/1‑149 (MQ=255)
        cccGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGAt                                                                                                                   >  1:166296/1‑148 (MQ=255)
         ccGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTa                                                                                                                 >  5:103613/1‑149 (MQ=255)
            tttGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCg                                                                                                              <  5:130530/149‑1 (MQ=255)
                tataTTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCgtgg                                                                                                                                       >  8:70595/1‑119 (MQ=255)
                tataTTAAAGCTAACCTAATTGACGAATTCCCTAACTTGCAGGAGATGGTCGCAGAACTGTCGCACCATTTCCCGCTGCTCTTACCTGCACCGCCACAGAAGGCTCTCTAAAGCCCGTCTGAGATCCGACATCTGGCGATTATCGTTaa                                                                                                          >  4:11852/1‑149 (MQ=255)
                 atatTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGccc                                                                                                                                           >  4:34305/1‑115 (MQ=255)
                         ggTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAata                                                                                                 >  2:19652/1‑149 (MQ=255)
                              cGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGAcc                                                                                            <  2:166296/149‑1 (MQ=255)
                                        gATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCCTTTCCCGCTGCGCTTCCCTGCACCGACACCGAAGGCGCTCTAACGCCCGTGTGAGATCAGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAgcggc                                                                                  >  7:47329/1‑149 (MQ=255)
                                                       cAGGAGATGGTCGCCGACGTGTCGCAACATTTCCCTCTGCGCTTACCTGCAACGACACCGAAGTAGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACACGAAAGCGGAGTTCCGCAATCAggt                                                                   >  6:77835/1‑149 (MQ=255)
                                                                gTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTCGaa                                                                                                   <  6:93160/108‑1 (MQ=255)
                                                                     cGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGCCCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAtt                                                     >  4:84817/1‑149 (MQ=255)
                                                                      gACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGAttt                                                    <  6:56586/149‑1 (MQ=255)
                                                                              gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCt                                                                                                       <  3:68529/90‑1 (MQ=255)
                                                                              gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCt                                                                                                       >  4:68529/1‑90 (MQ=255)
                                                                              gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGc                                            >  4:102035/1‑149 (MQ=255)
                                                                              gCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATAATCGTTAACCTGGAATATGACCCGACAGCGGCGGTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGc                                            >  3:78983/1‑149 (MQ=255)
                                                                               cACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCt                                           >  7:96590/1‑149 (MQ=255)
                                                                                aCCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTg                                          >  8:23301/1‑149 (MQ=255)
                                                                                  cATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGa                                        >  5:2173/1‑149 (MQ=255)
                                                                                   aTTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGAt                                       <  6:103613/149‑1 (MQ=255)
                                                                                    tttCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATg                                      >  2:8856/1‑149 (MQ=255)
                                                                                         cGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCtt                                  >  1:119375/1‑148 (MQ=255)
                                                                                          gCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCa                                <  1:19652/149‑1 (MQ=255)
                                                                                               gcTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCttt                           <  2:119375/149‑1 (MQ=255)
                                                                                                   aCCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGc                        >  2:163720/1‑148 (MQ=255)
                                                                                                    ccTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCga                      <  3:102035/149‑1 (MQ=255)
                                                                                                    ccTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGAGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCga                      <  5:77835/149‑1 (MQ=255)
                                                                                                           cgacaccgaAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCa                >  1:4515/1‑148 (MQ=255)
                                                                                                                 cgaAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGc         >  1:159091/1‑149 (MQ=255)
                                                                                                                      gcgcTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTggt    <  1:163720/149‑1 (MQ=255)
                                                                                                                          tctAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTggtag  >  6:77821/1‑147 (MQ=255)
                                                                                                                          tctAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTggtag  <  5:77821/147‑1 (MQ=255)
                                                                                                                              |                                                                                                                                              
CTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAG  >  NC_000913/3992658‑3992926

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGG  >  NC_000913/3992641‑3992927
                                                                                                                                               |                                                                                                                                               
gactggagttcagacgttgtgctcttccgatctTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCG                                                                                                                                           <  NB501016:228:HFK2FAFXY:2:11206:14560:19278/117‑1 (MQ=60)
     ggagttcagacgtggctcttccgatcTGTATATTCAAGGTAACGGCATTGGCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTG                                                                                                                                       <  NB501016:228:HFK2FAFXY:4:11612:23651:18687/123‑1 (MQ=60)
                        ACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATT                                                                                                                   >  NB501016:228:HFK2FAFXY:1:21307:22239:12882/1‑150 (MQ=60)
                         CCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCTTGTGAGATCCGCCATCTGGCGATTAT                                                                                                                 >  NB501016:228:HFK2FAFXY:3:21410:23533:17098/1‑151 (MQ=60)
                            GTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGT                                                                                                              <  NB501016:228:HFK2FAFXY:3:21511:6056:11502/151‑1 (MQ=60)
                                                        CGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCCTTTCCCGCTGCGCTTCCCTGCACCGACACCGAAGGCGCTCTAACGCCCGTGTGAGATCAGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCG                                                                                  >  NB501016:228:HFK2FAFXY:4:11512:2945:9263/1‑151 (MQ=60)
                                                                                GGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATAgatcggaagagcacacgtctgactccagtcacctgaagct                                                          >  NB501016:228:HFK2FAFXY:3:21405:12076:19940/1‑111 (MQ=60)
                                                                                              CGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTG                                                                                                       <  NB501016:228:HFK2FAFXY:2:11311:14610:5636/92‑1 (MQ=60)
                                                                                              CGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATAATCGTTAACCTGGAATATGACCCGACAGCGGCGGTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCT                                            >  NB501016:228:HFK2FAFXY:2:21104:16797:8475/1‑151 (MQ=60)
                                                                                               GCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTG                                           >  NB501016:228:HFK2FAFXY:4:21501:16055:6166/1‑151 (MQ=60)
                                                                                                  CCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGAT                                        >  NB501016:228:HFK2FAFXY:3:11401:11065:19193/1‑151 (MQ=60)
                                                                                                         CCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTC                                  >  NB501016:228:HFK2FAFXY:1:21112:12709:10127/1‑150 (MQ=60)
                                                                                                          CGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAG                                <  NB501016:228:HFK2FAFXY:1:11108:9479:11046/151‑1 (MQ=60)
                                                                                                                    ACCTGCAGCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAG                      <  NB501016:228:HFK2FAFXY:2:21204:11601:19914/151‑1 (MQ=60)
                                                                                                                    ACCTGCACCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGAGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAG                      <  NB501016:228:HFK2FAFXY:3:11611:2889:1528/151‑1 (MQ=60)
                                                                                                                           CCGACACCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAA                >  NB501016:228:HFK2FAFXY:1:11102:18070:15353/1‑150 (MQ=60)
                                                                                                                                 CCGAAGGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCC         >  NB501016:228:HFK2FAFXY:1:21304:20416:13772/1‑151 (MQ=60)
                                                                                                                                      GGCGCTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTA    <  NB501016:228:HFK2FAFXY:1:21306:1336:11881/151‑1 (MQ=60)
                                                                                                                                          CTCTAAAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGG  <  NB501016:228:HFK2FAFXY:3:11611:17083:1414/149‑1 (MQ=60)
                                                                                                                                               |                                                                                                                                               
TTTAACGGCCTGCTGACCTCGCGCACCCGTTTGTATATTAAAGGTAACGGCATTGTCGATTTGCCTAAGTTGCAGGAGATGGTCGCCGACGTGTCGCACCATTTCCCGCTGCGCTTACCTGCACCGACACCGAAGGCGCTCTACAGCCCGTGTGAGATCCGCCATCTGGCGATTATCGTTAACCTGGAATATGACCCGACAGCGGCGTTCCGCAATCAGGTGGTGCATTTCGATTTCCGTAAGCTGGATGTCTTCAGCTTTGGCGAGAATCAAAATTGCCTGGTAGG  >  NC_000913/3992641‑3992927

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