Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,264 A→C 10.0% intergenic (+273/+134) fdnI → / ← yddM formate dehydrogenase N subunit gamma/putative DNA‑binding transcriptional regulator YddM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,552,2640AC10.0% 243.0 / 19.1 90intergenic (+273/+134)fdnI/yddMformate dehydrogenase N subunit gamma/putative DNA‑binding transcriptional regulator YddM
Reads supporting (aligned to +/- strand):  ref base A (44/37);  new base C (7/2);  total (51/39)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGA  >  NC_000913/1552118‑1552413
                                                                                                                                                  |                                                                                                                                                       
ccAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTATCTCATCAGGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACg                                                                                                                                                      <  2:1000263/150‑1 (MQ=255)
 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGATTTGTTCATGCCGGATGCGGCGTGAACGc                                                                                                                                                     >  1:1095988/1‑150 (MQ=255)
 cAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAATCATTTTAATATTGTATACAACGTAGGCTTTGTTCATGCCGGATCCGGCTTGAAAGc                                                                                                                                                     >  1:436148/1‑150 (MQ=21)
  aaTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGcc                                                                                                                                                    <  2:3431/150‑1 (MQ=255)
    tAGATTGCAGTGAACGTGTTGGACTTATAACCGTAGCGCATGAGGCAATGTTTCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCtt                                                                                                                                                  <  2:821155/150‑1 (MQ=21)
          gCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCgg                                                                                                                                            <  1:793010/150‑1 (MQ=255)
           cactGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTGCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                           <  2:548135/147‑1 (MQ=255)
           cAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACTACGTAGGCTTTTTTCATTCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                                           <  2:1519933/150‑1 (MQ=255)
            aGTGAACGTGTAGGCCTGATAAGCGTATCGCATCAGGCAATGTTGTGTTTTTTATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCa                                                                                                                                          <  2:756241/150‑1 (MQ=255)
              ttaaCGTCTAGGACTCATAACGTTACCGCATCAGTCAATGTTTCTTTTTTATTCAGTTTAAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                                                                                                        <  2:94871/148‑1 (MQ=9)
              tgcaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCTGATGCGGCGTGAACGCCTTATCCGGAATg                                                                                                                                        <  2:84601/147‑1 (MQ=255)
              tGAACGTGTAGTCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTTTCATCAGTTTAAAATGGCTCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTTTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATg                                                                                                                                        <  2:1387178/150‑1 (MQ=21)
                aaCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGACATTTTAATATTGGAGCAAACGTAGACATTGTTCATGCCGGATGCGGAGGGAACGGCTTAAGCTGCATGa‑‑a                                                                                                                                    >  2:767213/1‑150 (MQ=18)
                   gtgtAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAAc                                                                                                                                 <  2:428464/150‑1 (MQ=255)
                      tAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCt                                                                                                                              >  2:855834/1‑150 (MQ=255)
                       aGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTAATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGGGTGAACGGCTTATACGGCATGA‑‑AAACCCtt                                                                                                                             >  2:132678/1‑150 (MQ=21)
                        ggCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTc                                                                                                                            <  2:128593/150‑1 (MQ=255)
                          ccTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTAAAATGGCGCTGTAAAAGGAGTCATTTTGTTATTGTATACAACGTAGGCTTTGTTAATGCCGGATGCGGCGTGAACGACTTATAAGGAATGA‑‑AAACCCTTCaa                                                                                                                          >  2:1255273/1‑150 (MQ=21)
                          ccTGATAACCGTAGCGCATCAGGCAATGTGGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCaa                                                                                                                          <  2:168516/150‑1 (MQ=255)
                           cTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCaaa                                                                                                                         >  1:1610179/1‑150 (MQ=255)
                            tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAAt                                                                                                                        <  2:1273770/150‑1 (MQ=255)
                            tGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTGTAAAAGGAGACATTTTAATATTGAATAAAAAGTAGGCTTTGTTCATACCCGATCCGGAGTGAACGCCTTATACTGCATGA‑‑AAACACTTCaaaa                                                                                                                        >  2:1338806/1‑149 (MQ=21)
                             gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAAcc                                                                                                                                >  2:72725/1‑141 (MQ=255)
                             gATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAAcc                                                                                                                                <  1:72725/141‑1 (MQ=255)
                              aTAAGCGTAGCGCATCAGGCAATGTTGCGTGTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATcc                                                                                                                      <  2:769128/150‑1 (MQ=255)
                                   cgtatatcATCAGGGAATGTTGATTTTGTCATCAGTTTGAAATGGCTCTGTAAAAGGCGTCATTTTCATATTGTAGTCAAGGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGTAATGA‑‑AAACCCTTCAAATCCAATAg                                                                                                                 <  2:1095988/143‑1 (MQ=17)
                                   cGTAGCTCATCATGCAATGTTGCGTTTTTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAg                                                                                                                 <  2:914816/150‑1 (MQ=255)
                                          cATCAGGCAATGTTGCCTTTGTCATCAGTTTCAAATGGCGCTGTAAAATTTTTCATTTTAATATTGTAGACAACTTAGGCTTTGTTCATGACTGATTAGGCGTTAACGCCTTTTCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAg                                                                                                          <  2:592217/150‑1 (MQ=17)
                                            tCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTg                                                                                                        <  1:855834/150‑1 (MQ=255)
                                               ggCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTGATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCAATATACGGCATGA‑‑AAACCCTTAAAATCCAATAGATTGCAGTGAAc                                                                                                     >  1:1370665/1‑150 (MQ=17)
                                                             tGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATa                                                                                       <  2:237990/150‑1 (MQ=255)
                                                               tcatcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAg                                                                                     >  2:686734/1‑150 (MQ=40)
                                                                catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                    >  2:1364911/1‑150 (MQ=40)
                                                                catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                    >  1:403751/1‑150 (MQ=40)
                                                                catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTACCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                    >  2:1110895/1‑150 (MQ=255)
                                                                   caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                 >  2:891233/1‑150 (MQ=38)
                                                                   caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                 >  2:256519/1‑150 (MQ=38)
                                                                   caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCATTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                 >  2:1656982/1‑150 (MQ=38)
                                                                   caGTTTCAAATGGAGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTAACCGGCATGA‑‑AAAACCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                 >  2:1656984/1‑150 (MQ=255)
                                                                     gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt                                                                                  >  2:360095/1‑147 (MQ=35)
                                                                     gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt                                                                                  <  1:360095/147‑1 (MQ=35)
                                                                         cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                           >  2:337398/1‑150 (MQ=32)
                                                                         cAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCAt                                                                           >  2:1110026/1‑150 (MQ=32)
                                                                           aaTGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGGAGCGCCTCa                                                                         >  2:1055424/1‑150 (MQ=17)
                                                                               gcgcTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                                     >  1:41284/1‑150 (MQ=32)
                                                                                cgcTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGAAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                                    >  2:224984/1‑150 (MQ=17)
                                                                                    gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGa‑‑a                                                                                                                                    <  1:643402/82‑1 (MQ=33)
                                                                                    gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                                >  2:1201375/1‑150 (MQ=32)
                                                                                    gTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGtt                                                                >  2:857891/1‑150 (MQ=32)
                                                                                     tAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTg                                                               >  1:202420/1‑150 (MQ=32)
                                                                                         aGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGtt                                                           <  1:1110026/150‑1 (MQ=32)
                                                                                                      atatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCt                                              >  1:2954/1‑150 (MQ=31)
                                                                                                       tatTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTa                                             >  1:1113155/1‑150 (MQ=32)
                                                                                                        atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaa                                            >  2:389975/1‑150 (MQ=32)
                                                                                                        atTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTaa                                            >  2:142999/1‑150 (MQ=32)
                                                                                                             aGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGgc                                       >  2:1550894/1‑150 (MQ=34)
                                                                                                                 aaCGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCtt                                   >  2:467500/1‑150 (MQ=35)
                                                                                                                   cGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                                 >  2:807201/1‑150 (MQ=37)
                                                                                                                   attAGGCTTTTTTCATTCAGTATCGTTCGTTAAGGCGTTATGCGGATTTA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTta                                 <  1:219377/148‑1 (MQ=11)
                                                                                                                         cTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGt                           <  1:224984/150‑1 (MQ=255)
                                                                                                                               tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                     >  1:1263604/1‑150 (MQ=255)
                                                                                                                               tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                     >  1:1336621/1‑150 (MQ=255)
                                                                                                                               tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAtt                     <  1:584202/150‑1 (MQ=255)
                                                                                                                                cATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttt                    >  1:1229766/1‑150 (MQ=14)
                                                                                                                                cATGCCCGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCACTTTATAAAGCGCCAttt                    <  1:847702/150‑1 (MQ=255)
                                                                                                                                 aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                     >  1:1602438/1‑150 (MQ=14)
                                                                                                                                 aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑ACACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                   >  1:692295/1‑150 (MQ=14)
                                                                                                                                 aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttt                   >  2:1389195/1‑150 (MQ=14)
                                                                                                                                  tGCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttt                  >  2:234942/1‑150 (MQ=14)
                                                                                                                                   gCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                 <  2:1215816/150‑1 (MQ=255)
                                                                                                                                   gCCGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑ACACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                 <  2:692295/150‑1 (MQ=17)
                                                                                                                                   gCCGGACGCGGCTTGATCTCCTAATCAGGCCTGA‑‑CAACCCTTCATTTTCACTTGACTGCAGTGCAGCCGTCGGCCCTACCTCCCTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtttttt                 <  1:1069937/150‑1 (MQ=12)
                                                                                                                                    ccGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTAAAAATGGCGCTTTATAAAGAGACAttttaat                >  2:899209/1‑147 (MQ=11)
                                                                                                                                    ccGGATAAGGCGTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAttttttt                >  1:728877/1‑150 (MQ=17)
                                                                                                                                        aTGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTg            <  2:64660/150‑1 (MQ=255)
                                                                                                                                          gcggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg          <  1:132678/150‑1 (MQ=255)
                                                                                                                                          gcggcgGGAACGCCTTATCCGGCATGA‑‑AAACCCTTCCACTCCACTAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg          <  1:767213/150‑1 (MQ=255)
                                                                                                                                          aaggcgTTCACGCCGCATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCg          >  1:225541/3‑150 (MQ=255)
                                                                                                                                            ggcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCACTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTa        >  2:1659685/1‑150 (MQ=255)
                                                                                                                                             gcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa       >  2:1551534/1‑150 (MQ=255)
                                                                                                                                             gcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa       >  1:189823/1‑150 (MQ=255)
                                                                                                                                             gcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa       >  1:1558661/1‑150 (MQ=255)
                                                                                                                                              cgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc      >  2:427420/1‑150 (MQ=255)
                                                                                                                                              cgTGAACGCCTTATCCCGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAc      <  1:1255273/150‑1 (MQ=255)
                                                                                                                                               gTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAcc     <  1:807201/150‑1 (MQ=255)
                                                                                                                                               gTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAAcc     <  1:1550894/150‑1 (MQ=255)
                                                                                                                                                tGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa               <  2:1476997/139‑1 (MQ=255)
                                                                                                                                                tGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa               >  1:1476997/1‑139 (MQ=255)
                                                                                                                                                tGAACGCCTAATCCCTCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGAAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCa    <  1:968908/150‑1 (MQ=255)
                                                                                                                                                 gAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAg   >  1:1102480/1‑150 (MQ=255)
                                                                                                                                                  aaCGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGa  <  1:337398/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                       
CCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGA  >  NC_000913/1552118‑1552413

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