Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 21.3% intergenic (‑181/+127) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA21.3% 620.1 / 147.9 231intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (85/97);  new base A (44/5);  total (129/102)
Fisher's exact test for biased strand distribution p-value = 1.57e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-02

Only 100 of 243 total aligned reads displayed.

GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070
                                                                                                  |                                                                                                         
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  2:2404183/1‑100 (MQ=255)
gggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcg                                                                                                          >  1:738494/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        >  1:3275436/1‑100 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  1:2987211/100‑1 (MQ=255)
  gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                        <  1:1965877/100‑1 (MQ=255)
    aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                      >  2:2939286/1‑100 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  1:2097709/100‑1 (MQ=255)
      gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                    <  1:86004/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   <  1:3577241/100‑1 (MQ=255)
       aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   >  1:1002174/1‑100 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                   >  2:3462030/1‑99 (MQ=255)
        aGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAAc                                                                                                  <  1:252699/100‑1 (MQ=255)
         gTATTGCCGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 >  2:1940974/1‑100 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                                 <  2:2532002/99‑1 (MQ=255)
          tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  1:3422502/100‑1 (MQ=255)
           aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCAGCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                                <  1:3209110/99‑1 (MQ=255)
            ttGCTGCCGGATGTCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  1:602182/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  2:1642827/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              >  1:2651115/1‑100 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  1:2744593/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              <  2:208250/100‑1 (MQ=255)
            ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                              >  1:1525855/1‑100 (MQ=255)
             tgcggcCGGATTGCGCTGCACGGAGAGGCACAGCGCCACCGGGTATAAAGATTGTTGACAAAGTTCGATTTATTTATGCCAGATGCGGAGTTAACggatt                                                                                             >  1:2296347/1‑96 (MQ=25)
                tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                          <  1:2080456/100‑1 (MQ=255)
                 gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                         <  2:2095945/100‑1 (MQ=255)
                   cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                         >  1:3058352/1‑98 (MQ=255)
                    ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGGCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                      <  2:677890/100‑1 (MQ=255)
                     gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                     <  2:1620752/100‑1 (MQ=255)
                          cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                <  2:2210724/100‑7 (MQ=25)
                             tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                            <  2:2583514/100‑1 (MQ=255)
                              gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                           <  2:441404/100‑1 (MQ=255)
                               cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                          >  1:1604397/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                        <  2:2634915/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                        >  2:2055967/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                        <  2:2370937/100‑1 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                        >  2:260996/1‑100 (MQ=255)
                                 cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                        <  1:3120527/100‑1 (MQ=255)
                                  gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                       <  1:3608290/100‑1 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                      >  1:3561829/1‑100 (MQ=255)
                                   cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCa                                                                      >  2:2339947/1‑100 (MQ=255)
                                      agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                   <  2:1633978/100‑1 (MQ=255)
                                       ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                                  <  2:3226075/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                 <  1:1387367/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                 <  1:2930836/100‑1 (MQ=255)
                                        gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                                 <  1:903891/100‑1 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                >  2:3608372/1‑100 (MQ=255)
                                         cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                >  2:1215504/1‑100 (MQ=255)
                                          acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCaa                                                               >  2:1905814/1‑100 (MQ=255)
                                           caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                              <  2:1107637/100‑1 (MQ=255)
                                            aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCTCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                              <  1:767953/99‑1 (MQ=255)
                                             gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                            >  2:504571/1‑100 (MQ=255)
                                               gcCACCGAGTATGCAGATTGTTGACAAAGTGCTCTTTGTGTATGCCCGATGCGGCGTTAACGCCTTAGCAGGCCTAC‑AAGCTCGTGCAAATTCAAt                                                              <  1:1725959/96‑1 (MQ=17)
                                                ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                         >  1:2629182/1‑100 (MQ=255)
                                                  aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGc                                                       <  1:316794/100‑1 (MQ=255)
                                                     gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGCCCTAC‑AAAATCGTGCAAATTCAATACATTGCAgc                                                    >  1:2570598/1‑99 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑‑t                                              >  1:3164643/1‑100 (MQ=255)
                                                        tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑At‑‑t                                              >  2:3545010/1‑100 (MQ=255)
                                                           aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTc                                           <  1:2951912/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑g                                        <  2:27927/100‑1 (MQ=255)
                                                            cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑g                                        <  1:660342/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑Gt                                       <  1:338780/100‑1 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑Gt                                       >  2:2750626/1‑100 (MQ=255)
                                                             aGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑Gt                                       <  1:98090/100‑1 (MQ=255)
                                                             aGATTGTTGAATAAGTGCGCTTTGTTTATGCCGGATTCGGCGTTAACGCCTTATCCGGACTAC‑AAAATCGTGCAAATTCAATACCTTGCAGA‑AT‑‑TTTC‑‑Gt                                       <  1:1663792/100‑1 (MQ=17)
                                                               aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAg                                     <  2:1525862/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCt                                 <  1:1838717/100‑1 (MQ=255)
                                                                   ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCt                                 >  2:3422694/1‑100 (MQ=255)
                                                                       cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATa                             <  1:2688472/100‑1 (MQ=255)
                                                                       cAAAGTGCGCTTGGTTTATGCCGGATGCGGCGTTCACGTCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATa                             <  2:74374/100‑1 (MQ=255)
                                                                           gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCg                         <  2:3276946/100‑1 (MQ=255)
                                                                              tccgatctttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCATA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAg                      <  2:2254921/93‑1 (MQ=255)
                                                                               gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     >  2:2311046/1‑100 (MQ=255)
                                                                                  atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATC‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCAt                 >  2:1236679/9‑100 (MQ=18)
                                                                                     tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               >  2:1459496/1‑100 (MQ=255)
                                                                                      ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              >  2:3065013/1‑100 (MQ=255)
                                                                                      gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGCA‑A‑TTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCa               >  2:873352/5‑100 (MQ=18)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑AT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg              >  2:115973/5‑100 (MQ=21)
                                                                                      gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGG‑AG‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:2676870/5‑100 (MQ=21)
                                                                                       taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:3142398/4‑100 (MQ=18)
                                                                                       tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAg              <  1:1409824/99‑1 (MQ=255)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGA‑AA‑‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2006942/3‑100 (MQ=18)
                                                                                        aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCAGA‑GT‑‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1558833/3‑100 (MQ=17)
                                                                                        aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAGAGAT‑‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:1026973/3‑100 (MQ=17)
                                                                                         ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATC‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAt          >  1:1062888/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGG‑AC‑‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1375208/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCAGT‑TT‑‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa           >  2:1118406/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGATCG‑AT‑‑CCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:1284612/2‑100 (MQ=21)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCAAT‑CT‑‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:3254941/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTGCAGG‑GA‑‑CTAT‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa           >  1:1164148/2‑100 (MQ=18)
                                                                                         ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCAGA‑GA‑‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2442943/2‑100 (MQ=17)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTGC‑GA‑ATC‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  2:365000/1‑100 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGC‑‑A‑AT‑‑TCTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg          >  1:2134425/1‑98 (MQ=18)
                                                                                          gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTGCATG‑GG‑‑CCAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:913008/1‑99 (MQ=17)
                                                                                          gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCAGG‑AG‑‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  2:1328885/1‑99 (MQ=17)
                                                                                           ccGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  2:3079328/1‑100 (MQ=255)
                                                                                               aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:2233374/100‑1 (MQ=255)
                                                                                                tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc    <  1:2659298/100‑1 (MQ=255)
                                                                                                 gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  1:970168/100‑1 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:1755694/1‑100 (MQ=255)
                                                                                                  cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  1:3470182/1‑100 (MQ=255)
                                                                                                  |                                                                                                         
GGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGA‑AT‑‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705873‑3706070

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