Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,385,848 G→A D367D (GAC→GAT yjeP ← putative periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,385,8480GA95.9% 727.9 / 6.5 220D367D (GAC→GATyjePputative periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (2/7);  new base A (99/111);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 223 total aligned reads displayed.

GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA  >  NC_000913/4385758‑4385938
                                                                                          |                                                                                          
gCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATc                                                                                          >  2:830132/1‑93 (MQ=255)
gCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATc                                                                                          <  2:572559/93‑1 (MQ=255)
 ccACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCg                                                                                         <  2:554483/93‑1 (MQ=255)
  cACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCgg                                                                                        >  1:4602498/1‑93 (MQ=255)
  cACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCgg                                                                                        <  1:3262023/93‑1 (MQ=255)
   aCCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCCTCGGc                                                                                       >  2:826923/1‑93 (MQ=255)
   aCCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGc                                                                                       >  1:1010888/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  1:2813400/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  2:3602972/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      <  1:1728880/93‑1 (MQ=255)
    cccTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  1:1409786/1‑93 (MQ=255)
    cccTGCAACAATGAGTTCAGCAAATCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGcc                                                                                      >  2:442846/1‑93 (MQ=255)
       tGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtga                                                                                   >  2:456720/1‑93 (MQ=255)
       tGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtga                                                                                   >  1:4048760/1‑93 (MQ=255)
        gCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtgat                                                                                  <  1:3119779/93‑1 (MQ=255)
        gCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGCCAGCGGCTGACCATCGGCCtgat                                                                                  >  1:4172361/1‑93 (MQ=255)
          aacaaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCtgatga                                                                                >  1:1887419/1‑93 (MQ=255)
           acaaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGaa                                                                               >  1:4642854/1‑93 (MQ=255)
           acaaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGaa                                                                               <  1:2703781/93‑1 (MQ=255)
            caaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAt                                                                              <  1:1352048/93‑1 (MQ=255)
            caaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAt                                                                              >  2:924536/1‑93 (MQ=255)
            caaTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGAACATCGGCCTGATGAAt                                                                              >  1:1378886/1‑93 (MQ=255)
              aTGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAAttt                                                                            >  2:1879480/1‑93 (MQ=255)
               tGAGTTCAGCAACTCACGCTGAGTCCGCAGTTGTGCTCCCAGAATACGGTTTTGCCCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTg                                                                           <  2:1018138/93‑1 (MQ=255)
                  gTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTTAGCGGCTGACCATCGGCCTGATGAATTTGCCg                                                                        <  1:3283229/93‑1 (MQ=255)
                    tcatcaACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTa                                                                      <  2:3313634/93‑1 (MQ=255)
                    tCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTa                                                                      <  1:2192791/93‑1 (MQ=255)
                      aGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTagc                                                                    <  1:3952895/93‑1 (MQ=255)
                       gCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTagca                                                                   <  1:3674675/93‑1 (MQ=255)
                       gCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTagca                                                                   >  2:3701412/1‑93 (MQ=255)
                          aCTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCg                                                                >  1:4652638/1‑93 (MQ=255)
                           cTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCgg                                                               >  2:3041260/1‑93 (MQ=255)
                           cTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCgg                                                               >  1:4375676/1‑93 (MQ=255)
                           cTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCgg                                                               >  1:4314465/1‑93 (MQ=255)
                            tCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGc                                                              >  2:1359834/1‑93 (MQ=255)
                             cACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCt                                                             <  2:4061282/93‑1 (MQ=255)
                                gCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGtt                                                          >  2:1492420/1‑93 (MQ=255)
                                 cTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGttt                                                         >  2:4687706/1‑93 (MQ=255)
                                 cTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGttt                                                         <  2:3458773/93‑1 (MQ=255)
                                   taGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTttat                                                       <  2:3651904/92‑1 (MQ=255)
                                    aGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTttatt                                                      <  1:1655271/93‑1 (MQ=255)
                                    aGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTttatt                                                      >  1:1146568/1‑93 (MQ=255)
                                      tGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTaa                                                    <  1:1079913/93‑1 (MQ=255)
                                      tGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTaa                                                    <  1:4193573/93‑1 (MQ=255)
                                       gcgcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAg                                                   <  2:2494830/93‑1 (MQ=255)
                                       gcgcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAg                                                   >  1:2817489/1‑93 (MQ=255)
                                       gcgcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGGTTATTAAg                                                   >  2:203381/1‑93 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  1:1630289/93‑1 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  2:1335080/93‑1 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  1:680261/93‑1 (MQ=255)
                                         gcAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCa                                                 <  2:822492/93‑1 (MQ=255)
                                          cAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAg                                                <  2:4264843/93‑1 (MQ=255)
                                          cAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAg                                                <  2:37237/93‑1 (MQ=255)
                                             ttGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATc                                             >  1:3699739/1‑93 (MQ=255)
                                             ttGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATc                                             >  1:1224811/1‑93 (MQ=255)
                                             ttGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATc                                             <  1:3979022/93‑1 (MQ=255)
                                              tgtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATcc                                            <  1:2842838/93‑1 (MQ=255)
                                               gtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCt                                           <  2:1897949/93‑1 (MQ=255)
                                               gtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCt                                           <  1:468288/93‑1 (MQ=255)
                                               gtgCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCACATCCt                                           >  2:684249/1‑93 (MQ=255)
                                                   ttCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATa                                       >  1:3486563/1‑93 (MQ=255)
                                                    tCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATaa                                      <  2:2944777/93‑1 (MQ=255)
                                                     ccAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAAc                                     <  2:1431386/93‑1 (MQ=255)
                                                      cAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACg                                    >  1:1797473/1‑93 (MQ=255)
                                                       aGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGt                                   >  1:2542280/1‑93 (MQ=255)
                                                       aGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGt                                   >  1:1994593/1‑93 (MQ=255)
                                                        gAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTa                                  >  1:2405596/1‑93 (MQ=255)
                                                         aaTACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTaa                                 <  1:4155987/93‑1 (MQ=255)
                                                         aaTACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTaa                                 >  2:504973/1‑93 (MQ=255)
                                                          aTACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAc                                <  1:1840370/93‑1 (MQ=255)
                                                          aTACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAc                                <  2:2631402/93‑1 (MQ=255)
                                                          aTACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAc                                <  2:211421/93‑1 (MQ=255)
                                                           tACGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAAcc                               >  1:3871530/1‑93 (MQ=255)
                                                             cGGTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGt                             <  2:2193868/93‑1 (MQ=255)
                                                              ggTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGtt                            <  1:215633/93‑1 (MQ=255)
                                                               gTTTTGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGGTg                           >  2:3261404/1‑93 (MQ=255)
                                                                ttttGCTCTGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGc                          <  2:1264732/93‑1 (MQ=255)
                                                                ttttGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGACGTAGCAGCGGCTGTTTATTAAGCAGATACTCATAACGTAACCGTTGc                          >  1:3398963/1‑93 (MQ=255)
                                                                   tGCTCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGcaca                       <  2:3394836/93‑1 (MQ=255)
                                                                      tCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGc                    >  2:1626311/1‑93 (MQ=255)
                                                                      tCGGCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGc                    <  2:3763511/93‑1 (MQ=255)
                                                                        ggCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCaa                  >  2:1629668/1‑93 (MQ=255)
                                                                         gCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATCAAGTAACCGTTTCACCCGCAAc                 >  2:3549336/1‑93 (MQ=255)
                                                                         gCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAc                 >  2:2272761/1‑93 (MQ=255)
                                                                         gCAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAAc                 <  1:2528306/93‑1 (MQ=255)
                                                                          cAGTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACt                <  1:2635094/93‑1 (MQ=255)
                                                                            gTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgc              <  2:2318150/93‑1 (MQ=255)
                                                                            gTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgc              <  1:4343663/93‑1 (MQ=255)
                                                                            gTCAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgc              <  2:268468/93‑1 (MQ=255)
                                                                              cAGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTgcgc            >  2:4402126/1‑93 (MQ=255)
                                                                               aGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGcc           <  1:36193/93‑1 (MQ=255)
                                                                               aGCGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGcc           <  1:2513614/93‑1 (MQ=255)
                                                                                 cGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAt         <  1:293155/93‑1 (MQ=255)
                                                                                 cGGCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAt         >  1:3344094/1‑93 (MQ=255)
                                                                                   gCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAttt       >  1:4391576/1‑93 (MQ=255)
                                                                                   gCTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCAttt       >  2:1078076/1‑93 (MQ=255)
                                                                                    cTGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTc      >  2:1468743/1‑93 (MQ=255)
                                                                                     tGACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCg     >  1:1269600/1‑93 (MQ=255)
                                                                                      gACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCgg    <  1:1488850/93‑1 (MQ=255)
                                                                                      gACCATCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCgg    <  2:1702087/93‑1 (MQ=255)
                                                                                          |                                                                                          
GCCACCCTGCAACAATGAGTTCAGCAACTCACGCTGAGTGCGCAGTTGTGCTTCCAGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTA  >  NC_000913/4385758‑4385938

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