Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,411,675 A→G 5.5% Y92H (TAT→CAT)  ttcA ← tRNA s(2)C32 thioltransferase, iron sulfur cluster protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,411,6750AG 5.5% 2679.9 / 88.3 805Y92H (TAT→CAT) ttcAtRNA s(2)C32 thioltransferase, iron sulfur cluster protein
Reads supporting (aligned to +/- strand):  ref base A (405/354);  new base G (24/21);  total (430/375)
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.35e-01

Only 100 of 814 total aligned reads displayed.

CCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATG  >  NC_000913/1411527‑1411823
                                                                                                                                                    |                                                                                                                                                    
     ttCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCgg                                                                                                                                                 >  1:1018449/1‑149 (MQ=255)
        cGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCa                                                                                                                                              <  4:5156043/149‑1 (MQ=255)
         gTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAg                                                                                                                                             >  3:1894969/1‑149 (MQ=255)
         gTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCACCTTTTCAAGATACTCGGCCAg                                                                                                                                             >  4:4155363/1‑149 (MQ=255)
           tGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGaa                                                                                                                                           <  2:4241187/149‑1 (MQ=255)
              ggTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATGCTCGGGCAGAACgt                                                                                                                                        >  4:597584/1‑149 (MQ=255)
                   gATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATGCTCGGGCAGAACGTGTTcc                                                                                                                                   <  1:5841410/149‑1 (MQ=255)
                     tAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCgg                                                                                                                                 <  4:3336976/149‑1 (MQ=255)
                        aGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGaa                                                                                                                              >  1:5146389/1‑149 (MQ=255)
                               cgcgACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGt                                                                                                                       <  4:5285874/149‑1 (MQ=255)
                               cgcgACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGt                                                                                                                       >  3:1699520/1‑149 (MQ=255)
                                 cgacgaAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTg                                                                                                                     >  3:2938665/1‑149 (MQ=255)
                                    cgaAGGCGAGAACACGGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACATTCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCtt                                                                                                                  <  1:1424398/149‑1 (MQ=255)
                                        ggCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGa                                                                                                              >  2:3656858/1‑149 (MQ=255)
                                         gCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATGCTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTCTGAt                                                                                                             <  4:2634652/149‑1 (MQ=255)
                                         gCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGAt                                                                                                             <  3:6609218/149‑1 (MQ=255)
                                           gagaACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGaa                                                                                                                              >  3:2344690/1‑130 (MQ=255)
                                           gagaACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCg                                                                                                           >  3:6829146/1‑149 (MQ=255)
                                            agaACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTGCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGa                                                                                                          <  3:6089392/149‑1 (MQ=255)
                                             gaACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAg                                                                                                         <  1:153051/149‑1 (MQ=255)
                                               acacAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGt                                                                                                       <  2:6092167/149‑1 (MQ=255)
                                                 acaGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTa                                                                                                     >  1:450081/1‑149 (MQ=255)
                                                    gTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGaa                                                                                                                              >  3:3103988/1‑121 (MQ=255)
                                                    gTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACa                                                                                                  <  4:6590381/149‑1 (MQ=255)
                                                     tGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAg                                                                                                         >  2:2444123/1‑141 (MQ=255)
                                                     tGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAg                                                                                                         <  1:2444123/141‑1 (MQ=255)
                                                        gCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAGCGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGcca                                                                                              <  4:2126510/149‑1 (MQ=255)
                                                        gCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGcca                                                                                              <  4:6378699/149‑1 (MQ=255)
                                                          aaGTGGTTTTGCCCGCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGccacc                                                                                            >  2:1838222/1‑149 (MQ=255)
                                                                ttttGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGaa                                                                                      <  3:815320/149‑1 (MQ=255)
                                                                ttttGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGaa                                                                                      >  1:6994844/1‑149 (MQ=255)
                                                                ttttGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGaa                                                                                      >  3:2068916/1‑149 (MQ=255)
                                                                     cccTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAatt                                                                                 <  3:3758018/149‑1 (MQ=255)
                                                                      ccTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACgtg                                                                                                                                       >  2:5780419/1‑94 (MQ=255)
                                                                       ctctGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattga                                                                               >  1:3733993/1‑149 (MQ=255)
                                                                        tctGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattgat                                                                              >  4:3961143/1‑149 (MQ=255)
                                                                          tGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATGCTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattgattg                                                                            <  3:5923620/149‑1 (MQ=255)
                                                                          tGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattgattg                                                                            >  1:4304745/1‑149 (MQ=255)
                                                                           ggAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTgg                                                                           >  3:4057881/1‑149 (MQ=255)
                                                                            gAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGc                                                                                         <  1:6677024/134‑1 (MQ=255)
                                                                            gAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGgc                                                                          <  1:4932804/149‑1 (MQ=255)
                                                                                 ttctCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTcc                                                                                                                                   <  4:1928140/87‑1 (MQ=255)
                                                                                  tctcTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCtt                                                                    >  4:2904239/1‑149 (MQ=255)
                                                                                    tctTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTtg                                                                  >  4:2361282/1‑149 (MQ=255)
                                                                                      tttCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTtgct                                                                >  2:6721898/1‑149 (MQ=255)
                                                                                      tttCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTtgct                                                                <  1:3424955/149‑1 (MQ=255)
                                                                                       ttCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTtgctg                                                               >  1:622286/1‑149 (MQ=255)
                                                                                           cGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCaaa                                                           <  3:6331443/149‑1 (MQ=255)
                                                                                               aCCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTgc                                                       <  4:6887867/149‑1 (MQ=255)
                                                                                                ccGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTgcg                                                      <  3:4761771/149‑1 (MQ=255)
                                                                                                 cGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTcc                                                                                                                                   <  4:4765184/71‑1 (MQ=255)
                                                                                                 cGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCTAAAAATTGATTGGCGCGCTTTGCTGCAAATTgcgc                                                     <  4:6143078/149‑1 (MQ=255)
                                                                                                 cGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTgcgc                                                     <  4:1448162/149‑1 (MQ=255)
                                                                                                 cGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTgcgc                                                     >  2:1102178/1‑149 (MQ=255)
                                                                                                  gTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCa                                                    >  2:267405/1‑149 (MQ=255)
                                                                                                   tAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAg                                                   >  2:739331/1‑149 (MQ=255)
                                                                                                     aGTATTCTCTTCAACAATCTTTTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTCTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAATTTGCGCAGaa                                                 <  3:4392973/149‑1 (MQ=255)
                                                                                                     aGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGaa                                                 <  4:5748050/149‑1 (MQ=255)
                                                                                                     aGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGaa                                                 <  1:831278/149‑1 (MQ=255)
                                                                                                           ctctTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCa                                           <  1:6435293/149‑1 (MQ=255)
                                                                                                           ctctTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCa                                           <  2:1565872/149‑1 (MQ=255)
                                                                                                            tctTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAg                                          >  2:3762249/1‑149 (MQ=255)
                                                                                                              ttCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGaa                                                                                      <  1:5691482/103‑1 (MQ=255)
                                                                                                                caacaaTCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATg                                      <  3:3259609/149‑1 (MQ=255)
                                                                                                                 aacaaTCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATgg                                     >  1:1924685/1‑149 (MQ=255)
                                                                                                                  acaaTCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGt                                    >  2:5286426/1‑149 (MQ=255)
                                                                                                                   caaTCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGta                                   <  3:2776824/149‑1 (MQ=255)
                                                                                                                    aaTCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATg                                      <  4:1115542/145‑1 (MQ=255)
                                                                                                                        ttGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCt                              >  1:7217942/1‑149 (MQ=255)
                                                                                                                        ttGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCt                              <  3:5439942/149‑1 (MQ=255)
                                                                                                                        ttGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCt                              >  1:2828648/1‑149 (MQ=255)
                                                                                                                           tACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattgat                                                                              <  3:1166337/98‑1 (MQ=255)
                                                                                                                            aCTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCa                                           >  2:6362375/1‑132 (MQ=255)
                                                                                                                             cTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCtt                         <  4:3703415/149‑1 (MQ=255)
                                                                                                                             cTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCtt                         >  4:826748/1‑149 (MQ=255)
                                                                                                                              tCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCttt                        >  2:2185045/1‑149 (MQ=255)
                                                                                                                              tCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCttt                        >  2:5578830/1‑149 (MQ=255)
                                                                                                                              tCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCttt                        >  2:6650402/1‑149 (MQ=255)
                                                                                                                              tCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCttt                        >  1:4032830/1‑149 (MQ=255)
                                                                                                                              tCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCttt                        >  3:4231571/1‑149 (MQ=255)
                                                                                                                              tCAACACCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattgat                                                                              <  1:4685604/95‑1 (MQ=255)
                                                                                                                              tCAACACCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAattgat                                                                              >  2:4685604/1‑95 (MQ=255)
                                                                                                                               cAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAg                                                                                                         >  3:3541933/1‑67 (MQ=255)
                                                                                                                                aaCGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTAc                      <  4:3070586/149‑1 (MQ=255)
                                                                                                                                aaCGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTAc                      <  1:2155433/149‑1 (MQ=255)
                                                                                                                                     ccAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCg                 >  4:4559977/1‑149 (MQ=255)
                                                                                                                                       aGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGga               <  2:3063768/149‑1 (MQ=255)
                                                                                                                                        gCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGgag              >  3:3919956/1‑149 (MQ=255)
                                                                                                                                        gCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTATCTTTACCCCCGgag              >  2:4461933/1‑149 (MQ=255)
                                                                                                                                         cTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGgaga             <  1:3157802/149‑1 (MQ=255)
                                                                                                                                          ttttCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGagag            <  4:5544039/149‑1 (MQ=255)
                                                                                                                                          ttttCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCGTTACCCCCGGagag            <  1:5454302/149‑1 (MQ=255)
                                                                                                                                               aaGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGagag            <  2:962336/144‑1 (MQ=255)
                                                                                                                                               aaGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCaaa       <  2:6966528/149‑1 (MQ=255)
                                                                                                                                               aaGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCaaa       >  2:1282026/1‑149 (MQ=255)
                                                                                                                                                aGATGCTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGATAGGCAAAc      >  4:4389926/1‑149 (MQ=255)
                                                                                                                                                 gATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAAcc     <  3:5308180/149‑1 (MQ=255)
                                                                                                                                                   tACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCat   <  2:487543/149‑1 (MQ=255)
                                                                                                                                                    aCTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCatg  >  3:5113877/1‑149 (MQ=255)
                                                                                                                                                    aCTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCatg  >  3:5686957/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCAGTTCCGTTGCGGTACGATAAAGGATACCGCGACGAAGGCGAGAACACAGTGAGCAAGTGGTTTTGCCCTCTGGAATCTTCTCTTTCACGATACCGTAAGTATTCTCTTCAACAATCTTGTACTCAACGCCCAGCTTTTCAAGATACTCGGGCAGAACGTGTTCCGGGAAGCCCGGTTGCTTTTGATCGAGGTTAACAGCCACCAGCGAAAAATTGATTGGCGCGCTTTGCTGCAAATTGCGCAGAATCTCCAGCATGGTATAGCTGTCTTTACCCCCGGAGAGGCAAACCATG  >  NC_000913/1411527‑1411823

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: