Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,605,059 T→C 5.2% F83F (TTT→TTC hyfD → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,605,0590TC 5.2% 238.9 / 2.2 76F83F (TTT→TTChyfDhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base C (2/2);  total (39/37)
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.87e-01

CGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCT  >  NC_000913/2604912‑2605204
                                                                                                                                                   |                                                                                                                                                 
cgccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGttt                                                                                                                                                   <  7:296850/148‑1 (MQ=255)
cgccgcGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGttt                                                                                                                                                   <  1:267686/148‑1 (MQ=255)
    gcgcTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTgtggt                                                                                                                                              >  1:26412/1‑149 (MQ=255)
           ccACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGGGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGttt                                                                                                                                       >  7:151807/1‑149 (MQ=255)
               gCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCAGTTTGTGGTGGTGTTTCTCg                                                                                                                                   >  2:143416/1‑149 (MQ=255)
                     cATGTTGTCACTGATCCCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGACAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGc                                                                                                                             >  7:215563/1‑149 (MQ=255)
                       ggttggCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGGCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTg                                                                                                                           <  7:188189/148‑1 (MQ=255)
                            tCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCAc                                                                                                                      >  4:81053/1‑149 (MQ=255)
                              aCTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTgtggtgg                                                                                                                                            <  5:272523/125‑1 (MQ=255)
                              aCTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTgtggtgg                                                                                                                                            >  6:272523/1‑125 (MQ=255)
                              aCTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGa                                                                                                                    >  3:117099/1‑149 (MQ=255)
                               cTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGGGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGAt                                                                                                                   <  5:103976/149‑1 (MQ=255)
                                tGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATc                                                                                                                  >  1:286770/1‑149 (MQ=255)
                                  aTCTCCGCGTTTTATCAGGCAGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATc                                                                                                                  <  1:156181/147‑1 (MQ=255)
                                     tcCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTc                                                                                                             >  1:6182/1‑149 (MQ=255)
                                     tcCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTc                                                                                                             >  8:118629/1‑149 (MQ=255)
                                       cgcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTc                                                                                                             <  4:117099/147‑1 (MQ=255)
                                        gcgTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCAc                                                                                                          <  6:129551/149‑1 (MQ=255)
                                        gcgGTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGGCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCAc                                                                                                          <  3:17557/149‑1 (MQ=255)
                                          gTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACgg                                                                                                        <  2:53951/149‑1 (MQ=255)
                                             ttATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGt                                                                                                      <  3:232265/148‑1 (MQ=255)
                                                tCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTAt                                                                                                   <  1:193097/148‑1 (MQ=255)
                                                       gATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGa                                                                                           >  4:311342/1‑148 (MQ=255)
                                                        aTAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                                                                                              >  6:225769/1‑145 (MQ=255)
                                                        aTAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGAc                                                                                              <  5:225769/145‑1 (MQ=255)
                                                         tAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGAt                                                                                          >  6:265005/1‑148 (MQ=255)
                                                                 ccGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgc                                                                                 <  3:81053/149‑1 (MQ=255)
                                                                  cGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATcgc                                                                                 >  5:4510/1‑148 (MQ=255)
                                                                       cgttgacgttgGTCAACGTGGGGGATGTGGCGTGGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAacac                                                                           <  5:250830/149‑1 (MQ=255)
                                                                        gttgacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGa                                                                                                                    <  5:126436/107‑1 (MQ=255)
                                                                        gttgacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGa                                                                                                                    >  6:126436/1‑107 (MQ=255)
                                                                           gacgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCg                                                                        >  4:321192/1‑148 (MQ=255)
                                                                            acgttgGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCa                                                                      >  4:39951/1‑149 (MQ=255)
                                                                                   tCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTTAGTACGCTGATTCAGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCTCATAACGGc                                                               <  1:143416/149‑1 (MQ=255)
                                                                                   tCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTtct                                                                                                                                     <  4:280390/79‑1 (MQ=255)
                                                                                   tCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTtct                                                                                                                                     >  3:280390/1‑79 (MQ=255)
                                                                                   tCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                                                                                         <  3:32394/107‑1 (MQ=255)
                                                                                   tCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACg                                                                                                         >  4:32394/1‑107 (MQ=255)
                                                                                   tCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACgg                                                                >  4:255756/1‑148 (MQ=255)
                                                                                    cAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCa                                                              <  3:39951/149‑1 (MQ=255)
                                                                                        gTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGaa                                                          >  7:59475/1‑149 (MQ=255)
                                                                                          gggggATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATc                                                        >  1:141760/1‑149 (MQ=255)
                                                                                          gggggATGTGGCGGTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATc                                                        <  2:286770/149‑1 (MQ=255)
                                                                                             ggATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGtt                                                     <  3:321192/149‑1 (MQ=255)
                                                                                              gATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGtta                                                    <  3:259715/149‑1 (MQ=255)
                                                                                                  tgGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGttatta                                                 <  5:119507/148‑1 (MQ=255)
                                                                                                  tgGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTAc                                                <  3:255756/149‑1 (MQ=255)
                                                                                                   ggCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACg                                               >  6:352240/1‑149 (MQ=255)
                                                                                                   ggCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTCGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACg                                               >  6:18653/1‑149 (MQ=255)
                                                                                                   ggCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTCGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACg                                               >  6:23878/1‑149 (MQ=255)
                                                                                                         gttTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCAtttt                                         >  6:328088/1‑149 (MQ=255)
                                                                                                             ggCCTGGTCATTGATCGCGGGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGGTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTAc                                      <  2:367454/148‑1 (MQ=255)
                                                                                                              gCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTg                                    <  2:26412/149‑1 (MQ=255)
                                                                                                              gCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTg                                    >  8:107333/1‑149 (MQ=255)
                                                                                                                cTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGgt                                  >  8:270393/1‑149 (MQ=255)
                                                                                                                  ggTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa                                                                               <  4:232655/102‑1 (MQ=255)
                                                                                                                  ggTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGa                                                                               >  3:232655/1‑102 (MQ=255)
                                                                                                                   gTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAGTCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGtt                               >  8:109962/1‑149 (MQ=255)
                                                                                                                          aTCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGtt                                                     <  2:389636/120‑1 (MQ=255)
                                                                                                                          aTCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGtt                                                     >  1:389636/1‑120 (MQ=255)
                                                                                                                          aTCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGgc                        <  7:270393/149‑1 (MQ=255)
                                                                                                                             gcgTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGgcgc                      >  2:232484/1‑148 (MQ=255)
                                                                                                                             gcgTGAGTACGCTGATTCTGTTCGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGgcgc                      <  5:23878/148‑1 (MQ=255)
                                                                                                                             gcgTGAGTACGCTGATTCTGTTCGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGGGTTTATCGgcgc                      <  5:18653/148‑1 (MQ=255)
                                                                                                                                tGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATg                  >  5:18695/1‑149 (MQ=255)
                                                                                                                                 gAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATgg                 >  2:15649/1‑149 (MQ=255)
                                                                                                                                   gTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGAACGAATCGTTTTTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCg               >  1:351020/1‑149 (MQ=255)
                                                                                                                                    tACGCTGATTCTGTTTGTGGTGGGGTTTCTCGGTTTGCTGGGCACGATCTACTGCACGGGTTATCTGTCGGATAAAAATCGCGAACAGCCGCATAACGGCACGAATGGTTATTACGCATTTTTTCTGGTGTTTATCGGCGCGATGGCgg              <  8:215563/149‑1 (MQ=255)
                                                                                                                                      cGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGa            <  1:232484/149‑1 (MQ=255)
                                                                                                                                         tGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATc                                                                                   <  2:97682/75‑1 (MQ=255)
                                                                                                                                         tGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATc                                                                                   >  1:97682/1‑75 (MQ=255)
                                                                                                                                         tGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTg         <  6:18695/149‑1 (MQ=255)
                                                                                                                                         tGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTg         >  7:107935/1‑149 (MQ=255)
                                                                                                                                           aTTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTgg        >  6:122800/1‑148 (MQ=255)
                                                                                                                                           aTTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGt       <  2:141760/149‑1 (MQ=255)
                                                                                                                                           aTTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGGCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGGGTTTATCGGCGCGATGGCGGGACTg         <  6:4510/147‑1 (MQ=255)
                                                                                                                                             tCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTa      <  8:59475/148‑1 (MQ=255)
                                                                                                                                                gTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTActct  >  4:183606/1‑149 (MQ=255)
                                                                                                                                                 tttGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTActct  <  7:107333/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CGCCGCGCTGACCACGCTGTGCATGTTGTCACTGATCTCCGCGTTTTATCAGGCCGATAAAGTTGCCGTCACGTTGACGTTGGTCAACGTGGGGGATGTGGCGTTGTTTGGCCTGGTCATTGATCGCGTGAGTACGCTGATTCTGTTTGTGGTGGTGTTTCTCGGTTTGCTGGTCACGATCTACTCCACGGGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCT  >  NC_000913/2604912‑2605204

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