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: