Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,331,983 G→C 6.8% I251M (ATC→ATG)  yhbE ← EamA family inner membrane putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,331,9830GC 6.8% 206.5 / 6.3 59I251M (ATC→ATG) yhbEEamA family inner membrane putative transporter
Reads supporting (aligned to +/- strand):  ref base G (31/24);  new base C (2/2);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.90e-01

GACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGA  >  NC_000913/3331835‑3332131
                                                                                                                                                    |                                                                                                                                                    
gACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatg                                                                                                                                                      >  3:85969/1‑149 (MQ=255)
 aCCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGtgatga                                                                                                                                                     >  2:23873/1‑149 (MQ=255)
          ggAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCa                                                                                                                                            <  4:285965/149‑1 (MQ=255)
            aaTACATCGCACCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGCGATCATCACGCTCAcc                                                                                                                                          <  1:50071/149‑1 (MQ=255)
                cATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcg                                                                                                                                      >  3:503390/1‑149 (MQ=255)
                 aTCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcgc                                                                                                                                     >  2:246129/1‑149 (MQ=255)
                    gcgcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgc                                                                                                                                  <  4:661644/149‑1 (MQ=255)
                     cgcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcc                                                                                                                                 <  2:224708/149‑1 (MQ=255)
                      gcCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGctgcct                                                                                                                                <  3:269685/149‑1 (MQ=255)
                                        cGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCtt                                                                                                              >  4:525504/1‑149 (MQ=255)
                                           aCCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGa                                                                                                                        >  2:327371/1‑136 (MQ=255)
                                           aCCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGa                                                                                                                        <  1:327371/136‑1 (MQ=255)
                                           aCCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCg                                                                                                           >  2:4967/1‑149 (MQ=255)
                                           aCCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCg                                                                                                           >  3:468567/1‑149 (MQ=255)
                                                    aCCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGgcgc                                                                                                  >  3:114355/1‑149 (MQ=255)
                                                       tAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGTCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCATGGCGCCat                                                                                               >  2:596302/1‑149 (MQ=255)
                                                       tAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                                                                               >  3:506443/1‑149 (MQ=255)
                                                       tAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCat                                                                                               >  4:152/1‑149 (MQ=255)
                                                         aaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCatat                                                                                             >  1:169731/1‑149 (MQ=255)
                                                         aaaGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCatat                                                                                             >  1:228233/1‑149 (MQ=255)
                                                           aGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATcc                                                                                           >  3:15892/1‑149 (MQ=255)
                                                           aGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATcc                                                                                           <  3:330915/149‑1 (MQ=255)
                                                              tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                                                                        >  1:526424/1‑149 (MQ=255)
                                                              tttAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAc                                                                                        <  4:503390/149‑1 (MQ=255)
                                                                 aaCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCaa                                                                                     <  2:115708/149‑1 (MQ=255)
                                                                  aCATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAg                                                                                    <  1:23873/149‑1 (MQ=255)
                                                                    aTCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGt                                                                                  >  2:292254/1‑149 (MQ=255)
                                                                     tCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTa                                                                                 <  3:152/149‑1 (MQ=255)
                                                                        gTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcgc                                                                                                                                     >  4:584352/1‑94 (MQ=255)
                                                                        gTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTgcgc                                                                                                                                     <  3:584352/94‑1 (MQ=255)
                                                                           tGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGt                                                                           >  4:352868/1‑149 (MQ=255)
                                                                            ggCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTc                                                                          <  1:255088/149‑1 (MQ=255)
                                                                             gCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTcc                                                                         <  3:364055/149‑1 (MQ=255)
                                                                               gaagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                       >  4:296704/1‑149 (MQ=255)
                                                                               gaagaaATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGc                                                                       >  1:376209/1‑149 (MQ=255)
                                                                                     aTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAGACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                                                                 >  2:408127/1‑149 (MQ=255)
                                                                                     aTCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCaaaaaa                                                                 <  4:85969/149‑1 (MQ=255)
                                                                                                    tGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCt                                                  >  4:152085/1‑149 (MQ=255)
                                                                                                      aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGc                                                <  4:56515/149‑1 (MQ=255)
                                                                                                      aTAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGc                                                >  1:254646/1‑149 (MQ=255)
                                                                                                        aaaaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCa                                              <  1:246129/149‑1 (MQ=255)
                                                                                                          aaGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGt                                            >  1:618776/1‑149 (MQ=255)
                                                                                                            gATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTgg                                          <  2:228233/149‑1 (MQ=255)
                                                                                                             aTCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                       >  4:516954/1‑103 (MQ=255)
                                                                                                             aTCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTAcc                                                                                       <  3:516954/103‑1 (MQ=255)
                                                                                                                tGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTaa                                      >  2:144748/1‑149 (MQ=255)
                                                                                                                  aaaaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTTAGc                                    <  1:596302/149‑1 (MQ=255)
                                                                                                                    aaaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTg                                  >  1:339449/1‑149 (MQ=255)
                                                                                                                     aaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTgc                                 >  3:542818/1‑149 (MQ=255)
                                                                                                                     aaaaCAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTgc                                 <  3:296704/149‑1 (MQ=255)
                                                                                                                            cGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCa                          <  3:352868/149‑1 (MQ=255)
                                                                                                                              tAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACt                        >  3:413167/1‑149 (MQ=255)
                                                                                                                               aaaCAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTc                       >  2:241196/1‑149 (MQ=255)
                                                                                                                                       gggTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCttt               <  1:446818/149‑1 (MQ=255)
                                                                                                                                        ggTGAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTg              >  2:5882/1‑149 (MQ=255)
                                                                                                                                         gTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGc             <  4:413167/149‑1 (MQ=255)
                                                                                                                                         gTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGc             >  3:284847/1‑149 (MQ=255)
                                                                                                                                          tGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGcc            >  1:256737/1‑149 (MQ=255)
                                                                                                                                           gAGCGCGATCATCACGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCa           >  3:147980/1‑149 (MQ=255)
                                                                                                                                               gTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAgg       <  4:506443/149‑1 (MQ=255)
                                                                                                                                               gTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAgg       <  2:254646/149‑1 (MQ=255)
                                                                                                                                                  atgatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGaa    >  3:686674/1‑149 (MQ=255)
                                                                                                                                                    gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAaga  >  4:508418/1‑149 (MQ=255)
                                                                                                                                                    gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAaga  <  1:4967/149‑1 (MQ=255)
                                                                                                                                                    gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAaga  <  4:542818/149‑1 (MQ=255)
                                                                                                                                                    gatCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAaga  <  1:241196/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACCAATGGCGGAATACATCGCGCCCGCAACCACGACAAACGCACCGAGATAACCTAAAAGGTTTAACATCGGTCTGGCGAAGAAATCGGGCCAGGCCAGTGATAAAAGATCTGAAAAAAACAGCGTAAACAGTGGGGTGAGCGTGATGATCGCGCTCACCTGCGCTGCCTGCCAGCGAGCCATCGCTTCCGCCAGGGCGCCATATCCTACCAAGGTATTCAGTCCGCAAAAAATTAAACATGCGAGCTGCCAGTGGCTAAGCTGCGCTATCACTCCAGGCTTTGCCAGAGGGAAGA  >  NC_000913/3331835‑3332131

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