Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,922 G→A 15.3% Q75* (CAA→TAA) ‡ 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,9220GA15.3% 368.7 / 48.7 149Q75* (CAA→TAA) higAantitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (57/69);  new base A (8/15);  total (65/84)
Fisher's exact test for biased strand distribution p-value = 3.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 151 total aligned reads displayed.

ACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAG  >  NC_000913/3233823‑3234021
                                                                                                   |                                                                                                   
  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:2939827/100‑1 (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)
           ttACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTCAAGCCATgcctg                                                                                          <  1:2623201/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)
                 ccAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAAtt                                                                                    >  1:3124797/1‑100 (MQ=255)
                  cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                                                                                   >  2:2046986/1‑100 (MQ=255)
                  cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                                                                                   <  1:3079766/100‑1 (MQ=255)
                  cAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                                                                                   >  2:2041088/1‑100 (MQ=255)
                   aaTTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTcc                                                                                  <  2:410997/100‑1 (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:1387124/100‑1 (MQ=255)
                              gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCggg                                                                       >  1:1003311/1‑100 (MQ=255)
                              gATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTAggg                                                                       >  1:961932/1‑100 (MQ=255)
                                tCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcg                                                                     <  1:2862246/100‑1 (MQ=255)
                                 cGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGgcgc                                                                    <  1:1893558/100‑1 (MQ=255)
                                  ggAAAGGGTTAAACCATATTGATCCATAAGGGTACGAAGCACGGCTAGACCGCCAGGCATGGGTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCt                                                                   >  2:2703769/1‑100 (MQ=255)
                                   gAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTg                                                                  >  1:2380111/1‑100 (MQ=255)
                                     aaGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAt                                                                <  1:328398/100‑1 (MQ=255)
                                        ggTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCAGGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttc                                                              <  2:333958/99‑1 (MQ=255)
                                         gTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttctt                                                            >  1:3387824/1‑100 (MQ=255)
                                         gTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGAttctt                                                            >  1:3168786/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:3246521/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                        <  2:602405/100‑1 (MQ=255)
                                                                             gCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCa                        <  1:821656/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: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)
                                                                                                   |                                                                                                   
ACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCTGATTCTTCCCACGCGGTTATTTTGGCACACACCAGATCCAGCAAGGGGTTTTCAGGATCGTTGAG  >  NC_000913/3233823‑3234021

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: