Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,945,977 T→G 5.2% intergenic (‑135/+104) tcdA ← / ← mltA tRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,945,9770TG 5.2% 389.8 / 10.2 135intergenic (‑135/+104)tcdA/mltAtRNA threonylcarbamoyladenosine dehydratase; sulfur acceptor for CsdA/membrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base T (63/65);  new base G (4/3);  total (67/68)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.32e-01

Only 100 of 140 total aligned reads displayed.

AATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACC  >  NC_000913/2945828‑2946126
                                                                                                                                                     |                                                                                                                                                     
aaTTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCTATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                       <  1:1545070/150‑1 (MQ=255)
aaTTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATACCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGtt                                                                                                                                                       >  1:1664145/1‑150 (MQ=255)
 aTTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                      <  2:851024/150‑1 (MQ=255)
   tACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCATACATACCTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                                    >  2:679609/1‑150 (MQ=255)
      cacaGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                                 >  1:1675161/1‑150 (MQ=255)
               aGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                        <  1:1379602/150‑1 (MQ=255)
                gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcc                                                                                                                                       <  2:855352/150‑2 (MQ=37)
                gCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                       >  1:1160137/1‑150 (MQ=255)
                 cACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCCCCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                      >  1:1646548/1‑150 (MQ=255)
                   ccTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGcctac                                                                                                                                    <  1:1112663/150‑5 (MQ=37)
                   ccTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAtac                                                                                                                                    <  2:1180898/150‑3 (MQ=255)
                   ccTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                                                                                                    >  1:546577/1‑150 (MQ=255)
                   ccTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                                                                                                    >  1:1121841/1‑150 (MQ=255)
                     tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                                  <  1:902239/150‑1 (MQ=255)
                     tCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaac                                                                                                                                  >  1:627152/1‑150 (MQ=255)
                         aaTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg                                                                                                                              <  2:1469856/150‑1 (MQ=255)
                          aTTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCAGGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcgc                                                                                                                             >  1:829599/1‑150 (MQ=255)
                              aaaaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCAc                                                                                                                         >  1:528842/1‑150 (MQ=255)
                                aaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACtt                                                                                                                       >  2:19824/1‑150 (MQ=255)
                                aaGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACtt                                                                                                                       >  1:1534231/1‑150 (MQ=255)
                                  gggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAg                                                                                                                     <  1:1150324/150‑1 (MQ=255)
                                  gggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAg                                                                                                                     >  1:101873/1‑150 (MQ=255)
                                   ggCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGt                                                                                                                    >  1:1439827/1‑150 (MQ=255)
                                     cGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCa                                                                                                                  >  1:823712/1‑150 (MQ=255)
                                       aGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                      >  1:826780/1‑128 (MQ=255)
                                       aGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAt                                                                                                                                      <  2:826780/128‑1 (MQ=255)
                                        ggTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCaata                                                                                                               >  1:272474/1‑150 (MQ=255)
                                              cccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                                                         <  2:1446458/150‑1 (MQ=255)
                                              cccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAgg                                                                                                         <  2:497476/150‑1 (MQ=255)
                                               ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                                        >  2:1574169/1‑150 (MQ=255)
                                               ccTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGa                                                                                                        <  1:1651927/150‑1 (MQ=255)
                                                 tcgcccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtg                                                                                                      >  1:713150/4‑150 (MQ=255)
                                                  cACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgt                                                                                                     >  2:1455125/1‑150 (MQ=255)
                                                   aCCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGAtgtg                                                                                                    <  2:272474/150‑1 (MQ=255)
                                                    cccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGa                                                                                                   >  1:1560254/1‑150 (MQ=255)
                                                     ccTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAg                                                                                                  >  2:640526/1‑150 (MQ=255)
                                                        aaaCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATATTCAGGATGTGAGGtt                                                                                               >  2:64548/1‑150 (MQ=255)
                                                          aCGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAc                                                                                             <  2:815210/150‑1 (MQ=255)
                                                             cATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcccc                                                                                          >  2:1516293/1‑150 (MQ=255)
                                                             cATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTAcccc                                                                                          <  2:998340/150‑1 (MQ=255)
                                                              aTAAATTCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCt                                                                                         <  1:603414/150‑1 (MQ=255)
                                                                  aTGCCTGATGCGCTACGCTTATCAGGCCTACATTACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAcc                                                                                     <  1:611218/150‑1 (MQ=255)
                                                                  aTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATTTTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAcc                                                                                     <  1:500243/150‑1 (MQ=255)
                                                                  aTGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAcc                                                                                     <  1:1257621/150‑1 (MQ=255)
                                                                   tGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCAccc                                                                                    <  1:1685207/150‑1 (MQ=255)
                                                                    gCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGGATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCt                                                                                   >  1:1313636/1‑150 (MQ=255)
                                                                    gCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCt                                                                                   >  1:1332953/1‑150 (MQ=255)
                                                                      cTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGt                                                                                 <  1:1067705/150‑1 (MQ=255)
                                                                            cgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATc                                                                           <  1:11714/150‑1 (MQ=255)
                                                                            cgcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATc                                                                           >  2:413716/1‑150 (MQ=255)
                                                                             gcTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTTGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATcc                                                                          >  2:920462/1‑150 (MQ=255)
                                                                              cTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCatagaaa                                                                         >  2:797459/1‑146 (MQ=255)
                                                                               tACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGAGTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGt                                                                        <  2:1693350/150‑1 (MQ=255)
                                                                                   cTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTaaaa                                                                    >  2:619696/1‑150 (MQ=255)
                                                                                     tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                                  >  2:408465/1‑150 (MQ=255)
                                                                                     tATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAAc                                                                  >  2:629907/1‑150 (MQ=255)
                                                                                        tatgCCTACATATCTCATGCAATATTTTGGATTTTCATGATTTGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGGCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                               <  1:1500084/147‑1 (MQ=37)
                                                                                        cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCCCTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATACGTAAAGCCGGc                                                               >  2:260557/1‑150 (MQ=255)
                                                                                        cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                               <  2:756489/150‑1 (MQ=255)
                                                                                        cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                               <  1:615946/150‑1 (MQ=255)
                                                                                        cAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGc                                                               <  2:1483899/150‑1 (MQ=255)
                                                                                         aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                              <  2:999151/150‑1 (MQ=255)
                                                                                         aGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCa                                                              <  1:138543/150‑1 (MQ=255)
                                                                                                 aTAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAat                                                      <  2:713150/150‑1 (MQ=255)
                                                                                                     cTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATca                                                  >  1:1627628/1‑150 (MQ=255)
                                                                                                     cTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATca                                                  >  2:1662372/1‑150 (MQ=255)
                                                                                                      tCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATcac                                                 <  2:733691/150‑1 (MQ=255)
                                                                                                         tGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATc                                              <  1:1506902/150‑1 (MQ=255)
                                                                                                         tGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACACCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATc                                              <  1:679609/150‑1 (MQ=255)
                                                                                                          gCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCa                                             >  2:1048469/1‑150 (MQ=255)
                                                                                                           aaaTATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAg                                            <  1:1629920/149‑1 (MQ=255)
                                                                                                             atatTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCTCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGcc                                          <  2:1086531/150‑1 (MQ=255)
                                                                                                                ttttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCt                                       <  2:1594542/150‑1 (MQ=255)
                                                                                                                 tttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTa                                      <  1:984305/150‑1 (MQ=255)
                                                                                                                  ttGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaa                                     >  2:1180672/1‑150 (MQ=255)
                                                                                                                   tGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTaaa                                    <  2:962541/150‑1 (MQ=255)
                                                                                                                    ggATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAg                                   <  2:823712/150‑1 (MQ=255)
                                                                                                                       tttGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACg                                <  1:239512/150‑1 (MQ=255)
                                                                                                                        ttGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGt                               >  1:349072/1‑150 (MQ=255)
                                                                                                                          gCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTa                             >  2:1100198/1‑150 (MQ=255)
                                                                                                                             tGAGTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATGAACAACGCCCACTTAGTCAAAAATCAGGATTTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGTCATACAGAATATCACATCAGCCGCTAAAGACGTTACCt                          <  2:140828/150‑1 (MQ=38)
                                                                                                                              gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTg                         >  1:904316/1‑150 (MQ=255)
                                                                                                                              gATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTg                         <  2:961873/150‑1 (MQ=255)
                                                                                                                               aTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACTGAAGATCACATCAGCCGCTAAAGACGTTACCTgc                        >  1:128344/1‑150 (MQ=255)
                                                                                                                                tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACGACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAGACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTgcg                       >  1:981194/1‑150 (MQ=255)
                                                                                                                                tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGATACCTgcg                       <  2:468950/150‑1 (MQ=255)
                                                                                                                                   tcatgCCGGATGCGGCGTGAACGCCTTATCCGGCCTGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGccc                    >  1:429094/5‑150 (MQ=38)
                                                                                                                                     atgCCGGATGCGGCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCgg                  >  2:1429262/3‑150 (MQ=39)
                                                                                                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCggg                 >  2:1129389/1‑150 (MQ=255)
                                                                                                                                      ggCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCggg                 >  1:1461322/1‑150 (MQ=255)
                                                                                                                                       gCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCgggg                >  2:417560/1‑150 (MQ=255)
                                                                                                                                        ccGGATAAGGCGTTCACGCCGCATCCGGCATGAAAAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTAATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGc               <  2:138145/150‑1 (MQ=255)
                                                                                                                                           gatgcggCGTGAACGCCTTATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGTCGCTAAAGACGTTACCTGCGCCCGGGGCGGt            >  1:1049819/6‑150 (MQ=255)
                                                                                                                                               aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc        <  2:101873/150‑1 (MQ=255)
                                                                                                                                               aGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTc        <  1:141027/150‑1 (MQ=255)
                                                                                                                                                ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCa       <  2:627152/150‑1 (MQ=255)
                                                                                                                                                ggCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCa       <  1:19824/150‑1 (MQ=255)
                                                                                                                                                 gCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAg      <  1:964679/150‑1 (MQ=255)
                                                                                                                                                    ttCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCAc   <  1:625811/150‑1 (MQ=255)
                                                                                                                                                     tCATGCCGGATGCGGCGAGAACAACGAGCGCTGAGTCAAATATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCAAAAGACGTTTCCTGCGCACGGGGAGGTTTTCAGCAcc  >  1:1269373/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATTACCACAGACATAGCACCTCAGAATTAAAAAGGGCGAGGTTATCCCTCACCCTAAACGCATAAATGCCTGATGCGCTACGCTTATCAGGCCTACATAACTCCTGCAATATTTTGGATTTGCATGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACTTAGTCAATAATCAGGATGTGAGGTTACCCCTCACCCTGTCATATCCGTAAAAACGGCATACAGAATATCACATCAGCCGCTAAAGACGTTACCTGCGCCCGGGGCGGTTTTCAGCACC  >  NC_000913/2945828‑2946126

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