Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,921 T→C 15.8% Q75R (CAA→CGA) ‡ higA ← antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,9210TC15.8% 346.7 / 52.9 147Q75R (CAA→CGA) higAantitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (55/68);  new base C (8/15);  total (64/83)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 151 total aligned reads displayed.

GACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGA  >  NC_000913/3233822‑3234020
                                                                                                   |                                                                                                   
gACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCtt                                                                                                     <  1:2761797/100‑1 (MQ=255)
   acCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGTTTCAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAg                                                                                                  <  2:2564043/100‑1 (MQ=255)
    ccATAGATTTACTGCCCATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGc                                                                                                 <  2:2129729/100‑1 (MQ=255)
    ccATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGc                                                                                                 >  2:2921394/1‑100 (MQ=255)
     cATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGcc                                                                                                <  2:786890/100‑1 (MQ=255)
     cATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGcc                                                                                                >  1:1433485/1‑100 (MQ=255)
     cATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGAGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGcc                                                                                                >  1:2569354/1‑100 (MQ=255)
      aTAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCa                                                                                               >  1:1773248/1‑100 (MQ=255)
       tAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCAt                                                                                              <  2:823265/100‑1 (MQ=255)
        aGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATg                                                                                             >  1:155985/1‑100 (MQ=255)
          aTTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGccct                                                                                              >  2:3638798/1‑95 (MQ=255)
           tttACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcct                                                                                          <  2:1711803/100‑4 (MQ=255)
            ttACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt                                                                                         <  1:3080873/100‑1 (MQ=255)
            ttACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctg                                                                                         <  1:2623201/100‑5 (MQ=255)
            ttACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCAGATTCATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctg                                                                                         <  1:3251149/100‑5 (MQ=255)
             tACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAtt                                                                                        <  1:1580818/100‑1 (MQ=255)
              aCTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTa                                                                                       <  2:1353961/100‑1 (MQ=255)
              aCTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctggc                                                                                       <  1:3084397/100‑7 (MQ=255)
              aCTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctggc                                                                                       <  1:665165/100‑7 (MQ=255)
               cTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTaa                                                                                      <  2:1489258/100‑1 (MQ=255)
               cTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctggcg                                                                                      <  2:2624457/100‑8 (MQ=255)
                tGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTaaa                                                                                     >  1:918331/1‑100 (MQ=255)
                tGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctggcgg                                                                                     <  1:2099276/100‑9 (MQ=255)
                 gCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAAt                                                                                    >  1:1405044/1‑100 (MQ=255)
                 gCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAAt                                                                                    >  1:2225489/1‑100 (MQ=255)
                 gCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAAt                                                                                    <  2:884336/100‑1 (MQ=255)
                 gCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAAt                                                                                    <  1:101498/100‑1 (MQ=255)
                 gCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctggcg                                                                                      <  1:1533409/98‑8 (MQ=255)
                 gCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGCCATGGCTTGAGCCATGGCATTAAAt                                                                                    >  2:2092958/1‑100 (MQ=255)
                   cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                                                                                  >  2:2046986/1‑100 (MQ=255)
                   cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                                                                                  >  2:2041088/1‑100 (MQ=255)
                    aaTTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc                                                                                 <  2:410997/100‑1 (MQ=255)
                    aaTTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc                                                                                 >  1:3146537/1‑100 (MQ=255)
                     aTTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCg                                                                                >  1:211571/1‑100 (MQ=255)
                      tttCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGc                                                                               <  2:3522221/100‑1 (MQ=255)
                         ccGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCaaa                                                                            >  1:2605452/1‑100 (MQ=255)
                         ccGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCACaaa                                                                            >  2:3233729/1‑100 (MQ=255)
                          cGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAAt                                                                           >  1:1498856/1‑100 (MQ=255)
                            gCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTc                                                                         <  1:1735212/100‑1 (MQ=255)
                             cAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCg                                                                        <  2:2987285/100‑1 (MQ=255)
                               gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCggg                                                                      <  1:2979384/100‑1 (MQ=255)
                               gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCggg                                                                      <  1:1387124/100‑1 (MQ=255)
                               gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCggg                                                                      >  1:1003311/1‑100 (MQ=255)
                               gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTAggg                                                                      >  1:961932/1‑100 (MQ=255)
                                  cGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                   <  1:1893558/100‑1 (MQ=255)
                                   ggAAAGGGTTAAACCATATTGATCCATAAGGGTACGAAGCACGGCTAGACCGCCAGGCATGGGTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCt                                                                  >  2:2703769/1‑100 (MQ=255)
                                    gAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg                                                                 >  1:2380111/1‑100 (MQ=255)
                                         ggTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCAGGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttc                                                             <  2:333958/99‑1 (MQ=255)
                                          gTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttctt                                                           >  1:3168786/1‑100 (MQ=255)
                                          gTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttctt                                                           >  1:3387824/1‑100 (MQ=255)
                                             aaaCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTccc                                                        <  1:2004267/100‑1 (MQ=255)
                                              aaCCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCa                                                       <  1:2778778/100‑1 (MQ=255)
                                                ccATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcg                                                     <  1:362344/100‑1 (MQ=255)
                                                  atatTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCAcgcg                                                   <  2:1287063/100‑1 (MQ=255)
                                                    atTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGt                                                 <  2:2569338/100‑1 (MQ=255)
                                                      tGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTa                                               >  1:3628138/1‑100 (MQ=255)
                                                      tGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTa                                               >  2:2348319/1‑100 (MQ=255)
                                                      tGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTa                                               <  1:3009580/100‑1 (MQ=255)
                                                      tGATCCATAAGGGTACGAAGCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTa                                               <  1:432215/100‑1 (MQ=255)
                                                       gATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAt                                              >  1:1995329/1‑100 (MQ=255)
                                                         tCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAttt                                            <  2:3189401/100‑1 (MQ=255)
                                                          ccATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAtttt                                           <  1:9301/100‑1 (MQ=255)
                                                          ccATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTAtttt                                           <  1:1215353/100‑1 (MQ=255)
                                                           cATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTg                                          >  1:2678527/1‑100 (MQ=255)
                                                           cATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTg                                          <  2:1434765/100‑1 (MQ=255)
                                                           cATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTg                                          <  1:2879062/100‑1 (MQ=255)
                                                            aTAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTgg                                         <  2:3239224/100‑1 (MQ=255)
                                                             tAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGc                                        <  1:3632570/100‑1 (MQ=255)
                                                             tAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGc                                        <  1:1786869/100‑1 (MQ=255)
                                                              aaGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGcn                                       <  1:3478402/100‑2 (MQ=255)
                                                                gggTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGcaca                                     >  2:2416357/1‑100 (MQ=255)
                                                                   tACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCacacac                                  >  1:1606643/1‑100 (MQ=255)
                                                                    aCGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACAcc                                 <  2:3586347/100‑1 (MQ=255)
                                                                     cGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCa                                <  2:929597/100‑1 (MQ=255)
                                                                     cGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCa                                >  2:1684582/1‑100 (MQ=255)
                                                                      gAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACAcc                                 >  2:36324/1‑98 (MQ=255)
                                                                        aTCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGAt                             >  1:1228003/1‑100 (MQ=255)
                                                                          cACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATcc                           >  1:1327960/1‑100 (MQ=255)
                                                                          cACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATcc                           >  1:1804922/1‑100 (MQ=255)
                                                                           aCGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCa                          >  1:193405/1‑100 (MQ=255)
                                                                            cGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAg                         <  2:3436122/100‑1 (MQ=255)
                                                                              gCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCa                       <  1:821656/100‑1 (MQ=255)
                                                                              gCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCa                       <  2:602405/100‑1 (MQ=255)
                                                                                tataCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAg                     <  1:210718/100‑1 (MQ=255)
                                                                                tataCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAg                     <  1:432565/100‑1 (MQ=255)
                                                                                 ataCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAgg                    <  2:81830/100‑1 (MQ=255)
                                                                                 ataCCGCCAGGAATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAgg                    <  2:1566394/100‑1 (MQ=255)
                                                                                  taCCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCCCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAggg                   <  2:1980123/100‑1 (MQ=255)
                                                                                        atgccATGGCTCAAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTc             >  1:2871655/5‑100 (MQ=255)
                                                                                        atgccATGGCTCAAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTc             >  1:628924/5‑100 (MQ=255)
                                                                                          ggCATGGCTTGAGCCATGGCATTAAATTCCGCATGTTCGGGCGCGGATTCTTCCCACGCGGTTATTTTGGCACACTCCAGATCCAGCAAGGGGTTTTccg           <  2:1920348/100‑3 (MQ=255)
                                                                                          ggCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTctg           <  1:2458930/100‑3 (MQ=255)
                                                                                          ggCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAg           <  2:2802343/100‑1 (MQ=255)
                                                                                          ggCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAg           <  1:3530638/100‑1 (MQ=255)
                                                                                            cATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGa         >  1:3281782/1‑100 (MQ=255)
                                                                                            cATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGa         <  1:3419375/100‑1 (MQ=255)
                                                                                            cATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGa         >  1:1320008/1‑100 (MQ=255)
                                                                                             aTGGCTCAAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGAt        >  2:513075/1‑100 (MQ=255)
                                                                                               ggCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCg      >  2:400350/1‑100 (MQ=255)
                                                                                                   tGAGTCATGGCATTAAATTCCGCAAATCCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCg      <  2:1152846/96‑1 (MQ=255)
                                                                                                   |                                                                                                   
GACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGA  >  NC_000913/3233822‑3234020

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