Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,956,638 C→A 5.5% V54V (GTG→GTT fliI ← flagellum‑specific ATP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,956,6380CA 5.5% 236.0 / 2.5 73V54V (GTG→GTTfliIflagellum‑specific ATP synthase
Reads supporting (aligned to +/- strand):  ref base C (39/30);  new base A (2/2);  total (41/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

CAGCGGCACGACCCGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTGGG  >  NC_002947/4956491‑4956786
                                                                                                                                                   |                                                                                                                                                    
cAGCGGCACGACCCGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCAc                                                                                                                                                    >  2:1886200/1‑150 (MQ=255)
 aGCGGCACGACCCGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCAct                                                                                                                                                   <  2:15865/150‑1 (MQ=255)
      cACGACCCGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACtgctgc                                                                                                                                              >  1:1456661/1‑150 (MQ=255)
       aCGACCCGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTgcg                                                                                                                                             >  1:1545298/1‑150 (MQ=255)
             cGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAgg                                                                                                                                       >  1:1521517/1‑150 (MQ=255)
                 cccccGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGGTCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCtt                                                                                                                                   >  2:1241689/1‑150 (MQ=255)
                   cccGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCg                                                                                                                                 <  1:1886200/150‑1 (MQ=255)
                    ccGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCgg                                                                                                                                <  2:1066525/150‑1 (MQ=255)
                          cgATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCa                                                                                                                          >  2:1809756/1‑150 (MQ=255)
                           gATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTNCAg                                                                                                                         <  2:1202617/150‑1 (MQ=255)
                            aTGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAgg                                                                                                                        <  1:1844286/150‑1 (MQ=255)
                              gCCGACGATACTGCCCACCGGCATCCGAAACACCTTGGGCCCGGCCCCATTTTTGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCCCCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGt                                                                                                                      <  1:1057978/150‑1 (MQ=255)
                              gCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGt                                                                                                                      <  1:167207/150‑1 (MQ=255)
                                cgacgaTACTGCCCACCGGCATCAGAAACGCCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGa                                                                                                                    <  2:1123840/150‑1 (MQ=255)
                                 gacgaTACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAg                                                                                                                   <  1:1489299/150‑1 (MQ=255)
                                    gaTACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGcc                                                                                                                >  2:1383960/1‑150 (MQ=255)
                                      tACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGa                                                                                                              >  2:1759746/1‑150 (MQ=255)
                                       aCTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGAc                                                                                                             <  2:1865058/150‑1 (MQ=255)
                                         tGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCa                                                                                                           >  1:524234/1‑150 (MQ=255)
                                                    aTCAGAAACACCTTGGGCCCGGCAAAGCACATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAg                                                                                                >  2:532477/1‑150 (MQ=255)
                                                     tCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACGTGCACCGGGGGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGc                                                                                               >  2:479479/1‑150 (MQ=255)
                                                       aGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCgg                                                                                             >  2:1702611/1‑150 (MQ=255)
                                                       aGAAAAACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCgg                                                                                             >  2:1702287/1‑150 (MQ=255)
                                                           acacCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCtt                                                                                         >  1:553950/1‑150 (MQ=255)
                                                           acacCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCAACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCAt                                                                                                          >  2:1179844/1‑133 (MQ=255)
                                                           acacCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCAACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCAt                                                                                                          >  2:1180122/1‑133 (MQ=255)
                                                           acacCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCAACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCAt                                                                                                          <  1:1179844/133‑1 (MQ=255)
                                                           acacCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCAACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCAt                                                                                                          <  1:1180122/133‑1 (MQ=255)
                                                               cTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTcgac                                                                                     >  2:763820/1‑150 (MQ=255)
                                                               cTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCCTCGTGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTcgac                                                                                     >  1:1453506/1‑150 (MQ=255)
                                                                   ggCCCGGCAAAGCCCATGACTGCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGAtg                                                                                 <  2:845678/150‑1 (MQ=255)
                                                                     cccGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAg                                                                                                                                        >  2:1510382/1‑93 (MQ=255)
                                                                     cccGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAg                                                                                                                                        <  1:1510382/93‑1 (MQ=255)
                                                                        tgCAAATCCCATGAATTAGGACTCGCCATCCACCGAGTGATATCCGCGAACGTCGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGAtgggctg                                                                            <  2:1453506/149‑1 (MQ=255)
                                                                        ggCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGAtgggctg                                                                            <  2:1793928/150‑1 (MQ=255)
                                                                         gCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCCAGCGGCCTTCGACGAtgggctgg                                                                           <  2:943471/150‑1 (MQ=255)
                                                                                       ttCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTAGACGATGGGCTGGGCCGGCAGCTCGAt                                                             >  1:84236/1‑150 (MQ=255)
                                                                                       ttCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGAAAGCGGCCTTCGACGATGGACTGGGCCGCCAGCTCGAt                                                             >  2:1545355/1‑150 (MQ=255)
                                                                                        tCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGGTCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATg                                                            <  1:1241689/150‑1 (MQ=255)
                                                                                        tCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATg                                                            <  1:1759746/150‑1 (MQ=255)
                                                                                          ggCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGACCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGc                                                          >  2:746924/1‑150 (MQ=255)
                                                                                            ccTCGCCCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCt                                                        >  2:852371/1‑150 (MQ=255)
                                                                                            ccTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCt                                                        <  2:1290557/150‑1 (MQ=255)
                                                                                             cTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTc                                                       >  1:1554904/1‑150 (MQ=255)
                                                                                                   cTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTGGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTa                                                 >  2:1090134/1‑150 (MQ=255)
                                                                                                     gCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAg                                               <  1:1545355/150‑1 (MQ=255)
                                                                                                             tGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGa                                       <  2:1456661/150‑1 (MQ=255)
                                                                                                                 aaCTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCt                                   >  1:199875/1‑150 (MQ=255)
                                                                                                                   cTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTg                                 >  1:1243229/1‑150 (MQ=255)
                                                                                                                    tGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGTTGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGc                                <  2:1866001/150‑1 (MQ=255)
                                                                                                                    tGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGTTGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGc                                >  1:1866001/1‑150 (MQ=255)
                                                                                                                    tGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGc                                                                  <  2:504767/116‑1 (MQ=255)
                                                                                                                    tGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGc                                                                  >  1:504767/1‑116 (MQ=255)
                                                                                                                    tGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCCGCGTAAGTGCCCAGACGCTTGc                                >  2:670868/1‑150 (MQ=255)
                                                                                                                      tCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTaa                                                <  2:197169/132‑1 (MQ=255)
                                                                                                                      tCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTaa                                                >  1:197169/1‑132 (MQ=255)
                                                                                                                        aTCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGaa                            >  1:1352423/1‑150 (MQ=255)
                                                                                                                         tCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAg                           >  2:1341676/1‑150 (MQ=255)
                                                                                                                           gTTGATCACCAGGCAGCGGCTGCACACTGCTGCGCGCAGGCCTTCGTCTTTCGGATTTAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCCGCGTAAGTGCCCAGACGCCTGCCGAAGcc                         >  2:738426/1‑149 (MQ=255)
                                                                                                                            ttGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAgttt                        <  1:1809756/150‑4 (MQ=255)
                                                                                                                                 cagtgggCAGCCGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCGGGTGCg                   >  1:937420/6‑150 (MQ=255)
                                                                                                                                 cACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGAtgggctg                                                                            <  1:477245/93‑1 (MQ=255)
                                                                                                                                 cACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGAtgggctg                                                                            >  2:477245/1‑93 (MQ=255)
                                                                                                                                  aCCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGt                  >  2:1346123/1‑150 (MQ=255)
                                                                                                                                    cAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGc                                                    >  1:1364718/1‑114 (MQ=255)
                                                                                                                                    cAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGc                                                    <  2:1364718/114‑1 (MQ=255)
                                                                                                                                    cAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCGGGGCCGGCAGCTCGATGGCCTCAGCGTGAGTGCCCAGACGCTTGCCGAAGCTGGTGCGttg                >  1:634702/1‑149 (MQ=255)
                                                                                                                                     aGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCa               >  1:1625473/1‑150 (MQ=255)
                                                                                                                                      ggCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCGGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCaa              >  2:174296/1‑150 (MQ=255)
                                                                                                                                      ggCAGCGGCTGCCCACTGCTGCGCGCAGGAATTCGGATTAAATTAGCAGCCCTATCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCaa              <  1:1374684/150‑1 (MQ=255)
                                                                                                                                       gCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAg             <  2:1625473/150‑1 (MQ=255)
                                                                                                                                        cAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCTACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAgg            <  1:592355/150‑1 (MQ=255)
                                                                                                                                         aGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGgc           <  1:1155946/150‑1 (MQ=255)
                                                                                                                                                 cccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   <  2:1554904/150‑1 (MQ=255)
                                                                                                                                                 cccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   <  2:1521517/150‑1 (MQ=255)
                                                                                                                                                 cccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   >  2:1423957/1‑150 (MQ=255)
                                                                                                                                                 cccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   <  2:1061741/150‑1 (MQ=255)
                                                                                                                                                 cccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGCGTGCGTGTGGCAGCTCAATGTGGTTAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   <  1:705200/150‑1 (MQ=255)
                                                                                                                                                 cccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACAAAGCCCAAAGCCGGCAGCTTTATGTCTTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   <  1:705210/150‑1 (MQ=255)
                                                                                                                                                 cccACTGATGGGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTgg   >  2:1423918/1‑150 (MQ=255)
                                                                                                                                                  ccACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTggg  <  1:746924/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CAGCGGCACGACCCGGGCCCCCGGCGCGATGCCGACGATACTGCCCACCGGCATCAGAAACACCTTGGGCCCGGCAAAGCCCATGACTTCGGCCTCGACCTGCACCGGGTGATAACTGTCATCGTTGATCACCAGGCAGCGGCTGCCCACTGCTGCGCGCAGGCCTTCGGCTTCCAGGGTGAGGCCGACCATACGCAGCAAGCGGCCTTCGACGATGGGCTGGGCCGGCAGCTCGATGGCCTCAGCGTAAGTGCCCAGACGCTTGCCGAAGCTGGTGCGTTCAAGGCGCATCTGGG  >  NC_002947/4956491‑4956786

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