Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,847,724 C→T 5.1% L506L (CTG→CTA hycC ← hydrogenase 3, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,847,7240CT 5.1% 259.9 / 2.1 78L506L (CTG→CTAhycChydrogenase 3, membrane subunit
Reads supporting (aligned to +/- strand):  ref base C (29/45);  new base T (2/2);  total (31/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCTCATC  >  NC_000913/2847577‑2847871
                                                                                                                                                   |                                                                                                                                                   
cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagca                                                                                                                                                    <  4:54670/149‑1 (MQ=255)
cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagc                                                                                                                                                     <  1:230323/148‑1 (MQ=255)
cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagc                                                                                                                                                     <  1:554333/148‑1 (MQ=255)
cGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagc                                                                                                                                                     <  2:338165/148‑1 (MQ=255)
 ggCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagcag                                                                                                                                                   <  4:595417/149‑1 (MQ=255)
 ggCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATcagcag                                                                                                                                                   <  4:191299/149‑1 (MQ=255)
       gCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGt                                                                                                                                             >  4:157763/1‑149 (MQ=255)
         aaaaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTg                                                                                                                                            <  1:942208/148‑1 (MQ=255)
         aaaaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGa                                                                                                                                           <  4:447730/149‑1 (MQ=255)
          aaaCCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGat                                                                                                                                          <  4:706613/149‑1 (MQ=255)
                tGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCg                                                                                                                                     <  2:555855/148‑1 (MQ=255)
                    cGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTga                                                                                                                                <  3:44700/149‑1 (MQ=255)
                       tAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaa                                                                                                                             >  4:500335/1‑149 (MQ=255)
                        aaTCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGaaa                                                                                                                            >  4:651344/1‑149 (MQ=255)
                            aCCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggt                                                                                                                        <  3:267699/149‑1 (MQ=255)
                            aCCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACgg                                                                                                                         >  1:156444/1‑148 (MQ=255)
                             ccATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtg                                                                                                                       >  3:246899/1‑149 (MQ=255)
                             ccATTGATTTTTCGTGGTCGTAACCGCACACCAAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtg                                                                                                                       <  3:533300/149‑1 (MQ=255)
                              cATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACggtg                                                                                                                       >  2:435118/1‑148 (MQ=255)
                                tggATTTTTCGTGGTCGTAAACGCACACCCAGGCCGCTCCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATCATGAAGGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGgtggtg                                                                                                                    <  4:211283/147‑1 (MQ=255)
                                 tGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAACGATGAATGGCAGCAGCGGGCAGGCAAACAGCAGCAAAGTGATCATCGGTTGAGAAACGGTGgtgt                                                                                                                   <  4:265749/149‑1 (MQ=255)
                                    tttttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAg                                                                                                                >  3:928511/1‑149 (MQ=255)
                                    tttttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAg                                                                                                                <  3:866658/149‑1 (MQ=255)
                                     ttttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAg                                                                                                                >  2:623220/1‑148 (MQ=255)
                                      tttCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGc                                                                                                               <  2:689134/148‑1 (MQ=255)
                                         cGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAgg                                                                                                            <  1:729198/148‑1 (MQ=255)
                                           tGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGTAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCt                                                                                                          >  1:820271/1‑148 (MQ=255)
                                           tGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGTAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCt                                                                                                          >  1:399545/1‑148 (MQ=255)
                                           tGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCt                                                                                                          <  2:353734/148‑1 (MQ=255)
                                             gTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCa                                                                                                       >  3:53297/1‑149 (MQ=255)
                                                 tAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCg                                                                                                    <  2:782629/148‑1 (MQ=255)
                                                  aaCCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGc                                                                                                  >  3:245265/1‑149 (MQ=255)
                                                      gCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCaga                                                                                               <  1:727885/148‑1 (MQ=255)
                                                      gCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCaga                                                                                               >  1:556717/1‑148 (MQ=255)
                                                      gCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCAATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCagag                                                                                              <  3:110469/149‑1 (MQ=255)
                                                         cacCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGt                                                                                            <  1:599187/148‑1 (MQ=255)
                                                          acCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAc                                                                                          <  4:933465/149‑1 (MQ=255)
                                                           cccAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAc                                                                                          >  1:901758/1‑148 (MQ=255)
                                                            ccAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAca                                                                                         >  2:960679/1‑148 (MQ=255)
                                                             cAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAAGGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcag                                                                                        <  2:901758/148‑1 (MQ=255)
                                                              aGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCaa                                                                                                                                                >  2:201878/1‑91 (MQ=255)
                                                              aGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCaa                                                                                                                                                <  1:201878/91‑1 (MQ=255)
                                                              aGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagca                                                                                      <  3:833789/149‑1 (MQ=255)
                                                                gCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcagc                                                                                    <  3:258376/149‑1 (MQ=255)
                                                                 ccGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcagca                                                                                   <  4:407917/149‑1 (MQ=255)
                                                                 ccGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcagca                                                                                   <  4:547495/149‑1 (MQ=255)
                                                                   gCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTAcagcagcag                                                                                  <  2:239920/148‑1 (MQ=255)
                                                                     aCCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGTAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                                                                <  2:399545/148‑1 (MQ=255)
                                                                     aCCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                                                                <  2:171634/148‑1 (MQ=255)
                                                                     aCCGCGGGAACGCGATGGCAAACGATCGCCTTTCCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGTAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCATAGGTACAGCAGCagag                                                                                <  2:820271/148‑1 (MQ=255)
                                                                       cgcgGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGc                                                                             <  3:264564/149‑1 (MQ=255)
                                                                       cgcgGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGc                                                                             <  4:719277/149‑1 (MQ=255)
                                                                         cgGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCa                                                                            >  2:266346/1‑148 (MQ=255)
                                                                          gggAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGat                                                                                                                                          <  3:326066/85‑1 (MQ=255)
                                                                          gggAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGat                                                                                                                                          >  4:326066/1‑85 (MQ=255)
                                                                                gcgATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAAc                                                                                                                                               >  4:427991/1‑74 (MQ=255)
                                                                                gcgATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAAc                                                                                                                                               <  3:427991/74‑1 (MQ=255)
                                                                                 cgATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGCGGTACAGCAGCAGAGAGCATCGGCAGt                                                                    <  2:842885/148‑1 (MQ=255)
                                                                                 cgATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGt                                                                    <  1:266346/148‑1 (MQ=255)
                                                                                  gATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAg                                                                  >  4:38083/1‑149 (MQ=255)
                                                                                        aaaCGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACgg                                                            >  3:672177/1‑149 (MQ=255)
                                                                                         aaCGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGgc                                                           <  4:245265/149‑1 (MQ=255)
                                                                                              tCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCgcag                                                       >  1:891208/1‑148 (MQ=255)
                                                                                                  cTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                                                                <  1:15869/119‑1 (MQ=255)
                                                                                                  cTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCagag                                                                                >  2:15869/1‑119 (MQ=255)
                                                                                                   tttGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAaca                                                  <  2:737829/148‑1 (MQ=255)
                                                                                                       cAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGc                                              >  1:757405/1‑148 (MQ=255)
                                                                                                         aaTCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGCCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCaa                                           <  4:246899/149‑1 (MQ=255)
                                                                                                          aTCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTCAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAAt                                          >  4:69640/1‑149 (MQ=255)
                                                                                                           tCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAAt                                          >  1:539088/1‑148 (MQ=255)
                                                                                                              ccATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTAc                                       <  2:851903/148‑1 (MQ=255)
                                                                                                               cATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACg                                      >  1:417202/1‑148 (MQ=255)
                                                                                                                aTAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACgc                                     <  2:392802/148‑1 (MQ=255)
                                                                                                                      atgaATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCaga                                                                                               >  1:205971/1‑84 (MQ=255)
                                                                                                                      atgaATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCaga                                                                                               <  2:205971/84‑1 (MQ=255)
                                                                                                                         aaTGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCGGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGc                           <  3:846132/149‑1 (MQ=255)
                                                                                                                           tGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCa                         <  3:157763/149‑1 (MQ=255)
                                                                                                                            ggCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGgtggtg                                                                                                                    <  2:960967/57‑1 (MQ=255)
                                                                                                                            ggCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGgtggtg                                                                                                                    >  1:960967/1‑57 (MQ=255)
                                                                                                                             gcagcaGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGt                       <  4:686084/149‑1 (MQ=255)
                                                                                                                                gcagcGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGc                     >  1:39846/1‑148 (MQ=255)
                                                                                                                                  agcGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTCAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCAc                  <  3:69640/149‑1 (MQ=255)
                                                                                                                                           gcaATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACtacgctta          <  1:758555/148‑1 (MQ=255)
                                                                                                                                               tCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCtc     >  3:38534/1‑149 (MQ=255)
                                                                                                                                                cagcagCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCtca    <  4:58529/149‑1 (MQ=255)
                                                                                                                                                  gcagcaACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCtcatc  <  3:500335/149‑1 (MQ=255)
                                                                                                                                                  gcagcaACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCtcat   <  1:813860/148‑1 (MQ=255)
                                                                                                                                                   cagcaACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCtcatc  >  2:688376/1‑148 (MQ=255)
                                                                                                                                                   cagcaACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCtcatc  >  2:11028/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CGGCATGGCAAAACCGTGAGCGGTAATCACCATTGATTTTTCGTGGTCGTAACCGCACACCCAGGCCGCACCGCGGGAACGCGATGGCAAACGATCGCCTTTGCAAATCGCCATAATGATGAATGGCAGCAGCGGGCAGGCAATCAGCAGCAACGTGATCATCGGTTGAGAAACGGTGGTGTTAGCAGGCTCCAGCGGCAGAGGTACAGCAGCAGAGAGCATCGGCAGTAGCCACGGCGCAGCAACACCGCCAATTACGCAGCAAATCGCCAGTGCCACTACGCTTACGCTCATC  >  NC_000913/2847577‑2847871

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