Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,552,271 T→C 10.7% intergenic (+280/+127) 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,2710TC10.7% 239.2 / 16.4 90intergenic (+280/+127)fdnI/yddMformate dehydrogenase N subunit gamma/putative DNA‑binding transcriptional regulator YddM
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base C (7/2);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

TAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTA  >  NC_000913/1552122‑1552420
                                                                                                                                                     |                                                                                                                                                       
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                                                                                           >  1:403751/1‑150 (MQ=40)
                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                           >  2:1364911/1‑150 (MQ=40)
                                                            catcaGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTACCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGc                                                                                           >  2:1110895/1‑150 (MQ=255)
                                                               caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                        >  2:256519/1‑150 (MQ=38)
                                                               caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                        >  2:891233/1‑150 (MQ=38)
                                                               caGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCATTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                        >  2:1656982/1‑150 (MQ=38)
                                                               caGTTTCAAATGGAGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCTTGAACGCCTTAACCGGCATGA‑‑AAAACCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTa                                                                                        >  2:1656984/1‑150 (MQ=255)
                                                                 gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt                                                                                         <  1:360095/147‑1 (MQ=35)
                                                                 gTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGt                                                                                         >  2:360095/1‑147 (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:1336621/1‑150 (MQ=255)
                                                                                                                           tCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCAtt                            >  1:1263604/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              >  1:189823/1‑150 (MQ=255)
                                                                                                                                         gcgTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTaa              >  2:1551534/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                      >  1:1476997/1‑139 (MQ=255)
                                                                                                                                            tGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTa                      <  2:1476997/139‑1 (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)
                                                                                                                                                   cTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg    <  1:1289111/150‑1 (MQ=255)
                                                                                                                                                   cTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg    <  2:41284/150‑1 (MQ=255)
                                                                                                                                                   cTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCg    <  1:899209/150‑1 (MQ=255)
                                                                                                                                                   cTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGAAGGCg    >  1:1091792/1‑150 (MQ=255)
                                                                                                                                                     tATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTa  >  1:474814/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTTCAAATGGCGCTGTAAAAGGCGTCATTTTCATATTGTAGACAACGTAGGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCATGA‑‑AAACCCTTCAAATCCAATAGATTGCAGTGAACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCAATGTTGCGTTTGTCATCAGTTCTAAATGGCGCTTTATAAAGTGCCATTTTTTTATTGCGTAACCAGACGGCGTA  >  NC_000913/1552122‑1552420

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