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 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,958,260 +C intergenic (‑829/‑32) ppiC ← / → rep peptidyl‑prolyl cis‑trans isomerase C/rep helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,958,2581.C79.4% 150.0 / 40.6 64intergenic (‑827/‑34)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base . (10/3);  new base C (22/28);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 5.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-02

Only 100 of 111 total aligned reads displayed.

CCGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑T‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCC  >  NC_000913/3958166‑3958341
                                                                                                 |                                                                                             
 cGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑c                                                                                             >  2:4081931/1‑93 (MQ=255)
 cGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑c                                                                                             <  2:2098638/93‑1 (MQ=255)
  gATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑ttc                                                                                              >  1:976842/1‑91 (MQ=255)
  gATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑cc                                                                                            >  1:3328748/1‑93 (MQ=255)
  gATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGATAGCGGCTTACACGCGGTCACCTTCAAATGCGATTCTGCTACAA‑‑t‑‑cc                                                                                            >  2:754235/1‑93 (MQ=255)
   aTTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑ccc                                                                                           >  2:3173151/1‑92 (MQ=255)
   aTTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑ccc                                                                                           >  1:1202453/1‑92 (MQ=255)
    tttCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑cccc                                                                                          >  1:480089/1‑91 (MQ=255)
    tttCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑cccc                                                                                          >  2:32204/1‑91 (MQ=255)
    tttCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑cccc                                                                                          >  2:673689/1‑91 (MQ=255)
      tCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑t‑‑cccacc                                                                                        >  2:1636796/1‑89 (MQ=255)
        ggTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑cc‑‑c‑‑ccccc                                                                                  >  1:4401209/1‑93 (MQ=255)
        ggTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑cc‑‑c‑‑ccccc                                                                                  <  2:4501298/93‑1 (MQ=255)
         gTCTTCTCTCTCTGCTTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGACTCTGCTACAC‑‑A‑‑Cc‑‑c‑‑cccccc                                                                                 >  2:3026830/1‑93 (MQ=255)
         gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑Cc‑‑c‑‑cccccc                                                                                 >  2:3936964/1‑93 (MQ=255)
         gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑Cc‑‑c‑‑cccccc                                                                                 >  1:3723554/1‑93 (MQ=255)
         gTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGCTAGCGGCATACACGCGGTCACATTCAAATGCGACTCTGCTACAA‑‑T‑‑CC‑‑T‑‑cccccc                                                                                 >  2:3791412/1‑93 (MQ=255)
          tcttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑C‑‑caccccc                                                                                >  2:3044620/1‑92 (MQ=255)
            ttctCTCTCTGATTTAAACATCCGCAGCCAACTTGTTTGCGGCTTACGCGCGGTCACATTCAAAGGCGATTCTGCTACAA‑‑T‑‑Cc‑‑c‑‑cccccccgt                                                                              <  2:4178842/93‑4 (MQ=255)
            ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑Cc‑‑c‑‑ccccccccg                                                                              >  2:1502607/1‑90 (MQ=255)
            ttctCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑Cc‑‑c‑‑ccccccccc                                                                              >  1:129859/1‑90 (MQ=255)
              ctctctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑T‑‑cccccccgta                                                                             >  2:2012588/1‑89 (MQ=255)
                 tctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGaa                                                                          <  1:585573/93‑1 (MQ=255)
                  ctctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAg                                                                        <  2:3807482/93‑1 (MQ=255)
                   tgtGATTTAAACCGCCGCGGCCAACCGGGTAGCGGCTTTCCCGGCGTGAAGTTCAAATCCATTTCCGCTACCAC‑‑C‑‑CC‑‑C‑‑CCCCCCGTTCGAAGa                                                                       <  2:2778192/91‑1 (MQ=255)
                   tctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGa                                                                       <  2:662099/93‑1 (MQ=255)
                   tctGATTTAAACATCCGCAGCCAACCGGGTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAg                                                                        <  1:846323/93‑1 (MQ=255)
                    ctGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGAt                                                                      <  1:3205071/93‑1 (MQ=255)
                     tGATTTAAACCTCCCGAGCCAACCGGTTAGCGGCTTACCCGCGGTCCCATTCAAATGCGATTCTGCCACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGAt                                                                      <  2:2763350/93‑1 (MQ=255)
                     tGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCCCATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGAt                                                                      <  2:3242446/93‑1 (MQ=255)
                     tGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGAt                                                                      <  1:512836/93‑1 (MQ=255)
                     tGAATTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCTGTTACATACAAATGCGAATCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGAt                                                                      <  2:3915100/93‑1 (MQ=255)
                        tttAAACCACCGGAACCAAGCGGTTTGGAGGTTCGACGCGGTGACATTCAAATTGGATTTGGGTGCAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGa                                                                   <  2:720491/93‑1 (MQ=255)
                        tttAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAACCC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGa                                                                   <  1:4623660/93‑1 (MQ=255)
                        tttAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAg                                                                  >  1:3104589/1‑93 (MQ=255)
                        attAAACATCCGCAGCCAACCGGTTAGCCGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGa                                                                   <  2:4081986/92‑1 (MQ=255)
                         ttAAACATCCGCAGCCAACCGGTTAGCGGCTTACTCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCCAAGTTTGAg                                                                  >  2:2992938/1‑93 (MQ=255)
                          tAAACATCCGCAGCCAACCGGTTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGc                                                                 <  2:906268/93‑1 (MQ=255)
                          tAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGc                                                                 <  2:4301858/93‑1 (MQ=255)
                          tAAACACCCCCAGCCAACCCGGTAGCCGCCTACCCGCCGTCACATTCAAATGCGATTCTGCTACAACCC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGc                                                                 <  2:4129563/93‑1 (MQ=255)
                              cATCCGCCGCCAACCGGGTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATa                                                            <  1:990902/93‑1 (MQ=255)
                               aTCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CCCAC‑‑CCCCCCGTTCGAAGATTGAGCAATac                                                           <  1:4232044/93‑1 (MQ=255)
                                accGCAGCCAACCGGTTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATaca                                                          <  1:1514394/92‑1 (MQ=255)
                                 ccGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACAcc                                                        <  2:3807158/93‑1 (MQ=255)
                                   gCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCCTTCCAAGATTAAGCAACACACCt                                                       >  2:2102449/1‑93 (MQ=255)
                                    cAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTAt                                                    <  2:1813702/93‑1 (MQ=255)
                                    cAGCCAACCCGTTAGCGTCTTACCCGCGGTCCCATTCAAAAGCGAGTCTGCCAAAAGCC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTa                                                     <  2:811686/93‑1 (MQ=255)
                                     cgCCAACCGGTTAGCGGCTTACACGCGGTCGAATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATg                                                   <  2:298492/92‑1 (MQ=255)
                                     cgCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTAt                                                    <  1:3277072/92‑1 (MQ=255)
                                      gCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTAtc                                                   >  1:584662/1‑92 (MQ=255)
                                       ccAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGc                                                  <  1:4074116/93‑1 (MQ=255)
                                        cAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATAGACCTATgggc                                                >  2:1904487/1‑90 (MQ=255)
                                          aCCGGTTAGGGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑Tc                                             <  1:4430935/93‑1 (MQ=255)
                                             cgTTAGCGGCTTTCCAGCGGGCCCATTCAAATGGGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTaa                                          <  1:4499007/92‑1 (MQ=255)
                                              gTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAAc                                        <  1:686124/93‑1 (MQ=255)
                                               ttAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTCCCAAGATTGACCAATACACTTATGCG‑‑TCTTAAc                                        >  2:3905541/1‑93 (MQ=255)
                                               ttAGCGGCTTACACGCGGTCACATTCAAATGCCATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACATCTATGCG‑‑TCTAAAc                                        <  2:2863190/93‑1 (MQ=255)
                                                   cGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCg                                    <  1:2467014/93‑1 (MQ=255)
                                                    ggCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCAATACG‑‑TCTAAACCCCgc                                   >  1:2390259/1‑92 (MQ=255)
                                                      cTTACACGCGGTCACATTCAAATGCGATTATGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATAAACCTATGGG‑‑CCTAAACCCCGGCca                                >  1:4102310/1‑93 (MQ=255)
                                                           acGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTCCGAAGAATGAGAAATACCACTATGCG‑‑TCTAAACCCCGGCcaaaca                            >  2:1441080/1‑90 (MQ=255)
                                                            cgcgGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCcaacaac                           <  1:4576805/93‑1 (MQ=255)
                                                            cgcgGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGCTCGAAGATTAAGCAAAACACCTATTCG‑‑TCTAAACCCCGCGcaacaaa                           >  1:192847/1‑92 (MQ=255)
                                                            cgaggTCTCATTCAAAATCGATTCGGACACAACCC‑‑CC‑‑C‑‑CCCCCCGTGCGAAGATTGAGCAATACACCTATGCG‑‑TATAAACCCCGGCcaacaac                           <  2:3609865/90‑1 (MQ=255)
                                                             gcgGTCACATTCAAATGCGATTTCGCTAAAACC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCcaacaacaa                         <  1:2117532/93‑1 (MQ=255)
                                                              cgGTCACATTCAAATTGCATTCTGCCACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCcaacaacaa                         <  2:4087260/93‑1 (MQ=255)
                                                               ggTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑CCCCCCCCCCCTCCAAGCTTGGGAAAAACACCTATGCGCCTCAAAACCCCCCacaacaaca                          >  1:2647443/1‑93 (MQ=255)
                                                               ggTCACATTCAAAAGCGATTCTGCTACAACCC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAg                        <  1:2318301/93‑1 (MQ=255)
                                                                gTCACATTCAAATGCGATTCTGCTACCACCC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCAGCCAACAACAAGc                       <  2:1030714/93‑1 (MQ=255)
                                                                gTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGc                       <  1:824143/93‑1 (MQ=255)
                                                                gTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGCTCCAAGATTGTGCAATACACCTATTGC‑‑TCCAAACCCCCCCccacaacaacc                       >  2:2773053/1‑91 (MQ=255)
                                                                gTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑CCTAAACCCCGGCcaacaacaaaat                      >  2:2989211/1‑90 (MQ=255)
                                                                gTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑CCCCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCt                      >  2:3472489/1‑93 (MQ=255)
                                                                gTCACATTCAAATGCCATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGc                       <  1:3341274/93‑1 (MQ=255)
                                                                cttacaTTCAAATTCAATCCCGCTTCAAA‑‑C‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCt                      <  1:631016/90‑1 (MQ=255)
                                                                 tCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTTAGCAATACACCTATGCG‑‑TCTAAACCCCCGCCAACCACAAGCtt                     >  2:3022422/1‑92 (MQ=255)
                                                                 tCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTg                     >  1:768291/1‑93 (MQ=255)
                                                                   cctttCAAATGCGTTTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGt                    <  2:910597/90‑1 (MQ=255)
                                                                   acaTTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑CCCCCCCCCGTTCGAAAATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGCCCAACAACAAACTgc                    >  2:695848/1‑92 (MQ=255)
                                                                   acaTTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑CCCCCCCCCCTTGGAAGATTGACCAATAACCTAATCCT‑‑CCTAAACCCCGCCCAACAAAAACCtcc                    >  2:731391/1‑91 (MQ=255)
                                                                    caTTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑CCCCCCCCCGTTCGAAGATTGGCAACTACACCCATGGG‑‑TCTAAAACCCGGCCAAAAACAAGCttca                   >  1:4370240/1‑90 (MQ=255)
                                                                     aTTCAACTGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGGAATACACCTATGCG‑‑TGTAAACCCCGGCCAACAACAAGCTGTCGa                 >  1:41481/1‑93 (MQ=255)
                                                                     aTTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCCTTCGACGATTGAGGAATATACCCACGCG‑‑TCTAAACCCCGGCCAACAAAAAACTGGCg                  >  1:2757985/1‑93 (MQ=255)
                                                                      ttCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTTCAGGATTTGGAAATAAACCTTTTCG‑‑GCTAAAACCCCCCCCACAACAACCTgctga                 >  2:2929074/1‑89 (MQ=255)
                                                                      ttCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAGGATTGACCAAAAAGCCTATGGT‑‑GCTAAACCCCGGCAAACAACAAACTATcta                 >  2:2951286/1‑91 (MQ=255)
                                                                      ttCAAATGCGATTCTGCTACAA‑‑TACCC‑‑C‑‑CCCCCCCCTCAAAGATATACCAATCAACCTATGCG‑‑CCTAAACCCCGACCAACACCAAGCTgctga                 >  2:1051803/1‑89 (MQ=255)
                                                                       taaaaTGCCGTTCCGCTACAA‑‑C‑‑CC‑‑CCCCCCCCCGTTCGAAGATTGAGCAATCCCCCTATGCG‑‑TCTAACCCCCGGCCAACAACAAGCTGTCGaa                <  2:2961613/91‑1 (MQ=255)
                                                                       tCAAATGCGATTTTGCTACAAACC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGaa                <  2:42708/93‑1 (MQ=255)
                                                                       tCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCATTACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGaa                <  2:1113773/93‑1 (MQ=255)
                                                                        taaaTGCGAATCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAAt               <  1:4533842/92‑1 (MQ=255)
                                                                         aaaTGCGATTCTGCTACAAGC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTc             <  1:2604362/93‑1 (MQ=255)
                                                                          aaggCGATCCTCCTACAC‑‑C‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGt           <  1:3689792/90‑1 (MQ=255)
                                                                          aaTGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCCGTCTAAGGATTTAGAGAAAACCCTATTGG‑‑CCTAAACCCCGGCCAGCAACACGATGTTGaagtt             >  2:4332518/1‑90 (MQ=255)
                                                                          aaTGCGATTCTGCTACAA‑‑T‑‑CC‑‑C‑‑CCCCCCGCTCGAAGATTTAGCAATACACCCATGCG‑‑TCTAAAAACCGGCCAACAAAAAGCTGTCGAATTCGt           >  2:3138943/1‑93 (MQ=255)
                                                                             gCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGtt          <  1:2173362/93‑1 (MQ=255)
                                                                              cGATTCTGCTACAAACC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTa         <  1:3428831/93‑1 (MQ=255)
                                                                                aTTCTGCCACAACC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCg      <  2:2355036/93‑1 (MQ=255)
                                                                                aTTCTGCCACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAcc       <  1:2266535/93‑1 (MQ=255)
                                                                                 ttCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACCCCTATGCG‑‑TCTAAACCCCGGCCAAAAACAAGCTGTCGAATTCGTTACCgg     >  2:1869179/1‑93 (MQ=255)
                                                                                    tGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGAGTGGGCCATACACCTATGCG‑‑TCTAACACCCGGCCCACAACAACCTGTCGAATTCGTTACCCGcac  >  2:3765369/1‑91 (MQ=255)
                                                                                                 |                                                                                             
CCGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑T‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCC  >  NC_000913/3958166‑3958341

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑TCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGCT  >  NC_000913/3958167‑3958353
                                                                                                  |                                                                                                
CGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑TC                                                                                                  >  SRR4302535.983777/1‑94 (MQ=60)
CGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC                                                                                                  >  SRR4302535.3351905/1‑94 (MQ=60)
 GATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑C                                                                                               >  SRR4302535.1211083/1‑94 (MQ=60)
 GATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑C                                                                                               >  SRR4302535.3715114/1‑94 (MQ=60)
  ATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CC                                                                                              >  SRR4302535.483457/1‑94 (MQ=60)
      CGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCC                                                                                          >  SRR4302535.4431327/1‑94 (MQ=60)
       GGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCC                                                                                         >  SRR4302535.3749349/1‑94 (MQ=60)
          CTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccc                                                                                      >  SRR4302535.130773/1‑91 (MQ=60)
           TTCTTTCTTTGATTTAAACACCCCCAACAAACCGGTTAGGGGCTTACCCGCGGTCCCATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccgt                                                                                     <  SRR4302535.3987308/94‑5 (MQ=37)
            TCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccctt                                                                                    >  SRR4302535.1900257/1‑89 (MQ=60)
              TCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccctcg                                                                                  >  SRR4302535.3376176/1‑87 (MQ=60)
               CTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccttcga                                                                                 >  SRR4302535.4611050/1‑86 (MQ=60)
                TCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAG                                                                                 <  SRR4302535.589677/94‑1 (MQ=60)
                 CTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccgctcgaca                                                                               >  SRR4302535.2905484/1‑84 (MQ=60)
                 CTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccttcgaaa                                                                               >  SRR4302535.2024370/1‑84 (MQ=60)
                  TCTGATTTAAACATCCGCAGCCAACCGGGTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGA                                                                                <  SRR4302535.852265/94‑1 (MQ=60)
                   TTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACCCGGGGTCACATTCAAATGCGATTCTGCTACAACCC‑‑CC‑‑CCCCCCCGTTCGAAGAT                                                                               <  SRR4302535.2138596/94‑1 (MQ=60)
                   CTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATT                                                                              <  SRR4302535.3227268/94‑1 (MQ=60)
                    TGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATT                                                                              <  SRR4302535.516420/94‑1 (MQ=60)
                      ATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTTCACGCGGGCACATTCAAATGCGATTCTGCTACAAACCC‑‑CC‑‑CCCCCCCGTTCGAAGATTG                                                                             <  SRR4302535.3175530/94‑1 (MQ=60)
                      ATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAG                                                                           >  SRR4302535.3126205/1‑94 (MQ=60)
                       TTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAACCC‑‑CC‑‑CCCCCCCGTTCGAAGATTGAG                                                                           <  SRR4302535.4655061/94‑1 (MQ=60)
                       TTTAAACACCCCCGGCCAACCGGTTAGCGCCTTACCCGCGGTCACCATCAAATGCGATTCTGCCACAACCC‑‑CC‑‑CCCCCCCGTTCGAAGATTGAG                                                                           <  SRR4302535.4488341/94‑1 (MQ=25)
                        TTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccgtccaagattgagg                                                                        >  SRR4302535.1897456/1‑77 (MQ=60)
                           AACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACCTTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CCcgccccccgcgtcggcgctataccaaa                                                                       >  SRR4302535.2479614/1‑67 (MQ=60)
                            ACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CCcaccccccccttcacacactgtcacata                                                                      >  SRR4302535.2820450/1‑66 (MQ=60)
                             CATCCGCCGCCAACCGGGTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATAC                                                                    <  SRR4302535.997935/94‑1 (MQ=60)
                              CTCCGCAGCCAACCGGTTAGCCGCTTACACGCGGTAACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACA                                                                   <  SRR4302535.1549322/94‑1 (MQ=60)
                              ATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CCCACCCCCCCGTTCGAAGATTGAGCAATACA                                                                   <  SRR4302535.4261008/94‑1 (MQ=60)
                              ATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccctccgcaaattagccataaaa                                                                  >  SRR4302535.3394621/1‑71 (MQ=60)
                               ACCGCAGCCAACCGGTTAGCGGCTTACCCGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACAC                                                                  <  SRR4302535.1525302/94‑1 (MQ=60)
                                    CGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATG                                                             <  SRR4302535.3299837/94‑1 (MQ=60)
                                    AGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATG                                                             <  SRR4302535.3980717/94‑1 (MQ=60)
                                    AGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATC                                                             >  SRR4302535.588762/1‑94 (MQ=60)
                                      CCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCG                                                           <  SRR4302535.4101991/94‑1 (MQ=60)
                                       CAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccctcctaaaaataagaccaaacacatatgct                                                         >  SRR4302535.4009678/1‑62 (MQ=60)
                                         ACCGGTTAGGGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCT                                                        <  SRR4302535.4461296/94‑1 (MQ=60)
                                            cgttagcggctttccagcgggccCATTCAAATGGGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAA                                                     <  SRR4302535.4529898/71‑1 (MQ=60)
                                            GGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccttcgagacttggcaataaaaccattcggcttcac                                                    >  SRR4302535.4237131/1‑57 (MQ=60)
                                             GTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACC                                                   <  SRR4302535.690931/94‑1 (MQ=60)
                                             GTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccggtcgaaagtttagaaatacacacagggttctaaac                                                   >  SRR4302535.928143/1‑56 (MQ=60)
                                              TTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCC                                                  <  SRR4302535.657527/94‑1 (MQ=60)
                                                  CGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGG                                               <  SRR4302535.2484797/94‑1 (MQ=60)
                                                  CGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCAATACGTCTAAACCCCGC                                               >  SRR4302535.2407445/1‑94 (MQ=60)
                                                    GCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccctcgcgaaataacgaaaaaacccatctcttcacacccccccc                                            >  SRR4302535.3396798/1‑49 (MQ=60)
                                                    GCTTACACGCGGTCACATTCAAATGCGATTATGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATAAACCTATGGGCCTAAACCCCGGCCA                                            >  SRR4302535.4130396/1‑94 (MQ=60)
                                                       taccgccggtacaATTCAAATGGGATTCTGCTCCAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACA                                         <  SRR4302535.3879931/81‑1 (MQ=60)
                                                          ACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGCTCGAAGATTAAGCAAAACACCTATTCGTCTAAACCCCGcgcaacaaa                                       >  SRR4302535.194176/1‑85 (MQ=60)
                                                           CGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACA                                      <  SRR4302535.4607990/94‑1 (MQ=60)
                                                            GCGGTCACATTCAAATGCGATTTCGCTAAAAC‑‑C‑‑CC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAG                                    <  SRR4302535.2132806/94‑1 (MQ=60)
                                                            GCGGTCACATTCAAATGCGATTCTGTTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAA                                     <  SRR4302535.4504898/94‑1 (MQ=60)
                                                             CGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccctccgaagatggggcaaaacctcacccgtctcacccccgcccacacaaaac                                   >  SRR4302535.2522563/1‑40 (MQ=60)
                                                             CGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccctccaagcttgggaaaaacacctatgcgcctcaaaaccccccacaacaaca                                   >  SRR4302535.2666406/1‑40 (MQ=60)
                                                             CGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccccccccccattcgcgaataaaaaaaatcattaaacccacaccccaaaaaaa                                   >  SRR4302535.3510603/1‑40 (MQ=60)
                                                              GGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCccggtcgaagatttagaaaaacacattatcgtttaaacgtcccccaataacacgc                                  >  SRR4302535.4284713/1‑39 (MQ=60)
                                                              GGTCACATTCAAAAGCGATTCTGCTACAACCC‑‑CC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGC                                   <  SRR4302535.2334895/94‑1 (MQ=60)
                                                               GTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCT                                  <  SRR4302535.829918/94‑1 (MQ=60)
                                                               GTCACATTCAAATGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTG                                 >  SRR4302535.773698/1‑94 (MQ=60)
                                                               GTCACATTCAAATGCCATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCT                                  <  SRR4302535.3364553/94‑1 (MQ=60)
                                                                TCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCcccttcgacgatgatgaataacaaccaacggtctaaacacccgccacccaccagagt                                >  SRR4302535.1947568/1‑37 (MQ=60)
                                                                 CACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGGGAAATCAACCAATGCGGCTAAACCCCGGCCAAaaaaaaagtgg                                >  SRR4302535.1380872/1‑83 (MQ=33)
                                                                  ACATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGGCAACTACACCCATGGGTCTAAAACCCGGCCAAAAACAAGCTtca                               >  SRR4302535.4400146/1‑91 (MQ=37)
                                                                   CATTCAACTGCGATTCTGCTACAA‑‑‑TCCC‑‑CCCCCCCGTTCGAAGATTGAGGAATACACCTATGCGTGTAAACCCCGGCCAACAACAAGCTGTCGA                             >  SRR4302535.41795/1‑94 (MQ=60)
                                                                   CATTCAAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCCTTCGACGATTGAGGAATATACCCACGCGTCTAAACCCCGGCCAACAAaaaactggcg                              >  SRR4302535.2777637/1‑84 (MQ=35)
                                                                      TCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAAATC                         >  SRR4302535.2788309/1‑94 (MQ=60)
                                                                       TAAATGCGAATCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATT                          <  SRR4302535.4564870/94‑1 (MQ=60)
                                                                        AAATGCGATTCTGCTACAAG‑‑C‑‑CC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCG                        <  SRR4302535.2623058/94‑1 (MQ=60)
                                                                        AAATGCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGGGCAATACACCcctgggccgaacccccagcccacacaaagctgccgaatcc                         >  SRR4302535.580347/1‑54 (MQ=38)
                                                                        AAAGGCGATTCTGCAACAAC‑‑C‑‑CC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCG                        <  SRR4302535.2611230/94‑1 (MQ=60)
                                                                            GCGATTCTGCTACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTA                     <  SRR4302535.2188909/94‑1 (MQ=60)
                                                                             CGATTCTGCTACAAACC‑‑CC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTAC                    <  SRR4302535.3452664/94‑1 (MQ=60)
                                                                               ATTCTGCCACAA‑‑‑TCCCC‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCG                  <  SRR4302535.2282688/94‑1 (MQ=60)
                                                                                                  C‑‑CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGCCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTCCCTGGGGCT  >  SRR4302535.3982444/1‑94 (MQ=60)
                                                                                                  |                                                                                                
CGATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑TCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGCT  >  NC_000913/3958167‑3958353

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