Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 517,380 A→T 10.3% intergenic (+18/+45) 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,3800AT10.3% 192.9 / 9.9 67intergenic (+18/+45)fetB/cnoXputative iron ABC exporter membrane subunit FetB/chaperedoxin
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base T (4/3);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

TATCAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCG  >  NC_000913/517231‑517527
                                                                                                                                                     |                                                                                                                                                   
tatCAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCgca                                                                                                                                                     >  2:416644/1‑150 (MQ=255)
tatCAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCgca                                                                                                                                                     <  1:261471/150‑1 (MQ=255)
tatCAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCgca                                                                                                                                                     >  2:624479/1‑150 (MQ=255)
  tCAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCgcagc                                                                                                                                                   >  2:301053/1‑150 (MQ=255)
     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)
                              tcAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCgcag                                                                                                                                                    <  2:1420657/121‑1 (MQ=255)
                              tcAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCgcag                                                                                                                                                    >  1:1420657/1‑121 (MQ=255)
                               cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                      >  1:830433/1‑150 (MQ=255)
                               cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                      <  2:198954/150‑1 (MQ=255)
                               cAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGcc                                                                                                                      >  2:302106/1‑150 (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                                                                                                          <  1:1077820/150‑1 (MQ=255)
                                           tGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAat                                                                                                          <  1:442148/150‑1 (MQ=255)
                                           tGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAat                                                                                                          <  2:725030/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:1225782/150‑1 (MQ=255)
                                                                aTTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACa                                                                                     <  1:302106/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                    <  1:236776/141‑1 (MQ=255)
                                                                                                                                       aacgcCGGATGCGATGCTGCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAAc                    >  2:236776/4‑144 (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        <  1:711298/150‑1 (MQ=255)
                                                                                                                                             ggATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTc        <  2:1334632/150‑1 (MQ=255)
                                                                                                                                                 gcgcAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGAc    <  2:938452/150‑1 (MQ=255)
                                                                                                                                                   gcagcaTCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCg  <  1:1186949/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TATCAGATTATGGTGACCTTTATGCTGCTCTCAACCGCCAGCTTGTCGACCATTATTGCCTGCTATTTAACCTATCGTAAGTTTTATAATTCGCGCCACCAGTTGGTGGTGACGCAATTGAAGAAGAAATGATGCGATGCCGGATGCGCAGCATCGCATCCGGCGTTGTGGTTCATGTGCCGATCAACGGAATATCAATACAACAATGCATACAGCTGGCGGCGATACTTCGACGCCAGTGCATCACCCGTACCCAGCGCAGCGAGGATCTCCTGGAACGTTTTACGCGTCTGACCG  >  NC_000913/517231‑517527

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: