Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,069,911 C→A 6.5% G148G (GGG→GGT)  rutF ← flavin:NADH reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,069,9110CA 6.5% 301.0 / 6.0 88G148G (GGG→GGT) rutFflavin:NADH reductase
Reads supporting (aligned to +/- strand):  ref base C (39/43);  new base A (2/4);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGC  >  NC_000913/1069768‑1070059
                                                                                                                                               |                                                                                                                                                    
aTCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTa                                                                                                                                                 <  1:206978/149‑1 (MQ=255)
 tCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAg                                                                                                                                                >  2:122198/1‑149 (MQ=255)
 tCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAg                                                                                                                                                <  2:332079/149‑1 (MQ=255)
    ggCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggt                                                                                                                                             <  4:187669/149‑1 (MQ=255)
    ggCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggt                                                                                                                                             >  2:491689/1‑149 (MQ=255)
     gcgcGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtg                                                                                                                                            >  1:364916/1‑149 (MQ=255)
       gcgACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATGAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtg                                                                                                                                          >  4:399085/1‑149 (MQ=255)
       gcgACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtg                                                                                                                                          <  1:628758/149‑1 (MQ=255)
          aCCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtg                                                                                                                                       >  1:386389/1‑149 (MQ=255)
           ccACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtgt                                                                                                                                      <  4:105837/149‑1 (MQ=255)
           ccACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGAGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggtgt                                                                                                                                      >  2:378433/1‑149 (MQ=255)
            cACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGgtgtggtgtg                                                                                                                                     <  4:221713/149‑1 (MQ=255)
             aCGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCACACCCTACGggctggtgtgg                                                                                                                                    <  4:137926/149‑2 (MQ=255)
              cgccgcTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCACACCCTACGggctggtgtggt                                                                                                                                   <  1:610570/149‑3 (MQ=255)
                    tctctGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGaa                                                                                                                             <  1:491689/149‑1 (MQ=255)
                        tGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGc                                                                                                                         >  4:431267/1‑149 (MQ=255)
                          tAGAGGTCGATTTTAACTGCCAGTGAGGAAAACAGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCtt                                                                                                                       <  2:425128/149‑1 (MQ=255)
                           agagGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                                                                                                      <  2:474101/149‑1 (MQ=255)
                            gagGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                                                                                                      <  3:157499/148‑1 (MQ=255)
                            gagGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                                                                                                      >  4:157499/1‑148 (MQ=255)
                                 cGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGgc                                                                                                                <  2:206336/149‑1 (MQ=255)
                                  gATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGgcg                                                                                                               >  3:132729/1‑149 (MQ=255)
                                  gATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGgcg                                                                                                               >  3:14934/1‑149 (MQ=255)
                                    ttttAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgg                                                                                                                                               <  3:300308/115‑1 (MQ=255)
                                    ttttAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAgg                                                                                                                                               >  4:300308/1‑115 (MQ=255)
                                    ttttAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCACACCCTACGggctggtgtg                                                                                                                                     >  1:489431/1‑125 (MQ=255)
                                    ttttAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCACACCCTACGggctggtgtg                                                                                                                                     <  2:489431/125‑1 (MQ=255)
                                       tAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCaaaa                                                                                                          <  2:556675/149‑1 (MQ=255)
                                         aCTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACa                                                                                                        <  4:132729/149‑1 (MQ=255)
                                            gCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGaa                                                                                                     <  2:364916/149‑1 (MQ=255)
                                            gCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGaa                                                                                                     <  3:157213/149‑1 (MQ=255)
                                                gTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTc                                                                                                 <  4:108816/149‑1 (MQ=255)
                                                gTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTc                                                                                                 >  4:377461/1‑149 (MQ=255)
                                                   aGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTg                                                                                              <  3:431267/149‑1 (MQ=255)
                                                    ggAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTgg                                                                                             <  3:130674/149‑1 (MQ=255)
                                                    ggAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTgg                                                                                             >  3:7906/1‑149 (MQ=255)
                                                    ggAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTgg                                                                                             <  3:399085/149‑1 (MQ=255)
                                                     gAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTggg                                                                                            >  2:536641/1‑149 (MQ=255)
                                                      aaaaCCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGt                                                                                           >  2:362116/1‑149 (MQ=255)
                                                        aaCCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGc                                                                                         <  2:386389/149‑1 (MQ=255)
                                                            gAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGAc                                                                                     >  1:389574/1‑149 (MQ=255)
                                                                    ccATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcacca                                                                             >  3:477385/1‑149 (MQ=255)
                                                                     cATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcaccac                                                                            >  2:535897/1‑149 (MQ=255)
                                                                     cATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTcaccac                                                                            >  3:65837/1‑149 (MQ=255)
                                                                       tCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCt                                                                          <  4:65837/149‑1 (MQ=255)
                                                                        cAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTg                                                                         >  3:250735/1‑149 (MQ=255)
                                                                              tctcCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTTCCGACGCTCACCACCTGGCTGAt                                                                   <  1:234746/149‑1 (MQ=255)
                                                                                 ccTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggt                                                                                                                                        <  4:109970/77‑1 (MQ=255)
                                                                                 ccTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGggtgtggt                                                                                                                                        >  3:109970/1‑77 (MQ=255)
                                                                                  cTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGaa                                                                                                     >  3:87165/1‑111 (MQ=255)
                                                                                  cTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGaa                                                                                                     <  4:87165/111‑1 (MQ=255)
                                                                                   ttAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGc                                                              >  4:432995/1‑149 (MQ=255)
                                                                                   ttAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCACACCCTACGggctggtgtg                                                                                                                                     <  1:189121/78‑1 (MQ=25)
                                                                                   ttAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCACACCCTACGggctggtgtg                                                                                                                                     >  2:189121/1‑78 (MQ=25)
                                                                                     aaGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAg                                                            <  3:377461/149‑1 (MQ=255)
                                                                                       ggaggTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAAt                                                                                                                            <  3:289493/83‑1 (MQ=255)
                                                                                       ggaggTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAAt                                                                                                                            >  4:289493/1‑83 (MQ=255)
                                                                                           gTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAAc                                                      <  4:224061/149‑1 (MQ=255)
                                                                                              aaCAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACa                                                                                                        <  2:331240/96‑1 (MQ=255)
                                                                                              aaCAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACa                                                                                                        >  1:331240/1‑96 (MQ=255)
                                                                                                  aGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCa                                               >  1:389253/1‑149 (MQ=255)
                                                                                                  aGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCa                                               >  4:311300/1‑149 (MQ=255)
                                                                                                    cAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCaga                                             <  3:551345/149‑1 (MQ=255)
                                                                                                     aGGGTGCATCAGCGCGTGATAACTGCGATCCAACCACACCAGCCCGTAGGGTGTGGTGTTACGCTGACTCCCTTCGAAGGCGCAAAACAGGACGTCGTGGGTGCAGACGCCCACCACCTGGCTGATACGGCAGTCAAAGTAAACCagag                                            <  4:500697/149‑1 (MQ=255)
                                                                                                     aGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCagag                                            >  2:538552/1‑149 (MQ=255)
                                                                                                      gggCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGc                                           >  4:215537/1‑149 (MQ=255)
                                                                                                          gcATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCctct                                       <  3:584148/149‑1 (MQ=255)
                                                                                                           cATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCtt                                      <  1:389160/149‑1 (MQ=255)
                                                                                                             tCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTcc                                    <  1:21361/149‑1 (MQ=255)
                                                                                                             tCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTcc                                    >  1:124962/1‑149 (MQ=255)
                                                                                                                            tGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCg                     <  4:7906/149‑1 (MQ=255)
                                                                                                                             gCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCgg                    >  2:685880/1‑149 (MQ=255)
                                                                                                                                aTCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGc                              >  4:253796/1‑136 (MQ=255)
                                                                                                                                aTCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGc                              <  3:253796/136‑1 (MQ=255)
                                                                                                                                   gAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAAc                                                      >  1:135865/1‑109 (MQ=255)
                                                                                                                                   gAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAAc                                                      <  2:135865/109‑1 (MQ=255)
                                                                                                                                   gAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCAccc              <  1:362116/149‑1 (MQ=255)
                                                                                                                                   gAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCAccc              <  4:14934/149‑1 (MQ=255)
                                                                                                                                    aaCCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTGGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCAcccc             <  4:109331/149‑1 (MQ=255)
                                                                                                                                      ccACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCg           >  4:655705/1‑149 (MQ=255)
                                                                                                                                       cacaCCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGt          >  2:472925/1‑149 (MQ=255)
                                                                                                                                         cacCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                                                                                                      >  1:16548/1‑39 (MQ=255)
                                                                                                                                         cacCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTc                                                                                                                      <  2:16548/39‑1 (MQ=255)
                                                                                                                                         cacCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCt        >  2:675410/1‑149 (MQ=255)
                                                                                                                                         cacCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCt        <  1:535897/149‑1 (MQ=255)
                                                                                                                                          acCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTg       >  2:434208/1‑149 (MQ=255)
                                                                                                                                            cAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCAt                        >  3:624722/1‑130 (MQ=255)
                                                                                                                                            cAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCAt                        <  4:624722/130‑1 (MQ=255)
                                                                                                                                               cccGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAGTGTCGTGGGTGCCGACGCTCACCTCCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCAAGTTGCAGGCATCCGGTCACCCCCGTCTGCCAgc  <  3:311300/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ATCCGGCGCGACCACGCCGCTCTCTGTAGAGGTCGATTTTAACTGCCAGTGAGGAAAACCGAACATTGCCATCAGCTGTCTCCTTAAGGAGGTTAACAAGCAGGGCGCATCAGCGCGTGATAACTGCGATCGAACCACACCAGCCCGTAGGGTGTGGTGTGACGATGAATCGCTTCGATGGCGCAAAACAGAATGTCGTGGGTGCCGACGCTCACCACCTGGCTGATACGGCAGTCAAACGAAACCAGAGCCTCTTCCAGTTGCGGGCATCCGGTCACCCCCGTCTGCCAGC  >  NC_000913/1069768‑1070059

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑