Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,801,056 A→T 10.7% N31I (AAT→ATT)  nrdI → NrdEF cluster assembly flavodoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,801,0560AT10.7% 258.8 / 20.2 84N31I (AAT→ATT) nrdINrdEF cluster assembly flavodoxin
Reads supporting (aligned to +/- strand):  ref base A (43/32);  new base T (6/3);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 7.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

GTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTT  >  NC_000913/2800910‑2801198
                                                                                                                                                  |                                                                                                                                              
gtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATg                                                                                                                                              >  3:8366/1‑149 (MQ=255)
 tttCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGa                                                                                                                                             >  1:522453/1‑149 (MQ=255)
   tccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGc                                                                                                                                           <  1:738842/149‑1 (MQ=255)
   tccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGc                                                                                                                                           <  4:22844/149‑1 (MQ=255)
       tccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGa                                                                                                                                       <  1:431962/149‑1 (MQ=255)
        ccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                      <  3:582899/149‑2 (MQ=255)
         cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                     <  1:637214/149‑1 (MQ=255)
           gACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACgg                                                                                                                                   >  3:219250/1‑149 (MQ=255)
            aCATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGa                                                                                                                                  <  1:112226/149‑1 (MQ=255)
            aCATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGa                                                                                                                                  <  3:334266/149‑1 (MQ=255)
                gATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                              >  2:618573/1‑149 (MQ=255)
                gATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                              >  3:434029/1‑149 (MQ=255)
                gATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                              >  1:452514/1‑149 (MQ=255)
                 aTTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAg                                                                                                                             >  1:182601/1‑149 (MQ=255)
                 aTTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAg                                                                                                                             >  4:93143/1‑149 (MQ=255)
                        gtttGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGGCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACg                                                                                                                      >  4:508491/4‑149 (MQ=255)
                         tCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGa                                                                                                                     <  4:273740/149‑1 (MQ=255)
                         tCTGCATCCAGCGCCACAAGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCACGCGGTGCGGATCCCGCTCAATGAGCGGGAACGGATTCAGGGAGACGa                                                                                                                     >  4:249011/1‑149 (MQ=255)
                           tGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGc                                                                                                                   <  3:390086/149‑1 (MQ=255)
                              aTCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCtt                                                                                                                >  2:81494/1‑149 (MQ=255)
                                  aGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACAt                                                                                                            >  1:615704/1‑149 (MQ=255)
                                   gcgcCACACCCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATc                                                                                                           <  2:576182/149‑1 (MQ=255)
                                        acacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGAt                                                                                                      <  4:561223/149‑1 (MQ=255)
                                         cacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                     <  2:535067/149‑1 (MQ=255)
                                            gcgGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTg                                                                                                  >  4:589481/1‑149 (MQ=255)
                                               gCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGccc                                                                                               >  4:78050/1‑149 (MQ=255)
                                                    tGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTa                                                                                          >  3:603223/1‑149 (MQ=255)
                                                     gCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAc                                                                                         <  3:249011/149‑1 (MQ=255)
                                                      cATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcg                                                                                        >  3:9958/1‑149 (MQ=255)
                                                               gCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGc                                                                                                                                           >  4:273798/1‑89 (MQ=255)
                                                               gCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGc                                                                                                                                           <  3:273798/89‑1 (MQ=255)
                                                                    tCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACg                                                                          >  4:8112/1‑149 (MQ=255)
                                                                    tCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACg                                                                          >  3:333501/1‑149 (MQ=255)
                                                                       aCTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCt                                                                       >  1:471596/1‑149 (MQ=255)
                                                                         ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATc                                                                                                           <  1:331981/111‑1 (MQ=255)
                                                                         ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATc                                                                                                           >  2:331981/1‑111 (MQ=255)
                                                                            tcCAGCAGCTCCGATAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGgcgc                                                                  >  3:58486/1‑149 (MQ=255)
                                                                            tcCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGgcgc                                                                  >  3:30607/1‑149 (MQ=255)
                                                                             ccAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGcgcg                                                                 >  3:376211/1‑149 (MQ=255)
                                                                               agcagcTCCGAAAACACGCCGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGCCGCGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGt                                                               <  4:315427/149‑1 (MQ=255)
                                                                                gcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                     >  4:628184/1‑110 (MQ=255)
                                                                                gcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATc                                                                                                     <  3:628184/110‑1 (MQ=255)
                                                                                    cTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCAc                                                          <  3:247995/149‑1 (MQ=255)
                                                                                        gAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACa                                                      >  3:553242/1‑149 (MQ=255)
                                                                                         aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGAc                                                       <  2:656120/147‑1 (MQ=255)
                                                                                         aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGAc                                                       >  1:656120/1‑147 (MQ=255)
                                                                                         aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAg                                                     <  2:182601/149‑1 (MQ=255)
                                                                                         aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAg                                                     <  2:472756/149‑1 (MQ=255)
                                                                                           aaCACGCAGCGTGTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGt                                                   >  1:36849/1‑149 (MQ=255)
                                                                                            acacGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTa                                                  >  4:531123/1‑149 (MQ=255)
                                                                                             cacGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTaa                                                 <  2:749660/149‑1 (MQ=255)
                                                                                                gCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATTGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                              >  2:167825/1‑149 (MQ=255)
                                                                                                gCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                              >  2:675955/1‑149 (MQ=255)
                                                                                                  aGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGc                                            >  2:68078/1‑149 (MQ=255)
                                                                                                  aGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGc                                            <  4:33598/149‑1 (MQ=255)
                                                                                                   gCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCt                                           <  1:364541/149‑1 (MQ=255)
                                                                                                    cGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtt                                          <  2:169051/149‑1 (MQ=255)
                                                                                                      ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTGCATCCTGCTCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCACtttt                                        >  2:551532/1‑149 (MQ=255)
                                                                                                       tttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCttttt                                       <  4:719132/149‑1 (MQ=255)
                                                                                                          aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                    >  4:331832/1‑149 (MQ=255)
                                                                                                             gAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                              <  1:131967/136‑1 (MQ=255)
                                                                                                             gAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                              >  2:131967/1‑136 (MQ=255)
                                                                                                             gAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacga                                 >  2:583567/1‑149 (MQ=255)
                                                                                                             gAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacga                                 >  3:563396/1‑149 (MQ=255)
                                                                                                              aaCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacgac                                >  1:155290/1‑149 (MQ=255)
                                                                                                                cGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAAcgacga                              <  4:8366/149‑1 (MQ=255)
                                                                                                                       gTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAAc                       <  4:739923/149‑1 (MQ=255)
                                                                                                                        tCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAAcc                      >  2:379460/1‑149 (MQ=255)
                                                                                                                            cccGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTg                                                              <  3:549704/105‑1 (MQ=255)
                                                                                                                            cccGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTg                                                              >  4:549704/1‑105 (MQ=255)
                                                                                                                             ccGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCg                 >  1:264278/1‑149 (MQ=255)
                                                                                                                              cgcgGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGt                <  2:507822/149‑1 (MQ=255)
                                                                                                                               gcgGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGtt               <  4:454945/149‑1 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGt                                                               <  4:418557/91‑1 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGt                                                               >  3:418557/3‑93 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg       >  2:35220/3‑149 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg       >  1:316770/3‑149 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg       >  3:63405/3‑149 (MQ=255)
                                                                                                                                       gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg       >  2:758725/3‑149 (MQ=255)
                                                                                                                                        ttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcgg      >  1:495352/2‑149 (MQ=255)
                                                                                                                                        aTCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcgg      >  3:223265/1‑149 (MQ=255)
                                                                                                                                         tCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcggc     <  3:9932/149‑1 (MQ=255)
                                                                                                                                           ccGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGt   >  3:459045/1‑149 (MQ=255)
                                                                                                                                            cGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCGTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGtt  <  4:110389/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
GTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTT  >  NC_000913/2800910‑2801198

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: