Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 886,947 A→G 5.1% V131V (GTT→GTC ybjJ ← putative drug efflux MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913886,9470AG 5.1% 244.3 / 3.1 78V131V (GTT→GTCybjJputative drug efflux MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base G (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

GATACCTACCAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATC  >  NC_000913/886800‑887091
                                                                                                                                                   |                                                                                                                                                
gATACCTACCAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACTCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCa                                                                                                                                                  <  1:210056/148‑1 (MQ=255)
    cctaccAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACg                                                                                                                                               <  3:370631/147‑1 (MQ=255)
     ctaccAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGgc                                                                                                                                             <  8:168771/148‑1 (MQ=255)
     ctaccAGCACCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGgc                                                                                                                                             <  5:19013/148‑1 (MQ=255)
       accAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGgcag                                                                                                                                           >  5:381696/1‑148 (MQ=255)
        ccAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGgcagc                                                                                                                                          <  4:336314/148‑1 (MQ=255)
                    aTAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGtt                                                                                                                              >  5:88345/1‑148 (MQ=255)
                    aTAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATAATCGGCAAAACCGTCTTATTAATTTCTAGCTCAACGGAATCACATTCAAAGtt                                                                                                                              >  5:379320/1‑148 (MQ=255)
                      aaaaTGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTACGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTa                                                                                                                            >  3:168650/1‑148 (MQ=255)
                       aaaTGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCGATAAAAACCGTGCATCATCGGCAAAACCGGCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTAt                                                                                                                           <  6:123657/148‑1 (MQ=255)
                        aaTGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATc                                                                                                                          <  2:234101/148‑1 (MQ=255)
                         aTGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAACACCGTGCATCATCGGCAACACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCg                                                                                                                         >  7:107086/1‑148 (MQ=255)
                                ccGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTc                                                                                                                  <  5:95595/148‑1 (MQ=255)
                                  gTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCGACGGCAGCACCTTCAACGTTTATCGcc                                                                                                                       >  8:308558/1‑141 (MQ=255)
                                  gTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCGACGGCAGCACCTTCAACGTTTATCGcc                                                                                                                       >  8:314092/1‑141 (MQ=255)
                                  gTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCGACGGCAGCACCTTCAACGTTTATCGcc                                                                                                                       <  7:314092/141‑1 (MQ=255)
                                  gTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCGACGGCAGCACCTTCAACGTTTATCGcc                                                                                                                       <  7:308558/141‑1 (MQ=255)
                                      ccGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAg                                                                                                             <  5:123277/147‑1 (MQ=255)
                                          aaCGCCAAAGGCCGTCAGTGCCATCCCGACACCTACGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCa                                                                                                        <  4:168650/148‑1 (MQ=255)
                                            cGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCaaa                                                                                                      <  4:187515/148‑1 (MQ=255)
                                                aaaGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACtt                                                                                                  <  3:196981/148‑1 (MQ=255)
                                                       gTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCCGCACTTCCGCAGAACCAAAACTTGCCGCa                                                                                            >  8:177365/1‑147 (MQ=255)
                                                          aGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGa                                                                                        <  8:11984/148‑1 (MQ=255)
                                                             gCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAACCCGTGCATCATCGGCAAAACCGTCTGATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGc                                                                                     >  1:229961/1‑148 (MQ=255)
                                                              ccATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGcc                                                                                    <  1:133888/148‑1 (MQ=255)
                                                              ccATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGcc                                                                                    >  2:293608/1‑148 (MQ=255)
                                                               cATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAACGACGCCg                                                                                   <  5:36812/148‑1 (MQ=255)
                                                                 tCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGcc                                                                                    >  3:290832/1‑145 (MQ=255)
                                                                  cccGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCaaaa                                                                                                     >  8:188447/1‑127 (MQ=255)
                                                                  cccGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCaaaa                                                                                                     <  7:188447/127‑1 (MQ=255)
                                                                  cccGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCgaga                                                                                >  4:212737/1‑148 (MQ=255)
                                                                  cccGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCgag                                                                                 >  2:338203/1‑147 (MQ=255)
                                                                  cccGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCgag                                                                                 <  1:338203/147‑1 (MQ=255)
                                                                     gACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGAcc                                                                              >  3:122531/1‑147 (MQ=255)
                                                                      acacCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCGCCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCaa                                                                            >  2:91089/1‑148 (MQ=255)
                                                                          cTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACgg                                                                         >  1:83257/1‑147 (MQ=255)
                                                                           tGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGc                                                                        >  3:170711/1‑147 (MQ=255)
                                                                            gcgcCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCaa                                                                      <  3:17682/148‑1 (MQ=255)
                                                                               ccTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAaga                                                                   <  4:290832/148‑1 (MQ=255)
                                                                                   ccAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGc                                                                 <  6:88345/146‑1 (MQ=255)
                                                                                   ccAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGc                                                                 <  6:132176/146‑1 (MQ=255)
                                                                                   ccAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGc                                                                 <  8:396524/146‑1 (MQ=255)
                                                                                     aGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAggg                                                             <  3:212737/148‑1 (MQ=255)
                                                                                               ggCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGa                                                         >  4:13161/1‑142 (MQ=255)
                                                                                               ggCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGa                                                         <  3:13161/142‑1 (MQ=255)
                                                                                                 cTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGcc                                                 <  1:293608/148‑1 (MQ=255)
                                                                                                 cTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGc                                                  >  4:124704/1‑147 (MQ=255)
                                                                                                    taAAAACCGTGCATCATCGGCTAAACCGTCTTATTGATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAATCTTGCCCCTAAGTCGCCGAGACGAACGGCAAAGAGCAGGGGAGATGTCAGCCaga                                              <  8:396780/148‑1 (MQ=255)
                                                                                                    taAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCaga                                              <  1:168778/148‑1 (MQ=255)
                                                                                                     aaaaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAgg                                                              >  1:377589/1‑131 (MQ=255)
                                                                                                     aaaaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAgg                                                              <  2:377589/131‑1 (MQ=255)
                                                                                                     aaaaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCg                                                          <  8:185909/135‑1 (MQ=255)
                                                                                                     aaaaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCg                                                          >  7:185909/1‑135 (MQ=255)
                                                                                                     aaaaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCagag                                             >  6:147305/1‑148 (MQ=255)
                                                                                                       aaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCaa                                                                      >  8:275322/1‑121 (MQ=255)
                                                                                                       aaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCaa                                                                      <  7:275322/121‑1 (MQ=255)
                                                                                                       aaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGTCAGAGTg                                           >  3:327072/1‑148 (MQ=255)
                                                                                                       aaaCCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTg                                           >  5:130889/1‑148 (MQ=255)
                                                                                                           cGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAg                                       <  4:122531/148‑1 (MQ=255)
                                                                                                                  catcGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGc                                                                 >  1:307797/1‑115 (MQ=255)
                                                                                                                  catcGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGc                                                                 <  2:307797/115‑1 (MQ=255)
                                                                                                                   atcGGCAAAACAGTCTTATTCATTTCTCGGTCACCGGCAGCACCTTCAACGTTTATCGGGACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAat                                <  6:379320/147‑1 (MQ=255)
                                                                                                                     cGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTg                                                                                                 >  7:100324/1‑80 (MQ=255)
                                                                                                                     cGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTg                                                                                                 <  8:100324/80‑1 (MQ=255)
                                                                                                                                    aTTCATTTCTCGCTCAACGGCAGCACCTTCATCGTTTATCGCCACTTCCGCAGAAACAAAAGTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCa                <  4:119090/146‑1 (MQ=255)
                                                                                                                                     ttCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCg                                                          <  5:242852/103‑1 (MQ=255)
                                                                                                                                     ttCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCg                                                          >  6:242852/1‑103 (MQ=255)
                                                                                                                                         tttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCa                                                                <  8:136694/93‑1 (MQ=255)
                                                                                                                                         tttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCa                                                                >  7:136694/1‑93 (MQ=255)
                                                                                                                                         tttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGAc                                     >  6:104316/1‑120 (MQ=255)
                                                                                                                                         tttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGAc                                     <  5:104316/120‑1 (MQ=255)
                                                                                                                                         tttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACCTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCa         <  2:227265/148‑1 (MQ=255)
                                                                                                                                          ttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCATCGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATatcatc                          <  4:50905/130‑1 (MQ=255)
                                                                                                                                          ttCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCATCGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATatcatc                          >  3:50905/1‑130 (MQ=255)
                                                                                                                                              cGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGAc     >  2:354716/1‑147 (MQ=255)
                                                                                                                                              cGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACa    >  1:57538/1‑148 (MQ=255)
                                                                                                                                               gCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGTCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACAt   <  4:327072/148‑1 (MQ=255)
                                                                                                                                                cTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATc  <  2:57538/148‑1 (MQ=255)
                                                                                                                                                  cAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCAAGAGTGCCAGACTTAATATAATCATCCCTAACAATGCGCAGGACATc  >  7:252693/1‑146 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GATACCTACCAGCGCCGCCAATAAAATGTGCACCGTTGCCGGAACGCCAAAGGCCGTCAGTGCCATCCCGACACCTGCGCCTGCCAGCGTGCCCAGGCTATAAAAACCGTGCATCATCGGCAAAACCGTCTTATTCATTTCTCGCTCAACGGCAGCACCTTCAACGTTTATCGCCACTTCCGCAGAACCAAAACTTGCCCCAAAGACGCCGAGACCAACGGCAAAGAGCAGGGGCGATGTCAGCCAGAGTGCCAGACTTAATATCATCATCCCGATCAATGCGCAGGACATC  >  NC_000913/886800‑887091

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