Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 710,173 A→G 12.3% intergenic (+57/+27) chiQ → / ← fur lipoprotein ChiQ/DNA‑binding transcriptional dual regulator Fur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913710,1730AG12.3% 235.2 / 20.0 91intergenic (+57/+27)chiQ/furlipoprotein ChiQ/DNA‑binding transcriptional dual regulator Fur
Reads supporting (aligned to +/- strand):  ref base A (46/33);  new base G (4/7);  total (50/41)
Fisher's exact test for biased strand distribution p-value = 2.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

AGTGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGACG  >  NC_000913/710025‑710319
                                                                                                                                                    |                                                                                                                                                  
agTGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAg                                                                                                                                                   >  1:131362/1‑150 (MQ=255)
 gtgtCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAgg                                                                                                                                                  <  2:115697/150‑1 (MQ=255)
      tGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                             <  2:1363334/150‑1 (MQ=255)
      tGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGc                                                                                                                                             <  1:292729/150‑1 (MQ=255)
       gCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCt                                                                                                                                            <  2:1169139/150‑1 (MQ=255)
       gCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCt                                                                                                                                            <  2:642642/150‑1 (MQ=255)
        cGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCtt                                                                                                                                           <  2:257375/150‑1 (MQ=255)
           gCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttt                                                                                                                                        <  1:272303/150‑2 (MQ=255)
           gCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtc                                                                                                                                        >  1:1485245/1‑150 (MQ=255)
            cAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtct                                                                                                                                       <  1:1494272/150‑1 (MQ=255)
            cAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTtct                                                                                                                                       >  2:508343/1‑150 (MQ=255)
              aCGCAAACCGGAAATGATCAGGCGGGGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAAAAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttttat                                                                                                                                     <  2:667106/150‑5 (MQ=255)
               cGCAAACCGGAAATTATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGt                                                                                                                                    <  2:454652/150‑1 (MQ=255)
               cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAGGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGt                                                                                                                                    >  2:301236/1‑150 (MQ=255)
               cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCATTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttttatg                                                                                                                                    <  2:254570/150‑6 (MQ=255)
               cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCGGGTTGGCttttttatg                                                                                                                                    <  1:408724/150‑6 (MQ=255)
               cGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGt                                                                                                                                    >  1:873687/1‑150 (MQ=255)
                 cAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTc                                                                                                                                  <  2:579885/150‑1 (MQ=255)
                 cAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTc                                                                                                                                  <  2:468793/150‑1 (MQ=255)
                  aaaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCa                                                                                                                                 >  1:573634/1‑150 (MQ=255)
                   aaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAg                                                                                                                                >  1:1359647/1‑150 (MQ=255)
                   aaCCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAg                                                                                                                                <  1:623063/150‑1 (MQ=255)
                     ccGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggc                                                                                                                              >  1:1196326/1‑150 (MQ=255)
                       ggAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctg                                                                                                                            >  1:905511/1‑150 (MQ=255)
                        gAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAggctgg                                                                                                                           <  2:163675/150‑1 (MQ=255)
                             gATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAt                                                                                                                      >  1:537314/1‑150 (MQ=255)
                             gATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTAt                                                                                                                      >  2:675453/1‑150 (MQ=255)
                                cAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTg                                                                                                                   <  2:1365397/150‑1 (MQ=255)
                                  ggCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGcc                                                                                                                 >  1:1010487/1‑150 (MQ=255)
                                    cGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCtt                                                                                                               >  1:67082/1‑150 (MQ=255)
                                      gTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCg                                                                                                             <  2:442683/150‑1 (MQ=255)
                                          aaGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTgcg                                                                                                         >  2:426816/1‑150 (MQ=255)
                                           aGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCTGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTgcgc                                                                                                        <  2:852288/150‑1 (MQ=255)
                                                  tgtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAATTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATATATAATCCAAGCCGCAGGTTGGCTTTTCTCGTTCAGGCTGTCTTATTTGCCTTCGTGCGCATGTTCa                                                                                                 <  2:1339588/150‑1 (MQ=255)
                                                  tgtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCa                                                                                                 <  1:450541/150‑1 (MQ=255)
                                                   gtgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAt                                                                                                >  2:321963/1‑150 (MQ=255)
                                                    tgATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCAta                                                                                               >  1:584627/1‑149 (MQ=255)
                                                     gATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCt                                                                                              >  1:840942/1‑150 (MQ=255)
                                                         aaGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgc                                                                                          >  2:411376/1‑150 (MQ=255)
                                                          aGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgcg                                                                                         <  2:180595/150‑1 (MQ=255)
                                                          aGTCTGCCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTATTCTTATTTCCATTGCATAAAAAAGCCAACCCGCACGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTcgcg                                                                                         >  1:1179325/1‑150 (MQ=255)
                                                           ttCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                                                                        <  2:1485245/149‑1 (MQ=255)
                                                           gTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCGAATTCGCTGTTTGGCTTTTCTCGTTCAGGCTGTCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                                                                        <  2:637940/150‑1 (MQ=255)
                                                           gTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACGCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCgg                                                                                        <  1:605047/150‑1 (MQ=255)
                                                            tCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAACAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGc                                                                                       <  2:241472/150‑1 (MQ=255)
                                                            tCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGc                                                                                       <  2:1280177/150‑1 (MQ=255)
                                                            tCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGc                                                                                       <  2:905511/150‑1 (MQ=255)
                                                                gCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATc                                                                                   >  2:1241372/1‑150 (MQ=255)
                                                                  aGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGc                                                                                 >  2:575160/1‑150 (MQ=255)
                                                                  aGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGc                                                                                 >  1:190235/1‑150 (MQ=255)
                                                                   ggAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGcc                                                                                <  2:67082/150‑1 (MQ=255)
                                                                        tACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTGGGCTGTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCgg                                                                           >  2:385529/1‑150 (MQ=255)
                                                                                aataaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGa                                                                   >  2:155303/1‑150 (MQ=255)
                                                                                  taaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAcc                                                                 >  1:279782/1‑150 (MQ=255)
                                                                                  taaCAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGAcc                                                                 >  2:450314/1‑150 (MQ=255)
                                                                                       aaTAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTaaa                                                            <  2:1032298/150‑1 (MQ=255)
                                                                                             tGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga                                                      <  1:265157/150‑1 (MQ=255)
                                                                                             tGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATaga                                                      <  1:301236/150‑1 (MQ=255)
                                                                                              gagaGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATagag                                                     >  2:1348028/1‑150 (MQ=255)
                                                                                                  gCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtg                                                 >  1:736258/1‑150 (MQ=255)
                                                                                                    tGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACtgtg                                               >  2:780728/1‑150 (MQ=255)
                                                                                                         ctctCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAg                                          >  1:115631/1‑150 (MQ=255)
                                                                                                         ctctCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAg                                          >  2:201461/1‑150 (MQ=255)
                                                                                                          tctcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGt                                         >  1:1083195/1‑150 (MQ=255)
                                                                                                           ctcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTc                                        >  2:962407/1‑150 (MQ=255)
                                                                                                           ctcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTc                                        >  2:567253/1‑150 (MQ=255)
                                                                                                            tcGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCa                                       <  1:938248/150‑1 (MQ=255)
                                                                                                                ttttCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCg                                   <  1:321963/150‑1 (MQ=255)
                                                                                                                ttttCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCg                                   <  2:515072/150‑1 (MQ=255)
                                                                                                                   tCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTCAGTCAGGCGAAt                                >  1:2778/1‑150 (MQ=255)
                                                                                                                    cTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATg                               <  1:352720/150‑1 (MQ=255)
                                                                                                                    cTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATg                               <  2:1010487/150‑1 (MQ=255)
                                                                                                                      tATTTCCCTTGCATAAAAAAGCCAACCCTCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGcc                             <  2:616150/150‑1 (MQ=255)
                                                                                                                           cccTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGtt                        >  2:1354589/1‑150 (MQ=255)
                                                                                                                             cTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGtttt                      >  1:841035/1‑150 (MQ=255)
                                                                                                                               tGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGc                    >  2:1232856/1‑150 (MQ=255)
                                                                                                                                gCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCg                   >  1:324683/1‑150 (MQ=255)
                                                                                                                                 cATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCgg                  >  1:802629/1‑150 (MQ=255)
                                                                                                                                 cATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCgg                  <  1:450314/150‑1 (MQ=255)
                                                                                                                                   tAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCa                >  1:1466297/1‑150 (MQ=255)
                                                                                                                                   tAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCa                <  1:567253/150‑1 (MQ=255)
                                                                                                                                    agaaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCaa               >  1:327625/3‑150 (MQ=255)
                                                                                                                                    aaaaaaGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCaa               >  1:1272868/1‑150 (MQ=255)
                                                                                                                                      aaaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAAtt             >  2:1266140/1‑150 (MQ=255)
                                                                                                                                       aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATg                               <  1:1440583/131‑1 (MQ=255)
                                                                                                                                       aaaGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATg                               >  2:1440583/1‑131 (MQ=255)
                                                                                                                                         aGCCAACCTGCGGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCa          >  1:849267/1‑150 (MQ=255)
                                                                                                                                           ccAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACg        >  1:1204674/1‑150 (MQ=255)
                                                                                                                                           ccAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGAAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACg        >  1:638209/1‑150 (MQ=255)
                                                                                                                                            cAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGc       >  2:989897/1‑150 (MQ=255)
                                                                                                                                                 cGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGAcg  >  1:1033590/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGTGTCTGCGTGCAACGCAAACCGGAAATGATCAGGCGGTGAAAGCCGATTGTGATAAAGTCTGGCAGGAAATACGCAGTAATAACAAATAAGTGAGAGCTGTAACTCTCGCTTTTCTTATTTCCCTTGCATAAAAAAGCCAACCCGCAGGTTGGCTTTTCTCGTTCAGGCTGGCTTATTTGCCTTCGTGCGCATGTTCATCTTCGCGGCAATCGCCTTCGGCACAGTGACCGTAAAGATAGAGACTGTGGTTAGTCAGGCGAATGCCATGTTTTGCGGCAATTTCACGCTGACG  >  NC_000913/710025‑710319

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