Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 878,446 G→A 7.0% E67K (GAA→AAA)  bssR → repressor of biofilm formation by indole transport regulation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913878,4460GA 7.0% 278.8 / 8.2 85E67K (GAA→AAA) bssRrepressor of biofilm formation by indole transport regulation
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base A (3/3);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCGTTGACCTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATCCA  >  NC_000913/878300‑878594
                                                                                                                                                  |                                                                                                                                                    
gcgTTGACCTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGaa                                                                                                                                                    <  4:456175/149‑1 (MQ=255)
    tGACCTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagt                                                                                                                                                <  3:300012/149‑1 (MQ=255)
     gACCTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagtg                                                                                                                                               <  3:405084/149‑1 (MQ=255)
       ccTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagtgga                                                                                                                                             <  2:687058/149‑1 (MQ=255)
        cTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagtggag                                                                                                                                            <  2:536293/149‑1 (MQ=255)
        cTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGAAAgagtggag                                                                                                                                            >  2:543253/1‑149 (MQ=255)
        cTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGAAAgagtggag                                                                                                                                            >  4:209267/1‑149 (MQ=255)
        cTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGAAAgagtggag                                                                                                                                            >  2:92029/1‑149 (MQ=255)
         tCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagtggagt                                                                                                                                           >  2:41814/1‑149 (MQ=255)
                      gCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagtgga                                                                                                                                             >  3:722923/1‑134 (MQ=255)
                      gCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagtgga                                                                                                                                             <  4:722923/134‑1 (MQ=255)
                       cAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGc                                                                                                                             <  3:227275/149‑1 (MQ=255)
                       cAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGAAAGAGTGGAGTGCTCTTCATCATGc                                                                                                                             <  1:543253/149‑1 (MQ=255)
                       cAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGAAAGAGTGGAGTGCTCTTCATCATGc                                                                                                                             <  3:209267/149‑1 (MQ=255)
                       cAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGAAAGAGTGGAGTGCTCTTCATCATGc                                                                                                                             <  1:92029/149‑1 (MQ=255)
                           tGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTgaa                                                                                                                         >  1:385101/1‑149 (MQ=255)
                              ggAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTgaagaa                                                                                                                      <  4:23995/149‑1 (MQ=255)
                                          gATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGctct                                                                                                                                      >  3:565374/1‑121 (MQ=255)
                                          gATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGctct                                                                                                                                      <  4:565374/121‑1 (MQ=255)
                                            tACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGc                                                                                                                    >  1:487027/1‑137 (MQ=255)
                                            tACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGc                                                                                                                    <  2:487027/137‑1 (MQ=255)
                                             aCATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGcc                                                                                                       >  1:306281/1‑149 (MQ=255)
                                             aCATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGcc                                                                                                       >  4:7090/1‑149 (MQ=255)
                                               aTGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTtca                                                                                                                                   <  3:225276/119‑1 (MQ=255)
                                               aTGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTtca                                                                                                                                   >  4:225276/1‑119 (MQ=255)
                                               aTGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGc                                                                                                     >  3:669146/1‑149 (MQ=255)
                                               aTGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGc                                                                                                     >  4:737476/1‑149 (MQ=255)
                                                tGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGcc                                                                                                       <  4:560614/146‑1 (MQ=255)
                                                tGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGcc                                                                                                       >  3:560614/1‑146 (MQ=255)
                                                tGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCa                                                                                                    >  4:584217/1‑149 (MQ=255)
                                                tGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCa                                                                                                    <  3:333424/149‑1 (MQ=255)
                                                       cAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGtt                                                                                             <  4:296783/149‑1 (MQ=255)
                                                           gAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTtgat                                                                                         <  2:339322/149‑1 (MQ=255)
                                                            aaaGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTtgat                                                                                         >  3:41263/1‑148 (MQ=255)
                                                            aaaGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTtgat                                                                                         <  4:41263/148‑1 (MQ=255)
                                                                          cGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGa                                                                          <  1:797305/149‑1 (MQ=255)
                                                                             gaTAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTg                                                                       >  2:177061/1‑149 (MQ=255)
                                                                             gaTAACTTTTTTAAACTGAATCGCGAACAGCACGATAAATCGTTGCGGTTGAATCTTCATCTGTATCAGGTAGAGTGGAGTGCTCTTCATCATGGTGAAGACGCATTAGCGACAGCCGCAGTATGTTCGATGAGTGGGCACCATGATTg                                                                       >  4:580008/1‑149 (MQ=255)
                                                                              aTAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCg                                                                                                      <  4:691031/117‑1 (MQ=255)
                                                                              aTAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCg                                                                                                      >  3:691031/1‑117 (MQ=255)
                                                                                aaCTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGccc                                                                    >  2:363967/1‑149 (MQ=255)
                                                                                aaCTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGccc                                                                    <  1:581365/149‑1 (MQ=255)
                                                                                  cTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGa                                                                  >  4:755984/1‑149 (MQ=255)
                                                                                   ttttttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGAc                                                                 >  4:451195/1‑149 (MQ=255)
                                                                                   ttttttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGAc                                                                 >  1:670846/1‑149 (MQ=255)
                                                                                   ttttttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGAc                                                                 >  2:241296/1‑149 (MQ=255)
                                                                                   ttttttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGAc                                                                 >  1:460966/1‑149 (MQ=255)
                                                                                       ttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCGTCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGtt                                                             <  4:270637/149‑1 (MQ=255)
                                                                                       ttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagt                                                                                                                                                >  3:738766/1‑66 (MQ=255)
                                                                                       ttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAgagt                                                                                                                                                <  4:738766/66‑1 (MQ=255)
                                                                                       ttAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGtt                                                             <  3:584217/149‑1 (MQ=255)
                                                                                            cTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTAc                                                        <  4:789458/149‑1 (MQ=255)
                                                                                                 tCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCa                                                   >  1:754763/1‑149 (MQ=255)
                                                                                                     gAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGc                                               <  3:416492/149‑1 (MQ=255)
                                                                                                        cTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGa                                            <  3:737476/149‑1 (MQ=255)
                                                                                                        cTGCACGATAAATAGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGa                                            <  3:758435/149‑1 (MQ=255)
                                                                                                              gATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGtta                                                            <  3:275639/127‑1 (MQ=255)
                                                                                                              gATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGtta                                                            >  4:275639/1‑127 (MQ=255)
                                                                                                              gATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCt                                      >  3:273861/1‑149 (MQ=255)
                                                                                                              gATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCt                                      >  1:592177/1‑149 (MQ=255)
                                                                                                                    tCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGaaaa                                <  1:241296/149‑1 (MQ=255)
                                                                                                                    tCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGaaaa                                <  2:618573/149‑1 (MQ=255)
                                                                                                                        tGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGt                            >  1:482882/1‑149 (MQ=255)
                                                                                                                             ttGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTtgatga                                                                                       >  2:532686/1‑85 (MQ=255)
                                                                                                                             ttGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtgat                       >  1:38289/1‑149 (MQ=255)
                                                                                                                             ttGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtgat                       >  2:792110/1‑149 (MQ=255)
                                                                                                                             ttGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtgat                       <  3:129977/149‑1 (MQ=255)
                                                                                                                             ttGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtgat                       <  4:483227/149‑1 (MQ=255)
                                                                                                                             ttGAATCTTCATCTGGATCAGGAAGAGTGAAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTtgatga                                                                                       <  1:532686/85‑1 (MQ=255)
                                                                                                                              tGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtgatg                      <  2:460966/149‑1 (MQ=255)
                                                                                                                                    ttCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTa                >  1:732437/1‑149 (MQ=255)
                                                                                                                                       aTCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACg             >  1:245544/1‑149 (MQ=255)
                                                                                                                                        tCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGACACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACg             >  2:215250/1‑148 (MQ=255)
                                                                                                                                        tCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGACACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACg             <  1:215250/148‑1 (MQ=255)
                                                                                                                                         cTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCt           <  2:245544/149‑1 (MQ=255)
                                                                                                                                          tGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTg          >  1:595663/1‑149 (MQ=255)
                                                                                                                                           ggATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGa         <  2:306281/149‑1 (MQ=255)
                                                                                                                                           ggATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGa         <  1:241880/149‑1 (MQ=255)
                                                                                                                                           ggATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGa         <  2:545927/149‑1 (MQ=255)
                                                                                                                                            gATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtga                        <  4:148218/133‑1 (MQ=255)
                                                                                                                                            gATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCtga                        >  3:148218/1‑133 (MQ=255)
                                                                                                                                             aTCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGCTTGCCCGCCTGTTATTACCGTCAACGc                                               <  4:434349/109‑1 (MQ=255)
                                                                                                                                             aTCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGc                                               >  3:434349/1‑109 (MQ=255)
                                                                                                                                                aGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATc    >  2:189508/1‑149 (MQ=255)
                                                                                                                                                aGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATc    <  2:592177/149‑1 (MQ=255)
                                                                                                                                                 ggAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGc                                                                                 >  4:51955/1‑71 (MQ=255)
                                                                                                                                                 ggAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGc                                                                                 <  3:51955/71‑1 (MQ=255)
                                                                                                                                                 ggAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATcc   >  4:489211/1‑149 (MQ=255)
                                                                                                                                                  gAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATCCa  >  2:764530/1‑149 (MQ=255)
                                                                                                                                                  gAAGAGTGGAGGGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATCCa  >  4:241840/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCGTTGACCTCGCCGCATACGTGCAACTGAGGAAGGCAAAAGGATACATGTCCGTCAGCGAAAGCAATCATCTACGAGATAACTTTTTTAAACTGAATCGCGAACTGCACGATAAATCGCTGCGGTTGAATCTTCATCTGGATCAGGAAGAGTGGAGTGCTCTTCATCATGCTGAAGAAGCATTAGCGACAGCCGCAGTATGTTTGATGAGTGGGCACCATGATTGCCCGACTGTTATTACCGTCAACGCCGATAAGCTTGAAAATTGTCTGATGAGCTTAACGCTGAGTATCCA  >  NC_000913/878300‑878594

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