Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,362 C→A 5.4% intergenic (‑33/+38) tufA ← / ← fusA translation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,471,3620CA 5.4% 256.6 / 2.8 74intergenic (‑33/+38)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (38/31);  new base A (2/2);  total (40/34)
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.94e-01

GTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGC  >  NC_000913/3471214‑3471508
                                                                                                                                                    |                                                                                                                                                  
gTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                                                                    <  1:355353/149‑1 (MQ=255)
gTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                                                                    <  4:110278/149‑1 (MQ=255)
  tttAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                                                                                  <  1:459837/149‑1 (MQ=255)
         aGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGggagag                                                                                                                                           <  4:478675/149‑1 (MQ=255)
               gtggtgATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                                                                     >  1:33104/1‑149 (MQ=255)
                     aTTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                                                                               <  2:771793/149‑1 (MQ=255)
                      ttGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                                                                        <  3:474847/139‑1 (MQ=21)
                      ttGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                                                                              >  3:95530/1‑149 (MQ=255)
                      ttGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATGTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                                                                        >  4:474847/1‑139 (MQ=21)
                       tGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                                                                             <  1:591237/149‑1 (MQ=255)
                            cGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTaa                                                                                                                        <  2:471829/149‑1 (MQ=255)
                              gTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                      >  4:91488/1‑149 (MQ=255)
                               tCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                                                                     >  4:569508/1‑149 (MQ=255)
                               tCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAccc                                                                                                                     >  1:29844/1‑149 (MQ=255)
                                   agtagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGa                                                                                                                                    <  4:45153/130‑1 (MQ=255)
                                   agtagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGa                                                                                                                                    >  3:45153/1‑130 (MQ=255)
                                   agtagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                                                                             <  2:736479/137‑1 (MQ=255)
                                   agtagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                                                                             >  1:736479/1‑137 (MQ=255)
                                   agtagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAg                                                                                                                 <  2:21506/149‑1 (MQ=255)
                                   agtagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAg                                                                                                                 <  4:114885/149‑1 (MQ=255)
                                     tagtTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGc                                                                                                               <  1:252693/149‑1 (MQ=255)
                                        ttttACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTa                                                                                                            <  2:29844/149‑1 (MQ=255)
                                        ttttACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTa                                                                                                            >  2:556653/1‑149 (MQ=255)
                                             ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                                                                       >  4:352743/1‑149 (MQ=255)
                                             ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAc                                                                                                       >  1:406905/1‑149 (MQ=255)
                                                   tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACggg                                                                                                 >  4:258166/1‑149 (MQ=255)
                                                    cAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTTCTATAGTGCTCTCCCCTTTAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGc                                                                                                <  2:794742/149‑1 (MQ=255)
                                                     aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                      <  2:668201/126‑1 (MQ=255)
                                                     aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                      <  2:447604/126‑1 (MQ=255)
                                                     aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                      >  1:668201/1‑126 (MQ=255)
                                                     aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAcc                                                                                                                      >  1:447604/1‑126 (MQ=255)
                                                         tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCaa                                                                                           <  1:556653/149‑1 (MQ=255)
                                                                  gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCt                                                                                  <  2:406905/149‑1 (MQ=255)
                                                                      cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTc                                                                                             <  4:16137/134‑1 (MQ=255)
                                                                      cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTc                                                                                             >  3:16137/1‑134 (MQ=255)
                                                                      cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCACGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGc                                                                              >  3:787225/1‑149 (MQ=255)
                                                                        aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCATGAGCaa                                                                            <  2:442897/149‑1 (MQ=255)
                                                                                gTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTc                                                                    >  3:605855/1‑149 (MQ=255)
                                                                                   cGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc                                                                 <  4:95530/149‑1 (MQ=255)
                                                                                    ggTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGggagag                                                                                                                                           >  4:738297/1‑74 (MQ=37)
                                                                                          gTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat                                                          <  4:787225/149‑1 (MQ=255)
                                                                                           tACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                                                                             >  4:775669/1‑81 (MQ=255)
                                                                                           tACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGt                                                                                                                             <  3:775669/81‑1 (MQ=255)
                                                                                           tACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                                                                            >  1:124757/1‑82 (MQ=255)
                                                                                           tACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                                                                            <  2:124757/82‑1 (MQ=255)
                                                                                              gTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa                                                      >  1:481277/1‑149 (MQ=255)
                                                                                                   aaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCa                                                 >  3:286744/1‑149 (MQ=255)
                                                                                                        ttCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAAt                                            >  1:462838/1‑149 (MQ=255)
                                                                                                           tttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATAc                                                     <  2:512874/137‑1 (MQ=255)
                                                                                                           tttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATAc                                                     >  1:512874/1‑137 (MQ=255)
                                                                                                           tttAGACACGGCTATATTCCGTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTcc                                         >  4:793418/1‑149 (MQ=255)
                                                                                                              aGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATa                                      <  3:91488/149‑1 (MQ=255)
                                                                                                              aGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATa                                      <  3:352743/149‑1 (MQ=255)
                                                                                                                  acGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAgtgt                                  >  2:494707/1‑149 (MQ=255)
                                                                                                                    ggCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTat                                >  4:766124/1‑149 (MQ=255)
                                                                                                                     gCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGa                                                                                <  2:51224/100‑1 (MQ=255)
                                                                                                                     gCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGa                                                                                >  1:51224/1‑100 (MQ=255)
                                                                                                                      cTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTatga                              <  3:569508/149‑1 (MQ=255)
                                                                                                                           ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCAc                         >  4:605261/1‑149 (MQ=255)
                                                                                                                               ttACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa                                                      >  3:75192/1‑116 (MQ=255)
                                                                                                                               ttACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa                                                      <  4:75192/116‑1 (MQ=255)
                                                                                                                                   tataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 <  4:286744/149‑1 (MQ=255)
                                                                                                                                   tataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 >  1:5160/1‑149 (MQ=255)
                                                                                                                                     taGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGAGGCACGACCTTTGGt               >  4:459124/1‑149 (MQ=255)
                                                                                                                                       gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGAc                      <  3:688447/140‑1 (MQ=255)
                                                                                                                                       gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGAc                      >  4:688447/1‑140 (MQ=255)
                                                                                                                                           tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAgaga         >  2:295016/1‑149 (MQ=255)
                                                                                                                                           tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAgaga         >  4:398905/1‑149 (MQ=255)
                                                                                                                                             tcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAAc       >  2:551390/1‑149 (MQ=255)
                                                                                                                                               cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 >  1:373340/1‑137 (MQ=255)
                                                                                                                                               cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTg                 <  2:373340/137‑1 (MQ=255)
                                                                                                                                               cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGc     <  1:551390/149‑1 (MQ=255)
                                                                                                                                                ccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCa    <  3:493373/149‑1 (MQ=255)
                                                                                                                                                ccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCa    >  1:772966/1‑149 (MQ=255)
                                                                                                                                                  ttCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGc  >  4:651214/1‑149 (MQ=255)
                                                                                                                                                  ttCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGc  >  4:269731/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTTTAGCCAGTACGGTGGTGATTGCAGCGGTCAGAGTAGTTTTACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATACTTCAGGAATTCCATAGTGTATGATGCACGACCTTTGGTCAGAGAACGCAGC  >  NC_000913/3471214‑3471508

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: