Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,093,463 A→G 5.3% 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 5.3% 330.0 / 6.1 94S270S (AGT→AGCrhaDrhamnulose‑1‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (44/45);  new base G (3/2);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

TATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTT  >  NC_000913/4093321‑4093605
                                                                                                                                              |                                                                                                                                              
tATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGc                                                                                                                                           >  3:837400/1‑148 (MQ=255)
             agaTTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACg                                                                                                                              <  1:967904/148‑1 (MQ=255)
                         tGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTg                                                                                                                  >  4:930101/1‑148 (MQ=255)
                               cAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCg                                                                                                                                          <  3:413278/118‑1 (MQ=255)
                               cAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCg                                                                                                                                          >  4:413278/1‑118 (MQ=255)
                               cAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAgc                                                                                                            <  4:837400/148‑1 (MQ=255)
                                   tCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCTCAATGGCGAGTGGCGTAACGCCGAAACGCt                                                                                                                    <  4:700325/136‑1 (MQ=255)
                                   tCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTaa                                                                                                                                <  4:767102/124‑1 (MQ=255)
                                   tCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTaa                                                                                                                                >  3:767102/1‑124 (MQ=255)
                                   tCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACACCGAAACGCt                                                                                                                    >  3:700325/1‑136 (MQ=255)
                                          aGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACTCTTGCCGAGCGCTATCAActc                                                                                                 >  1:9880/1‑148 (MQ=255)
                                             tGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTc                                                                                              <  1:204850/148‑1 (MQ=255)
                                                 aTTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCg                                                                                                               >  4:837302/1‑127 (MQ=255)
                                                 aTTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCg                                                                                                               <  3:837302/127‑1 (MQ=255)
                                                   tATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTaa                                                                                                                                >  3:818663/1‑108 (MQ=255)
                                                   tATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTaa                                                                                                                                <  4:818663/108‑1 (MQ=255)
                                                     tGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                             >  4:568195/1‑93 (MQ=255)
                                                     tGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTg                                                                                                                                             <  3:568195/93‑1 (MQ=255)
                                                      gCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGAt                                                                                     <  4:214968/148‑1 (MQ=255)
                                                       ccGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATg                                                                                    >  4:780371/1‑148 (MQ=255)
                                                             cTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                              >  3:137883/1‑148 (MQ=255)
                                                                   ctctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAg                                                                                                             >  4:421539/1‑111 (MQ=255)
                                                                   ctctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAg                                                                                                             <  3:421539/111‑1 (MQ=255)
                                                                   ctctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATg                                                                        <  3:780371/148‑1 (MQ=255)
                                                                   ctctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATg                                                                        >  4:460259/1‑148 (MQ=255)
                                                                    tctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGt                                                                                  >  3:379722/1‑137 (MQ=255)
                                                                    tctTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGt                                                                                  <  4:379722/137‑1 (MQ=255)
                                                                      ttCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccg                                                                     <  3:930101/148‑1 (MQ=255)
                                                                                tCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGa                                                                                                                                         >  2:450848/1‑70 (MQ=255)
                                                                                tCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGa                                                                                                                                         <  1:450848/70‑1 (MQ=255)
                                                                                   gggTTTACCTGGATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                        >  3:81168/1‑148 (MQ=255)
                                                                                      tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGc                                                                                         >  2:431521/1‑112 (MQ=255)
                                                                                      tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGc                                                                                         <  1:431521/112‑1 (MQ=255)
                                                                                      tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGa                                                                                      >  3:383410/1‑115 (MQ=255)
                                                                                      tttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGa                                                                                      <  4:383410/115‑1 (MQ=255)
                                                                                       ttACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCtt                                                    <  4:81168/148‑1 (MQ=255)
                                                                                         aCCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCa                                                  <  2:9880/148‑1 (MQ=255)
                                                                                            tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCa                                                                                                                                                <  3:228558/51‑1 (MQ=255)
                                                                                            tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCa                                                                                                                                                >  4:228558/1‑51 (MQ=255)
                                                                                            tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                             >  4:977153/1‑134 (MQ=255)
                                                                                            tGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                             <  3:977153/134‑1 (MQ=255)
                                                                                             gCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccg                                                                     >  2:912117/1‑125 (MQ=255)
                                                                                             gCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATgccg                                                                     <  1:912117/125‑1 (MQ=255)
                                                                                             gCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTaa                                              >  3:137288/1‑148 (MQ=255)
                                                                                              cATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTc                                                   <  4:556544/142‑1 (MQ=255)
                                                                                              cATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTc                                                   >  3:556544/1‑142 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCGTGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGTCTGATGGTCTGTTTCATGCCGCCCATCGAAt                                                           >  4:732022/1‑129 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGc                                                                                                                 >  3:842487/1‑75 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGc                                                                                                                 <  4:842487/75‑1 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAAt                                                           <  1:1005031/129‑1 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAAt                                                           >  2:1005031/1‑129 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAAt                                                           <  3:732022/129‑1 (MQ=255)
                                                                                                   cAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                        <  4:137883/148‑1 (MQ=255)
                                                                                                     aaGGCTTTGCGGATCATGTTTTCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                             >  3:308769/1‑125 (MQ=255)
                                                                                                     aaGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGa                                                             <  4:308769/125‑1 (MQ=255)
                                                                                                       ggCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACt                                          >  4:32942/1‑142 (MQ=255)
                                                                                                       ggCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACt                                          <  3:32942/142‑1 (MQ=255)
                                                                                                              cGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGAttttt                             <  3:460259/148‑1 (MQ=255)
                                                                                                               ggATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTc                            >  3:1126138/1‑148 (MQ=255)
                                                                                                                 aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGt                                                                                                                                       >  1:739203/1‑39 (MQ=255)
                                                                                                                 aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGt                                                                                                                                       <  2:739203/39‑1 (MQ=255)
                                                                                                                 aTCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTg                          <  4:137288/148‑1 (MQ=255)
                                                                                                                  tCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                >  4:671637/1‑109 (MQ=255)
                                                                                                                  tCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                <  3:671637/109‑1 (MQ=255)
                                                                                                                   cATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                        >  2:950996/1‑132 (MQ=255)
                                                                                                                   cATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                        <  1:950996/132‑1 (MQ=255)
                                                                                                                        ttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                <  2:1065763/103‑1 (MQ=255)
                                                                                                                        ttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                >  1:1065763/1‑103 (MQ=255)
                                                                                                                        ttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCa                                                  <  1:6406/117‑1 (MQ=255)
                                                                                                                        ttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCa                                                  >  2:6406/1‑117 (MQ=255)
                                                                                                                        ttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCg                        <  2:955473/143‑1 (MQ=255)
                                                                                                                        ttGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCg                        >  1:955473/1‑143 (MQ=255)
                                                                                                                          gCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGa                 <  4:1126138/148‑1 (MQ=255)
                                                                                                                            tCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                              >  1:714496/1‑85 (MQ=255)
                                                                                                                            tCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGt                                                                              <  2:714496/85‑1 (MQ=255)
                                                                                                                             ccTTATAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCg                  >  2:800525/1‑144 (MQ=255)
                                                                                                                             ccTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACt                                          >  4:223664/1‑120 (MQ=255)
                                                                                                                             ccTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACt                                          <  3:223664/120‑1 (MQ=255)
                                                                                                                             ccTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCg                  <  1:800525/144‑1 (MQ=255)
                                                                                                                             ccTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTa              >  1:708948/1‑148 (MQ=255)
                                                                                                                               ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                        >  1:378812/1‑104 (MQ=255)
                                                                                                                               ttACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATaaa                                                        <  2:378812/104‑1 (MQ=255)
                                                                                                                                tACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGa                 >  4:531114/1‑142 (MQ=255)
                                                                                                                                tACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGa                 <  3:531114/142‑1 (MQ=255)
                                                                                                                                  cagcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGtt                                                                             <  2:411653/80‑1 (MQ=255)
                                                                                                                                  cagcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGtt                                                                             >  1:411653/1‑80 (MQ=255)
                                                                                                                                   ctcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                        >  4:346621/3‑116 (MQ=255)
                                                                                                                                   ctcgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                        <  3:346621/114‑1 (MQ=255)
                                                                                                                                    gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                >  2:387121/1‑91 (MQ=255)
                                                                                                                                    gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                <  1:387121/91‑1 (MQ=255)
                                                                                                                                    gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATa                                                          >  1:322769/1‑97 (MQ=255)
                                                                                                                                    gcgccagcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATa                                                          <  2:322769/97‑1 (MQ=255)
                                                                                                                                     cgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                <  1:634023/90‑1 (MQ=255)
                                                                                                                                     cgccagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCAt                                                                >  2:634023/1‑90 (MQ=255)
                                                                                                                                        cagtgcGCTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTtg                                        >  1:69504/1‑111 (MQ=255)
                                                                                                                                         agcgcACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGttt  >  1:72849/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                              
TATGCGCACCGTGAGATTTCAGCACTGCTTCCAGTTCTGGCCAGATGGGATTATGCCGACGCTGATACTCTTCGTGGGCGTCGGGGTTTACCTGCATGACAAAGGCTTTGCGGATCATGTTTGCTCCTTACAGCGCCAGCGCACTGGCGAGTGGCGTAACGCCGAAACGCTTGCCGAGCGCTATCAACTCTTCACGGCTGATGGTCTGTTTCATGCCGCCCATCGAATAAACCTTCACTAATACTTGTGCTGATTTTTCTGCGGTGTCGATTAAACCGAAGGTTT  >  NC_000913/4093321‑4093605

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