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,261 T→C intergenic (‑830/‑31) 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,2610TC98.4% 185.6 / ‑5.6 62intergenic (‑830/‑31)ppiC/reppeptidyl‑prolyl cis‑trans isomerase C/rep helicase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (28/33);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑T‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCC  >  NC_000913/3958174‑3958341
                                                                                               |                                                                                       
ggTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑cc‑‑c‑‑ccccc                                                                                  <  2:4501298/93‑1 (MQ=255)
ggTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑cc‑‑c‑‑ccccc                                                                                  >  1:4401209/1‑93 (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)
                          cGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACAcc                                                        <  2:2260953/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)
                                 aaCCGGTTAGGTGCGTACACGCGGTCAAATTCAAAAGCGATTCTGCTACAAGCC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑t                                              <  2:3402559/93‑1 (MQ=255)
                                  aCCGGTTAGGGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑Tc                                             <  1:4430935/93‑1 (MQ=255)
                                  aCCGGTTAGCTGCTTACACGCAGTCAAACTTAAACGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑Tc                                             <  2:946589/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)
                                       ttAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAAcc                                       <  1:652972/93‑1 (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)
                                                tACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCcaa                               <  2:1185634/93‑1 (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)
                                                     gcgGTCACATTCAAATGCGATTCTGTTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCcaacaaca                          <  1:4474197/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)
                                                           acaTTCAAATGCGATTCTTCTACAA‑‑TCCCC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGt                    <  2:4551667/93‑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)
                                                                 aaaGGCGATTCTGCAACAACC‑‑CC‑‑C‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTc             <  1:2592607/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)
                                                                                               |                                                                                       
GGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑T‑‑CC‑‑T‑‑CCCCCCGTTCGAAGATTGAGCAATACACCTATGCG‑‑TCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCC  >  NC_000913/3958174‑3958341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ 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
GATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑TCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGCTG  >  NC_000913/3958168‑3958354
                                                                                                    |                                                                                              
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)
                            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)
                                                                                                    CCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGCTG  <  SRR4302535.2089538/94‑1 (MQ=60)
                                                                                                    |                                                                                              
GATTTCGGTCTTCTCTCTCTGATTTAAACATCCGCAGCCAACCGGTTAGCGGCTTACACGCGGTCACATTCAAATGCGATTCTGCTACAA‑‑‑T‑‑CC‑‑TCCCCCCGTTCGAAGATTGAGCAATACACCTATGCGTCTAAACCCCGGCCAACAACAAGCTGTCGAATTCGTTACCGGCCCCTGCCTGGTGCTG  >  NC_000913/3958168‑3958354

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