Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,093,463 A→G 7.4% S270S (AGT→AGC rhaD ← rhamnulose‑1‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,093,4630AG 7.4% 335.5 / 13.7 108S270S (AGT→AGCrhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (55/45);  new base G (4/4);  total (59/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.99e-01

Only 100 of 108 total aligned reads displayed.

GTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCA  >  NC_000913/4093316‑4093607
                                                                                                                                                   |                                                                                                                                                
gTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTTGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGAATCATGTTTGCTCCTTACAGCGCCAGCGCAc                                                                                                                                                 >  4:22722/1‑149 (MQ=255)
      aTGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGa                                                                                                                                           <  3:299030/149‑1 (MQ=255)
       tGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGTGCGCTGGCgc                                                                                                                                           <  2:196831/148‑2 (MQ=255)
         cgcACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGt                                                                                                                                         >  1:260624/1‑148 (MQ=255)
              cGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTa                                                                                                                                   >  4:756505/1‑149 (MQ=255)
                tGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGCCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAAc                                                                                                                                 >  4:567714/1‑149 (MQ=255)
                  agaTTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACg                                                                                                                                <  1:92170/148‑1 (MQ=255)
                          gCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACg                                                                                                                        >  1:940981/1‑148 (MQ=255)
                             cTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTg                                                                                                                    <  4:145776/149‑1 (MQ=255)
                              tGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGc                                                                                                                   >  3:624848/1‑149 (MQ=255)
                              tGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGc                                                                                                                   <  3:27940/149‑1 (MQ=255)
                              tGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGc                                                                                                                   >  4:332351/1‑149 (MQ=255)
                                 ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGGACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                                >  4:624080/1‑149 (MQ=255)
                                 ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                     <  3:608729/144‑1 (MQ=255)
                                 ttCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                     >  4:608729/1‑144 (MQ=255)
                                  tCCAGTTCTGGCCAGACGCGATTATGCCGACGCTGATACTCTTCGTGGGCGTCCGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGa                                                                                                                <  1:748914/148‑1 (MQ=255)
                                      gTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCt                                                                                                           <  3:56616/149‑1 (MQ=255)
                                        tCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTa                                                                                                          >  2:385390/1‑148 (MQ=255)
                                           ggCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCa                                                                                                       <  2:759418/148‑1 (MQ=255)
                                             ccAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAAc                                                                                                     >  1:680818/1‑148 (MQ=255)
                                              cAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAActc                                                                                                   <  4:767526/149‑1 (MQ=255)
                                              cAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAActc                                                                                                   <  4:364075/149‑1 (MQ=255)
                                               aGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAActc                                                                                                   >  2:871388/1‑148 (MQ=255)
                                                gATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCtt                                                                                                 <  3:712241/149‑1 (MQ=255)
                                                     gATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACg                                                                                             <  1:784074/148‑1 (MQ=255)
                                                      aTTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACgg                                                                                            >  1:931760/1‑148 (MQ=255)
                                                      aTTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTATACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCTAGCGCTATCAACTCTTCACGGc                                                                                           >  4:786253/1‑149 (MQ=255)
                                                       ttATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACATCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCt                                                                                          >  3:302812/1‑149 (MQ=255)
                                                       ttATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTg                                                                                                                                        >  3:437477/1‑103 (MQ=255)
                                                       ttATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTg                                                                                                                                        <  4:437477/103‑1 (MQ=255)
                                                            ccGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATgg                                                                                     <  4:291496/149‑1 (MQ=255)
                                                              gACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATTGTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGt                                                                                    <  1:883898/148‑1 (MQ=255)
                                                              gACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGt                                                                                    <  2:931760/148‑1 (MQ=255)
                                                               aCGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCGTGACAAAGGCTTTGCGGAGCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGCAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTc                                                                                   >  1:301991/1‑148 (MQ=255)
                                                                  cTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                                >  2:876599/1‑148 (MQ=255)
                                                                      tACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCa                                                                            <  1:379220/148‑1 (MQ=255)
                                                                      tACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                           >  3:300635/1‑149 (MQ=255)
                                                                       aCTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCAAGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCAt                                                                           >  2:483169/1‑148 (MQ=255)
                                                                         tctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgcc                                                                        >  3:597530/1‑149 (MQ=255)
                                                                          ctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                     <  2:777068/103‑1 (MQ=255)
                                                                          ctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCtt                                                                                                                     >  1:777068/1‑103 (MQ=255)
                                                                          ctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgct                                                                        >  2:697751/1‑147 (MQ=255)
                                                                          ctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgcc                                                                        >  1:546423/1‑148 (MQ=255)
                                                                          ctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgcc                                                                        >  2:136510/1‑148 (MQ=255)
                                                                           ttCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTAGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccg                                                                       >  2:732484/1‑148 (MQ=255)
                                                                             cGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGccc                                                                    <  3:82255/149‑1 (MQ=255)
                                                                                 ggCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATc                                                                 >  1:401625/1‑148 (MQ=255)
                                                                                       ggggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                          <  3:453060/149‑1 (MQ=255)
                                                                                        gggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                          <  1:705829/148‑1 (MQ=255)
                                                                                        gggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                         >  3:269448/1‑149 (MQ=255)
                                                                                        gggTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAAc                                                         >  4:266869/1‑149 (MQ=255)
                                                                                            ttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTc                                                     <  3:656213/149‑1 (MQ=255)
                                                                                              aCCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCa                                                    <  2:253426/148‑1 (MQ=255)
                                                                                                cTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTa                                                 >  4:413687/1‑149 (MQ=255)
                                                                                                 tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTa                                                 >  2:809999/1‑148 (MQ=255)
                                                                                                   cATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAAt                                               >  2:826679/1‑148 (MQ=255)
                                                                                                    aTGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATa                                              >  2:668781/1‑148 (MQ=255)
                                                                                                      gACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACtt                                           >  4:571403/1‑149 (MQ=255)
                                                                                                      gACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACt                                            <  2:260624/148‑1 (MQ=255)
                                                                                                        cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                          >  2:646367/1‑148 (MQ=255)
                                                                                                        cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                          >  2:896649/1‑148 (MQ=255)
                                                                                                         aaaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtgtg                                        >  3:686086/1‑149 (MQ=255)
                                                                                                          aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGc                                       <  4:369840/149‑1 (MQ=255)
                                                                                                           aGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCt                                      >  4:522633/1‑149 (MQ=255)
                                                                                                                ttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAtt                                  <  1:136510/148‑1 (MQ=255)
                                                                                                                ttGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAtt                                  <  1:646367/148‑1 (MQ=255)
                                                                                                                   cGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTc                              <  3:637275/149‑1 (MQ=255)
                                                                                                                    ggATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCt                             <  3:571403/149‑1 (MQ=255)
                                                                                                                       tCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGc                           >  1:933984/1‑148 (MQ=255)
                                                                                                                        cATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaa                                                           <  1:550827/115‑1 (MQ=255)
                                                                                                                        cATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaa                                                           >  2:550827/1‑115 (MQ=255)
                                                                                                                         aTGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCgg                         <  2:678380/148‑1 (MQ=255)
                                                                                                                            tttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGATTAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                     <  3:624080/149‑1 (MQ=255)
                                                                                                                              tGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCg                    <  3:468813/148‑1 (MQ=255)
                                                                                                                              tGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCg                    >  4:468813/1‑148 (MQ=255)
                                                                                                                               gCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGAt                  <  3:332351/149‑1 (MQ=255)
                                                                                                                                 tCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTa                >  3:866802/1‑149 (MQ=255)
                                                                                                                                 tCCTTACAGCGCCAGCGCACTGGCGAGTGGAGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGCCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTa                >  4:330547/1‑149 (MQ=255)
                                                                                                                                  ccTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTa                >  2:825381/1‑148 (MQ=255)
                                                                                                                                   cTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGtt                                                                               >  4:842990/1‑84 (MQ=255)
                                                                                                                                   cTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGtt                                                                               <  3:842990/84‑1 (MQ=255)
                                                                                                                                    ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTaaa              >  1:105155/1‑148 (MQ=255)
                                                                                                                                    ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAc             >  4:282443/1‑149 (MQ=255)
                                                                                                                                     tACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTc                              <  2:220664/131‑1 (MQ=255)
                                                                                                                                     tACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTc                              >  1:220664/1‑131 (MQ=255)
                                                                                                                                      cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                     <  4:301050/135‑1 (MQ=255)
                                                                                                                                      cactcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTc                     >  3:301050/5‑139 (MQ=255)
                                                                                                                                      aCAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAAcc            <  1:192373/148‑1 (MQ=255)
                                                                                                                                        agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa         >  4:148522/1‑149 (MQ=255)
                                                                                                                                        agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGa          <  2:933984/148‑1 (MQ=255)
                                                                                                                                        agcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGa          <  1:825381/148‑1 (MQ=255)
                                                                                                                                         tcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa         >  2:195876/2‑148 (MQ=255)
                                                                                                                                         tcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGa          >  1:164428/2‑147 (MQ=255)
                                                                                                                                         gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGaa         >  2:220392/1‑148 (MQ=255)
                                                                                                                                           gccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCAc                                                                                              >  4:463041/1‑61 (MQ=255)
                                                                                                                                           gccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCAc                                                                                              <  3:463041/61‑1 (MQ=255)
                                                                                                                                           gccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGt      >  3:538551/1‑149 (MQ=255)
                                                                                                                                               gcgcACTGGCGAGTGGCGTACCGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACCGCTGCTGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACCTGTGCTGATTTTTCTGCGGTGTCGTTTAAACCCAAGGttgc   <  2:525313/148‑3 (MQ=255)
                                                                                                                                               gcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCa  <  4:302812/149‑1 (MQ=255)
                                                                                                                                                cgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCa  >  1:942647/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTAGTTATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTTCA  >  NC_000913/4093316‑4093607

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