Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 517,384 T→G 10.9% intergenic (+22/+41) 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,3840TG10.9% 177.2 / 13.5 64intergenic (+22/+41)fetB/cnoXputative iron ABC exporter membrane subunit FetB/chaperedoxin
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base G (4/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTGCG  >  NC_000913/517236‑517533
                                                                                                                                                    |                                                                                                                                                     
gATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatc                                                                                                                                                      <  2:1224418/150‑1 (MQ=255)
 aTTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcg                                                                                                                                                     <  2:1140903/150‑1 (MQ=255)
 aTTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcg                                                                                                                                                     <  2:525280/150‑1 (MQ=255)
 aTTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAAACTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAgcatcg                                                                                                                                                     <  2:530351/150‑1 (MQ=255)
        tGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGATGCTGCGCATCCgg                                                                                                                                              <  1:1087868/150‑1 (MQ=255)
         gACCTATATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGACGCCGGATGCGATGCTGCGCATCCGGc                                                                                                                                             >  2:1325594/1‑150 (MQ=255)
          cccTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGCTGGTGGTGACGCAATAGAAGAAGAAATGATGCGATGCCGGAGGCGCAGCATCGCATCCGGCg                                                                                                                                            <  1:895887/149‑1 (MQ=255)
            cTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGACGCCGGATGCGATGCTGCGCATCCGGcatc                                                                                                                                          <  1:1325594/150‑4 (MQ=255)
             tttATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTTGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTtg                                                                                                                                         <  2:1233293/150‑1 (MQ=255)
                 tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                                     >  1:615160/1‑150 (MQ=255)
                 tgctgcTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGt                                                                                                                                     >  1:725030/1‑150 (MQ=255)
                      ctctCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCAtg                                                                                                                                >  2:1099918/1‑150 (MQ=255)
                         tcAACCGCCAGGTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGc                                                                                                                             <  2:1082629/150‑1 (MQ=255)
                          cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                            >  2:302106/1‑150 (MQ=255)
                          cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                            >  1:830433/1‑150 (MQ=255)
                          cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                            <  2:198954/150‑1 (MQ=255)
                             ccGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGAt                                                                                                                         >  2:273768/1‑150 (MQ=255)
                               gCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCa                                                                                                                       >  2:365461/1‑150 (MQ=255)
                               gCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCa                                                                                                                       >  1:1266515/1‑150 (MQ=255)
                                  aGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACg                                                                                                                    >  2:934774/1‑150 (MQ=255)
                                   gCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACgg                                                                                                                   >  2:1118078/1‑150 (MQ=255)
                                      tGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAat                                                                                                                <  2:725030/150‑1 (MQ=255)
                                      tGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAat                                                                                                                <  1:1077820/150‑1 (MQ=255)
                                      tGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAat                                                                                                                <  1:442148/150‑1 (MQ=255)
                                          gACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCa                                                                                                            >  1:415035/1‑150 (MQ=255)
                                              attattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATac                                                                                                        >  1:633200/1‑150 (MQ=255)
                                                tattGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATacaa                                                                                                      >  1:959440/1‑150 (MQ=255)
                                                 attGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATacaac                                                                                                     <  1:16989/150‑1 (MQ=255)
                                                   tGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATAcaacaa                                                                                                   <  1:365461/150‑1 (MQ=255)
                                                    gCCTGCTATTTAAACTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGAAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAAt                                                                                                  >  2:1186949/1‑150 (MQ=255)
                                                      cTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATg                                                                                                 >  1:936823/1‑149 (MQ=255)
                                                      cTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATg                                                                                                 <  2:936823/149‑1 (MQ=255)
                                                      cTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGc                                                                                                >  1:1334632/1‑150 (MQ=255)
                                                          tATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATAc                                                                                            <  2:315455/150‑1 (MQ=255)
                                                           aTTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACa                                                                                           <  1:302106/150‑1 (MQ=255)
                                                           aTTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACa                                                                                           <  1:1225782/150‑1 (MQ=255)
                                                               aaCCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTg                                                                                       >  1:999479/1‑150 (MQ=255)
                                                                  cTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcg                                                                                    >  1:880364/1‑150 (MQ=255)
                                                                    aTCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTggcggc                                                                                  <  1:496637/150‑1 (MQ=255)
                                                                       gTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGAt                                                                               <  1:533529/150‑1 (MQ=255)
                                                                        tAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATa                                                                              <  2:615160/150‑1 (MQ=255)
                                                                             tttATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCg                                                                         >  1:877226/1‑150 (MQ=255)
                                                                              ttATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCAt                                                                                                                                 >  2:1169494/1‑93 (MQ=255)
                                                                              ttATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCAt                                                                                                                                 <  1:1169494/93‑1 (MQ=255)
                                                                                 taATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGc                                                                     >  2:333730/1‑150 (MQ=255)
                                                                                      cgcgCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTg                                                                >  2:395923/1‑150 (MQ=255)
                                                                                        cgcCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCa                                                              <  1:808242/150‑1 (MQ=255)
                                                                                           caccaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCa                                                           >  1:941022/1‑150 (MQ=255)
                                                                                            accaGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCAc                                                          >  1:71531/1‑150 (MQ=255)
                                                                                                       gtgACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagc                                               <  1:273768/150‑1 (MQ=255)
                                                                                                         gACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgc                                             <  1:236267/150‑1 (MQ=255)
                                                                                                           cGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgcag                                           >  1:703235/1‑150 (MQ=255)
                                                                                                             cAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCcagcgcagcg                                         >  1:49617/1‑150 (MQ=255)
                                                                                                                 tGAAGAAGAACTGCTGCGATGCCGGAGGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATTCATACTGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCTCAGCGAggt                                     <  2:633200/150‑2 (MQ=255)
                                                                                                                         aaaTGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTgg                             <  1:1026721/150‑1 (MQ=255)
                                                                                                                              ccacaacgcCGGATGCGATGCTGCGCATCCGGCATTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGt                        >  2:199903/8‑150 (MQ=255)
                                                                                                                                  gatgcCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCACTACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTa                    <  2:135874/150‑1 (MQ=255)
                                                                                                                                  aacgcCGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAAc                          >  2:236776/4‑144 (MQ=255)
                                                                                                                                  aacgcCGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAAc                          <  1:236776/141‑1 (MQ=255)
                                                                                                                                    tgcCGGATGCGCAGCATCGCATCCGGCGTTGTCCTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCTATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAcg                  <  2:741591/150‑1 (MQ=255)
                                                                                                                                     gcCGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATTTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATTCTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAcgc                 >  2:1107976/1‑150 (MQ=255)
                                                                                                                                      ccGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTAcgcg                >  1:336731/1‑150 (MQ=255)
                                                                                                                                        ggATGCGCAGCATCTCATCCGGCGTTTTTGTTCATGTGCCGATCAAAGGAATATCAATACATCAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCTCAGCGAGGATCTCCTGGAACGTTTTACGCGTc              <  2:71531/150‑1 (MQ=255)
                                                                                                                                        ggATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTc              <  2:49617/150‑1 (MQ=255)
                                                                                                                                        ggATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTc              <  2:1334632/150‑1 (MQ=255)
                                                                                                                                        ggATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTc              <  1:711298/150‑1 (MQ=255)
                                                                                                                                            gcgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAc          <  2:938452/150‑1 (MQ=255)
                                                                                                                                              gcagcaTCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCg        <  1:1186949/150‑1 (MQ=255)
                                                                                                                                                  catcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTg    <  2:880364/150‑1 (MQ=255)
                                                                                                                                                   atcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTgc   >  1:895005/1‑150 (MQ=255)
                                                                                                                                                    tcgcatcCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATCCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTgcg  <  1:166482/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCGTCTGCG  >  NC_000913/517236‑517533

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: