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: