Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,801,056 A→T 7.6% 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,0560AT 7.6% 391.1 / 15.7 116N31I (AAT→ATT) nrdINrdEF cluster assembly flavodoxin
Reads supporting (aligned to +/- strand):  ref base A (55/52);  new base T (3/6);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

GGTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTG  >  NC_000913/2800909‑2801202
                                                                                                                                                   |                                                                                                                                                  
ggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAAt                                                                                                                                                   <  3:143761/149‑1 (MQ=255)
ggtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCGGCCCGCGGTGCGCATCCCGCTCAAt                                                                                                                                                   <  3:835944/149‑1 (MQ=255)
 gtTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATg                                                                                                                                                  <  3:214257/149‑1 (MQ=255)
    tccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGc                                                                                                                                               <  1:390321/149‑1 (MQ=255)
     ccgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCg                                                                                                                                              <  4:532459/149‑1 (MQ=255)
      cgtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATGGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCgg                                                                                                                                             <  3:616767/149‑1 (MQ=255)
       gtccgGACATGGTTAACCGACTGCATCGGGCCCCGTACGCGGGCAGTGCATGAGCAAGCTCGTCTACTTCTCCAGCAGCTTGGAAAACACGCAGCGTTTTATCGAGCGTTTAGGTCTGACCGCGGTGCGCATCGCGCTCATTAAGAggg                                                                                                                                            <  4:113926/149‑1 (MQ=255)
       gtccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCggg                                                                                                                                            <  1:968714/149‑1 (MQ=255)
        tccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGa                                                                                                                                           <  2:435652/149‑1 (MQ=255)
         ccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                          <  1:840059/149‑2 (MQ=255)
         ccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGat                                                                                                                                          <  1:649270/149‑2 (MQ=255)
         ccgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGGGCGGGaa                                                                                                                                          <  3:1077464/149‑1 (MQ=255)
          cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                         <  4:737896/149‑1 (MQ=255)
          cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                         <  2:720956/149‑1 (MQ=255)
          cgGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                         <  1:596800/149‑1 (MQ=255)
           ggACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACg                                                                                                                                        <  1:366994/149‑1 (MQ=255)
            gACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACgg                                                                                                                                       >  4:585791/1‑149 (MQ=255)
             aCATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGa                                                                                                                                      <  4:16172/149‑1 (MQ=255)
             aCATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGa                                                                                                                                      >  3:16172/1‑149 (MQ=255)
               aTGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGAtt                                                                                                                                    >  4:729875/1‑149 (MQ=255)
                 gATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                                  >  2:80663/1‑149 (MQ=255)
                  aTTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAg                                                                                                                                 >  2:566599/1‑149 (MQ=255)
                   ttAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAgg                                                                                                                                <  1:578872/149‑1 (MQ=255)
                      aCCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                                  >  4:898242/1‑144 (MQ=255)
                      aCCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCa                                                                                                                                  <  3:898242/144‑1 (MQ=255)
                       ccGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGa                                                                                                                            >  3:965106/1‑149 (MQ=255)
                          tCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGa                                                                                                                         >  1:181888/1‑149 (MQ=255)
                            tGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGc                                                                                                                       >  4:1016184/1‑149 (MQ=255)
                             gCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGcc                                                                                                                      >  2:52235/1‑149 (MQ=255)
                               aTCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCtt                                                                                                                    <  2:117482/149‑1 (MQ=255)
                                tCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTa                                                                                                                   >  2:758696/1‑149 (MQ=255)
                                 ccAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTAc                                                                                                                  <  3:24284/149‑1 (MQ=255)
                                 ccAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTAc                                                                                                                  >  1:900393/1‑149 (MQ=255)
                                    gcgcCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATc                                                                                                               <  2:14006/149‑1 (MQ=255)
                                       ccACACGCGGCCCGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTg                                                                                                            >  4:397052/1‑149 (MQ=255)
                                         acacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                         >  1:802616/1‑118 (MQ=255)
                                         acacGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAAc                                                                                                                                         <  2:802616/118‑1 (MQ=255)
                                            cgcgGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGt                                                                                                       <  3:634117/149‑1 (MQ=255)
                                               ggCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAGCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACTAGACTTACATCCTGCTCGTggc                                                                                                    >  4:637474/1‑147 (MQ=255)
                                               ggCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGcc                                                                                                    <  4:452628/149‑1 (MQ=255)
                                               ggCCAGTGCATGAGCCAGCTCGTCGACTTCTCCAGCCGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGCACGGATTAAGGTAGACGAGCCTTCCATCCTGATCGTGcc                                                                                                    <  4:679025/149‑1 (MQ=255)
                                                 ccAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCct                                                                                                  >  4:223295/1‑149 (MQ=255)
                                                    gTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCtt                                                                                               <  3:464222/149‑1 (MQ=255)
                                                    gTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCtt                                                                                               >  3:728850/1‑149 (MQ=255)
                                                       cATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcg                                                                                            >  1:843185/1‑149 (MQ=255)
                                                              cAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTATGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGt                                                                                                                               >  2:257039/1‑107 (MQ=255)
                                                              cAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGt                                                                                                                               <  1:257039/107‑1 (MQ=255)
                                                               aGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcggcggcggc                                                                                    >  3:558225/1‑149 (MQ=255)
                                                                    gTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTAcggcggc                                                                                       >  4:219583/1‑141 (MQ=255)
                                                                    gTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTAc                                                                               >  1:137565/1‑149 (MQ=255)
                                                                    gTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTAc                                                                               >  2:1019838/1‑149 (MQ=255)
                                                                     tCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGt                                                                                                       >  3:770067/1‑124 (MQ=255)
                                                                     tCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGt                                                                                                       <  4:770067/124‑1 (MQ=255)
                                                                       tACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGc                                                                            >  3:326712/1‑149 (MQ=255)
                                                                         cTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTg                                                                          >  3:597307/1‑149 (MQ=255)
                                                                          ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGa                                                                                                                                           >  3:182617/1‑83 (MQ=255)
                                                                          ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCATTGAGCGGGa                                                                                                                                           <  4:182617/83‑1 (MQ=255)
                                                                          ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGAtt                                                                                                                                    >  2:937914/1‑90 (MQ=255)
                                                                          ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGAtt                                                                                                                                    <  1:937914/90‑1 (MQ=255)
                                                                          ttCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTg                                                                          >  4:1076748/1‑148 (MQ=255)
                                                                           tctcCAGCAGCACCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGgc                                                                        >  3:614391/1‑149 (MQ=255)
                                                                            ctcCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGgcg                                                                       >  4:904835/1‑149 (MQ=255)
                                                                               cagcagCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCgg                                                                    >  3:194039/1‑149 (MQ=255)
                                                                                 gcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCgg                                                                    <  2:347895/147‑1 (MQ=255)
                                                                                 gcagcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCgg                                                                    >  1:347895/1‑147 (MQ=255)
                                                                                   agcTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGcc                                                                >  1:1071870/1‑149 (MQ=255)
                                                                                          aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAg                                                         <  4:674082/149‑1 (MQ=255)
                                                                                          aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAg                                                         <  1:1056158/149‑1 (MQ=255)
                                                                                          aaaaCACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTACCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAg                                                         <  3:729875/149‑1 (MQ=255)
                                                                                               acGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAAt                                                    <  4:501490/149‑1 (MQ=255)
                                                                                               acGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAAt                                                    >  1:529325/1‑149 (MQ=255)
                                                                                               acGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAAt                                                    >  3:793979/1‑149 (MQ=255)
                                                                                                 gCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                                  >  1:808618/1‑149 (MQ=255)
                                                                                                 ccAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTc                                                  >  2:934923/2‑149 (MQ=255)
                                                                                                   aGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGc                                                >  3:469715/1‑149 (MQ=255)
                                                                                                       ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                            >  4:379856/1‑149 (MQ=255)
                                                                                                       ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                            >  2:405159/1‑149 (MQ=255)
                                                                                                       ttttATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCGTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCtttt                                            >  4:680031/1‑149 (MQ=255)
                                                                                                           aTCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTaaa                                        <  1:955546/149‑1 (MQ=255)
                                                                                                            tCGAACGTTTAGGTCTGCGCGCGGGGGGAATGCCGCTCAATGTGCGGGAACGGAGGCATGTAGACGAGCCGTACGTCCTGATCGGGCCGTCTTACGTCGGCGGCGTTACGGCTGGCGCGGTGCCACGCCTGGTAATTCACTTTATacac                                       <  2:674458/149‑1 (MQ=255)
                                                                                                            tCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAac                                       <  4:864101/149‑1 (MQ=255)
                                                                                                             cGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacg                                      >  3:105671/1‑149 (MQ=255)
                                                                                                              gAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacga                                     >  4:832590/1‑149 (MQ=255)
                                                                                                              gAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacga                                     >  4:319924/1‑149 (MQ=255)
                                                                                                               aaCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAacgac                                    >  2:682230/1‑149 (MQ=255)
                                                                                                                aCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGt                                                                                 <  1:230547/103‑1 (MQ=255)
                                                                                                                aCGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGt                                                                                 >  2:230547/1‑103 (MQ=255)
                                                                                                                      aGGTCGGCCCGCGGTGCGCATCCCGCGCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGcaca                             >  4:923594/1‑149 (MQ=255)
                                                                                                                               cgcgGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGt                    <  4:105671/149‑1 (MQ=255)
                                                                                                                                gcgGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGtt                   <  3:508980/149‑1 (MQ=255)
                                                                                                                                    tGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCtt               >  1:962182/1‑149 (MQ=255)
                                                                                                                                        gttCCCGCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg           >  4:893033/3‑149 (MQ=255)
                                                                                                                                        cATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTcgcg           <  3:379856/149‑1 (MQ=255)
                                                                                                                                         aTCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCgcgg          >  4:494485/1‑149 (MQ=255)
                                                                                                                                            ccGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGt       <  3:397052/149‑1 (MQ=255)
                                                                                                                                              gCTCATTGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTa     >  4:231554/1‑149 (MQ=255)
                                                                                                                                              gCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTa     >  3:465351/1‑149 (MQ=255)
                                                                                                                                               cTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTAt    <  2:567320/149‑1 (MQ=255)
                                                                                                                                                tCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTAtt   <  1:682230/149‑1 (MQ=255)
                                                                                                                                                 cAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTg  <  1:1019838/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GGTTTCCGTCCGGACATGATTAACCGTCTGCATCCAGCGCCACACGCGGCCAGTGCATGAGCCAGCTCGTCTACTTCTCCAGCAGCTCCGAAAACACGCAGCGTTTTATCGAACGTTTAGGTCTGCCCGCGGTGCGCATCCCGCTCAATGAGCGGGAACGGATTCAGGTAGACGAGCCTTACATCCTGATCGTGCCCTCTTACGGCGGCGGCGGTACGGCTGGCGCGGTGCCACGACAGGTAATTCGCTTTTTAAACGACGAGCACAACCGGGCGTTGCTTCGCGGCGTTATTG  >  NC_000913/2800909‑2801202

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