Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 517,383 A→T 7.7% intergenic (+21/+42) fetB → / ← cnoX putative iron ABC exporter membrane subunit FetB/chaperedoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913517,3830AT 7.7% 234.5 / 6.5 77intergenic (+21/+42)fetB/cnoXputative iron ABC exporter membrane subunit FetB/chaperedoxin
Reads supporting (aligned to +/- strand):  ref base A (35/36);  new base T (2/4);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCT  >  NC_000913/517236‑517530
                                                                                                                                                   |                                                                                                                                                   
gATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatc                                                                                                                                                   >  1:1052291/1‑150 (MQ=255)
gATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatc                                                                                                                                                   <  2:1356821/150‑1 (MQ=255)
     tGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcgcatc                                                                                                                                              >  1:407564/1‑150 (MQ=255)
      ggTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATcc                                                                                                                                             >  2:1370145/1‑150 (MQ=255)
       gTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCg                                                                                                                                            <  1:1131305/150‑1 (MQ=255)
        tGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCgg                                                                                                                                           <  1:1393748/150‑1 (MQ=255)
            cTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTAATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCGGcatc                                                                                                                                       <  1:1394669/150‑4 (MQ=255)
              ttATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTtgt                                                                                                                                     <  1:7639/150‑1 (MQ=255)
                aTGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTgg                                                                                                                                   <  2:757547/150‑1 (MQ=255)
                 tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                                  >  2:137261/1‑150 (MQ=255)
                  gctgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGtt                                                                                                                                 <  2:265109/150‑1 (MQ=255)
                   ctgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcg                                                                                                                                                  <  1:880871/132‑1 (MQ=255)
                   ctgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcg                                                                                                                                                  >  2:880871/1‑132 (MQ=255)
                   ctgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTc                                                                                                                                >  1:163440/1‑150 (MQ=255)
                    tgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCa                                                                                                                               <  2:407564/150‑1 (MQ=255)
                    tgctCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCa                                                                                                                               <  1:295471/150‑1 (MQ=255)
                        ctcAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCAtgtg                                                                                                                           >  1:808822/1‑150 (MQ=255)
                          cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                         <  2:185953/150‑1 (MQ=255)
                          cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                         >  1:252906/1‑150 (MQ=255)
                             ccGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGAt                                                                                                                      <  1:1274469/150‑1 (MQ=255)
                                 cAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAAc                                                                                                                  <  2:808822/150‑1 (MQ=255)
                                   gCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACgg                                                                                                                <  1:667256/150‑1 (MQ=255)
                                     ttGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGaa                                                                                                              <  1:238649/150‑1 (MQ=255)
                                      tGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAat                                                                                                             <  1:1454315/150‑1 (MQ=255)
                                         cGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATc                                                                                                          >  2:205781/1‑150 (MQ=255)
                                          gACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCa                                                                                                         >  2:1342703/1‑150 (MQ=255)
                                            acATTATTGCTTGCTATTTAACCTATAGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAAt                                                                                                       <  2:677978/149‑1 (MQ=255)
                                             cATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATa                                                                                                      <  1:59946/150‑1 (MQ=255)
                                              attattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATac                                                                                                     <  1:1342703/150‑1 (MQ=255)
                                               ttattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATaca                                                                                                    <  1:205781/150‑1 (MQ=255)
                                                    gCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAAt                                                                                               >  1:1069225/1‑150 (MQ=255)
                                                     ccTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATg                                                                                              <  2:1053515/150‑1 (MQ=255)
                                                      cTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGc                                                                                             <  1:1370145/150‑1 (MQ=255)
                                                      cTGCTATTTAACCTATCCTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGc                                                                                             <  2:770351/150‑1 (MQ=255)
                                                        gCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCAt                                                                                           >  1:831115/1‑150 (MQ=255)
                                                          tATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATAc                                                                                         >  2:324981/1‑150 (MQ=255)
                                                          tATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATAc                                                                                         <  1:777900/150‑1 (MQ=255)
                                                            tttAACCTATAGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAg                                                                                       >  2:226979/1‑150 (MQ=255)
                                                              tAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCt                                                                                     >  2:555166/1‑150 (MQ=255)
                                                              tAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCt                                                                                     >  1:1435196/1‑150 (MQ=255)
                                                                 ccTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggc                                                                                  >  1:1317151/1‑150 (MQ=255)
                                                                    aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGGTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                               <  1:812294/150‑1 (MQ=255)
                                                                      cGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGa                                                                             >  1:76349/1‑150 (MQ=255)
                                                                       gTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGTGGCGAt                                                                            >  2:1259015/1‑150 (MQ=255)
                                                                       gTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGAt                                                                            >  1:484323/1‑150 (MQ=255)
                                                                           gTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACtt                                                                        >  1:1360867/1‑150 (MQ=255)
                                                                            ttttATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTc                                                                       <  1:547168/150‑1 (MQ=255)
                                                                               tataATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGAc                                                                    >  1:81321/1‑150 (MQ=255)
                                                                                ataATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGc                                                                                      <  1:327067/131‑1 (MQ=255)
                                                                                ataATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGc                                                                                      >  2:327067/1‑131 (MQ=255)
                                                                                ataATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACg                                                                   >  1:808792/1‑150 (MQ=255)
                                                                                       tcaacaccaACTGGAGGTAACGCAATTGAAAAAGAAAAGATCCCATGCCTGATCCGCACCATAACAAACAAAGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGc                                                            >  2:665951/5‑150 (MQ=255)
                                                                                        cgcCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCa                                                           <  1:1322814/150‑1 (MQ=255)
                                                                                          ccaccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATc                                                         >  2:971270/1‑150 (MQ=255)
                                                                                             ccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCAcc                                                      <  2:41042/150‑1 (MQ=255)
                                                                                               aGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCg                                                    >  2:1518102/1‑150 (MQ=255)
                                                                                                      ggtgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcag                                             <  2:1500648/150‑1 (MQ=255)
                                                                                                         gACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgc                                          <  1:324981/150‑1 (MQ=255)
                                                                                                           cGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgcag                                        >  2:376095/1‑150 (MQ=255)
                                                                                                           cGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgcag                                        >  1:350062/1‑150 (MQ=255)
                                                                                                             cAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgcagcg                                      >  1:1366925/1‑150 (MQ=255)
                                                                                                                ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                   <  1:235013/150‑1 (MQ=255)
                                                                                                                ttGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAgg                                   <  2:1052291/150‑1 (MQ=255)
                                                                                                                    agaagaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGAtct                               >  2:732681/1‑150 (MQ=255)
                                                                                                                      aagaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTcc                             <  1:137261/150‑1 (MQ=255)
                                                                                                                       agaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCt                            >  2:724503/1‑150 (MQ=255)
                                                                                                                       agaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCt                            >  1:1374517/1‑150 (MQ=255)
                                                                                                                       agaaATGATGCGATGACGGATGCTCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGTATCTCCt                            <  2:1193512/150‑1 (MQ=255)
                                                                                                                        gaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTg                           >  1:743560/1‑150 (MQ=255)
                                                                                                                        gaaATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGGTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTg                           <  1:104278/150‑1 (MQ=255)
                                                                                                                          aaTGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGa                         >  2:1330605/1‑150 (MQ=255)
                                                                                                                                   atgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAc                <  1:1340025/150‑1 (MQ=255)
                                                                                                                                      ccGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAcgcg             <  2:1486280/150‑1 (MQ=255)
                                                                                                                                      ccGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCt                            >  2:1445207/1‑135 (MQ=255)
                                                                                                                                      ccGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCt                            <  1:1445207/135‑1 (MQ=255)
                                                                                                                                      ccGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAcgcg             >  1:563644/1‑150 (MQ=255)
                                                                                                                                       cGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGt            <  1:1330605/150‑1 (MQ=255)
                                                                                                                                        ggATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTc           <  2:81321/150‑1 (MQ=255)
                                                                                                                                          aTGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATTAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTg         <  2:265053/150‑1 (MQ=255)
                                                                                                                                            gcgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAc       >  1:670295/1‑150 (MQ=255)
                                                                                                                                             cgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAcc      >  1:1329846/1‑150 (MQ=255)
                                                                                                                                             cgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAcc      >  1:1252533/1‑150 (MQ=255)
                                                                                                                                               cagcaTCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGt    <  1:176974/150‑1 (MQ=255)
                                                                                                                                                 gcatcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGCCGCCAGTGCATCCCCCGTACCCCGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCt  >  1:416453/1‑150 (MQ=255)
                                                                                                                                                 gcatcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCt  >  1:386939/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCT  >  NC_000913/517236‑517530

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: