Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,467,969 G→T 6.6% T115T (ACG→ACT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,467,9690GT 6.6% 216.1 / 2.7 69T115T (ACG→ACT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGCGCGTAATTAAC  >  NC_000913/4467827‑4468116
                                                                                                                                              |                                                                                                                                                   
aaGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGat                                                                                                                                               <  8:157182/149‑1 (MQ=255)
aaGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGat                                                                                                                                               <  3:95411/149‑1 (MQ=255)
   tGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGattta                                                                                                                                            <  8:41914/149‑1 (MQ=255)
   tGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGattta                                                                                                                                            <  6:203970/149‑1 (MQ=255)
    ggAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGatttat                                                                                                                                           <  5:149991/149‑1 (MQ=255)
     gAATCTGCCCGCGAACAACATGGTTAAAATAAATTACCCGCACAACAACCGTCGCCGTTGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGatttat                                                                                                                                           <  5:70675/148‑1 (MQ=255)
     gAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGatttatt                                                                                                                                          >  3:36144/1‑149 (MQ=255)
     gAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGatttat                                                                                                                                           <  5:146460/148‑1 (MQ=255)
          tGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTtgctg                                                                                                                                     >  7:33438/1‑149 (MQ=255)
             ccGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCgg                                                                                                                                  >  2:186274/1‑149 (MQ=255)
                cgAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCg                                                                                                                               >  3:122675/1‑149 (MQ=255)
                   acaacaTGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTa                                                                                                                            >  1:88733/1‑149 (MQ=255)
                           gTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTAAGCCACGGAAGATTTATTTGATGCGGGCGTTATCGCGCt                                                                                                                     >  1:174247/1‑148 (MQ=255)
                               aaataaatTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAAt                                                                                                                  <  2:176146/147‑1 (MQ=255)
                               aaataaatTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAAt                                                                                                                  >  1:176146/1‑147 (MQ=255)
                                aataaatTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGt                                                                                                               >  8:205188/1‑149 (MQ=255)
                                 ataaatTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTc                                                                                                              >  4:97741/1‑149 (MQ=255)
                                     atTACCAGCAAAACAACCGTAGCCGTTGTCGGTACATTTATGGGTCAGCTAGCGCAACCCCATTAATATTTAACTCACCATTATCGGCGCTATTTCTTACGCCACTGACGATTTATTTGATGCGGGCGTTATCGCGCTAATGGTCGCTa                                                                                                          <  5:177733/149‑1 (MQ=255)
                                          ccGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTc                                                                                                      >  3:77126/1‑148 (MQ=255)
                                            gCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTAc                                                                                                   <  7:237751/149‑1 (MQ=255)
                                               caacaaCCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGtt                                                                                                <  8:186009/149‑1 (MQ=255)
                                                aacaacCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTg                                                                                               >  7:39940/1‑149 (MQ=255)
                                                   aacCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTg                                                                                            <  8:176462/149‑1 (MQ=255)
                                                       gTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACt                                                                                        <  4:122675/149‑1 (MQ=255)
                                                          gCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTa                                                                                     <  1:185688/149‑1 (MQ=255)
                                                          gCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTa                                                                                     <  4:77126/149‑1 (MQ=255)
                                                              tgttgtgTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGACACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCa                                                                                 <  1:65085/143‑1 (MQ=255)
                                                                  gggTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAgg                                                                               >  5:156717/1‑147 (MQ=255)
                                                                  cggTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTCACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAgg                                                                               >  5:156718/2‑147 (MQ=255)
                                                                         tttATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGtt                                                                      <  6:156717/149‑1 (MQ=255)
                                                                         tttATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTCACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGtt                                                                      <  6:156718/149‑1 (MQ=255)
                                                                           tATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTc                                                                     <  1:36131/148‑1 (MQ=255)
                                                                            aTGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCa                                                                   >  6:25345/1‑149 (MQ=255)
                                                                             tGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCa                                                                   <  4:151791/148‑1 (MQ=255)
                                                                                 tCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACTGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTa                                                                >  3:21950/1‑147 (MQ=255)
                                                                                 tCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACTGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTa                                                                >  7:239277/1‑147 (MQ=255)
                                                                                 tCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTa                                                                >  3:81013/1‑147 (MQ=255)
                                                                                    gCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCg                                                           >  7:231966/1‑149 (MQ=255)
                                                                                     cTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCgg                                                          >  7:171007/1‑149 (MQ=255)
                                                                                      tATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCgg                                                          >  1:242576/1‑148 (MQ=255)
                                                                                           cAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATg                                                    >  4:36847/1‑149 (MQ=255)
                                                                                           cAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATg                                                    <  7:131931/149‑1 (MQ=255)
                                                                                              ccccTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACTGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCTCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGccc                                                 <  8:239277/149‑1 (MQ=255)
                                                                                              ccccTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACTGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGccc                                                 <  4:21950/149‑1 (MQ=255)
                                                                                               cccTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCt                                                >  1:134320/1‑149 (MQ=255)
                                                                                                     aaTATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAg                                           <  2:242576/148‑1 (MQ=255)
                                                                                                     aaTATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGc                                          >  5:205322/1‑149 (MQ=255)
                                                                                                      atatTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGc                                          >  7:13953/1‑148 (MQ=255)
                                                                                                       tatTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGa                                        <  8:33438/149‑1 (MQ=255)
                                                                                                         ttttACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGATCCACTAAAGCGGCAGATGCCCTGAAAGCGATg                                      <  8:39940/149‑1 (MQ=255)
                                                                                                          tttACTCACCATTCTAGGCTCTATTTATTACGCAACGGAAGATTTATTTGCTGCGGGAGTTATCGCGCTAATGGTCGCTATTTCTAGGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATg                                      <  5:242729/148‑1 (MQ=255)
                                                                                                               tCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAg                                 >  8:52720/1‑148 (MQ=255)
                                                                                                                 aCCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAATCGGCAGATGCCCTGAAAGCGATGGTCAGCaa                              >  3:62805/1‑149 (MQ=255)
                                                                                                                 aCCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGc                                <  1:186274/147‑1 (MQ=255)
                                                                                                                    aTTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATAc                           <  8:231966/149‑1 (MQ=255)
                                                                                                                     ttCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACt                          >  8:187129/1‑149 (MQ=255)
                                                                                                                        tcGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCg                       >  1:186141/1‑149 (MQ=255)
                                                                                                                          ggCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGTCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGa                      <  5:25345/148‑1 (MQ=255)
                                                                                                                          ggCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGAc                     <  8:171007/149‑1 (MQ=255)
                                                                                                                          gccgcTATTTCTTACGCCAGCGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCTATGGTCAGCAATACTGCGa                      <  2:60702/146‑1 (MQ=255)
                                                                                                                              cTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGtg                 >  2:162213/1‑149 (MQ=255)
                                                                                                                               tATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGtgc                >  5:211621/1‑149 (MQ=255)
                                                                                                                                     ttACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGcgcg          >  1:86139/1‑149 (MQ=255)
                                                                                                                                       aCGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTgcgc           <  2:134320/146‑1 (MQ=255)
                                                                                                                                        cGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATAGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGCGCGTaa       >  2:131809/1‑149 (MQ=255)
                                                                                                                                          ccACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGCGCGTAAtt     <  2:86139/149‑1 (MQ=255)
                                                                                                                                           cACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGTCAGATGCCCTGAAAGCGATGGTCAGCAATACTGAGACTGTTCTGCGCGTAAtt     >  5:5029/1‑148 (MQ=255)
                                                                                                                                             cGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGCGCGTAATTAAc  <  4:81013/149‑1 (MQ=255)
                                                                                                                                             cGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGCGCGTAATTAAc  <  6:211621/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
AAGTGGAATCTGCCCGCGAACAACATGGTGAAAATAAATTACCCGCACAACAACCGTCGCCGTGGTGGGTACATTTATGGGTCTGCTATCGCAACCCCTTTAATATTTTACTCACCATTCTCGGCGCTATTTCTTACGCCACGGAAGATTTATTTGCTGCGGGCGTTATCGCGCTAATGGTCGCTATTTCTACGTTGCTGAACTTTATTCAGGAAGCACGTTCCACTAAAGCGGCAGATGCCCTGAAAGCGATGGTCAGCAATACTGCGACGGTGCTGCGCGTAATTAAC  >  NC_000913/4467827‑4468116

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