Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,959 T→G 12.7% N110T (AAC→ACC)  hslV ← peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,9590TG12.7% 306.2 / 39.7 102N110T (AAC→ACC) hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (45/44);  new base G (7/6);  total (52/50)
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.44e-01

Only 100 of 103 total aligned reads displayed.

CTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAA  >  NC_000913/4121819‑4122106
                                                                                                                                            |                                                                                                                                                   
cTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                              >  3:223647/1‑148 (MQ=255)
 tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg                                                                                                                                               >  1:762338/1‑146 (MQ=255)
 tGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCgg                                                                                                                                               <  2:762338/146‑1 (MQ=255)
         ccAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAg                                                                                                                                     >  4:304500/1‑148 (MQ=255)
            aCGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGa                                                                                                                                  <  3:304500/148‑1 (MQ=255)
             cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGAt                                                                                                                                 <  4:223647/148‑1 (MQ=255)
                      aGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCa                                                                                                                        <  2:356788/148‑1 (MQ=255)
                        cAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATc                                                                                                                      >  1:126690/1‑148 (MQ=255)
                                        aaGTTCAGTGTTTTCTAACAGCGCGCGCGCGTCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGAACGTTGTCTGGCTGCCCCACGGCACAGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGt                                                                                                      >  1:508427/1‑148 (MQ=255)
                                        aaGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg                                                                                                                                             >  1:648938/1‑109 (MQ=255)
                                        aaGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtg                                                                                                                                             <  2:648938/109‑1 (MQ=255)
                                          gTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGc                                                                                                    <  2:126690/148‑1 (MQ=255)
                                           ttCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGc                                                                                                                                    <  3:857077/115‑1 (MQ=255)
                                           ttCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGc                                                                                                                                    >  4:857077/1‑115 (MQ=255)
                                                   tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAg                                                                                                                                     >  3:807544/1‑106 (MQ=255)
                                                   tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAg                                                                                                                                     <  4:807544/106‑1 (MQ=255)
                                                   tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCg                                                                                                                                   <  1:840642/108‑1 (MQ=255)
                                                   tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCg                                                                                                                                   >  2:840642/1‑108 (MQ=255)
                                                         gcggcgcgcgcgcTGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGAt                                                                                                                                 <  2:304774/101‑1 (MQ=255)
                                                         gcggcgcgcgcgcTGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGAt                                                                                                                                 >  1:304774/4‑104 (MQ=255)
                                                          cAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                                                                                                                         >  4:127517/1‑95 (MQ=255)
                                                          cAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                                                                                                                         <  3:127517/95‑1 (MQ=255)
                                                          cAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcgcg                                                                                                                  >  1:1071538/1‑118 (MQ=255)
                                                          cAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcgcg                                                                                                                  <  2:1071538/118‑1 (MQ=255)
                                                           aGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGt                                                                                                      >  3:23812/1‑129 (MQ=255)
                                                           aGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGt                                                                                                      <  4:23812/129‑1 (MQ=255)
                                                            gcgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTc                                                                                                                         >  1:672861/1‑109 (MQ=255)
                                                            gcgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTc                                                                                                                         <  2:672861/109‑1 (MQ=255)
                                                            gcgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTc                                                                                                                         <  2:712643/109‑1 (MQ=255)
                                                            gcgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTc                                                                                                                         >  1:712643/1‑109 (MQ=255)
                                                                gcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCtt                                                                                                  >  3:1067840/1‑128 (MQ=255)
                                                                gcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCtt                                                                                                  <  4:1067840/128‑1 (MQ=255)
                                                                 cgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGAt                                                                             >  4:111384/1‑148 (MQ=255)
                                                                     ccgcAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGt                                                                         <  3:111384/148‑1 (MQ=255)
                                                                                     aGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             <  1:665460/144‑1 (MQ=255)
                                                                                     aGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             >  2:665460/1‑144 (MQ=255)
                                                                                      gggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTa                                                                        >  3:179808/1‑132 (MQ=255)
                                                                                      gggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTa                                                                        <  4:179808/132‑1 (MQ=255)
                                                                                        gccgccGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                                                                                                                         >  3:934195/1‑65 (MQ=255)
                                                                                        gccgccGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatga                                                                                                                                         <  4:934195/65‑1 (MQ=255)
                                                                                           gccgGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                    <  2:396637/83‑1 (MQ=255)
                                                                                           gccgGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                    >  1:396637/1‑83 (MQ=255)
                                                                                                     aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCa                                                                                                                        >  3:628693/1‑69 (MQ=255)
                                                                                                     aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCa                                                                                                                        <  4:628693/69‑1 (MQ=255)
                                                                                                     aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCt                                          >  1:108517/1‑147 (MQ=255)
                                                                                                     aTAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCt                                          <  2:108517/147‑1 (MQ=255)
                                                                                                        gCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGa                                                                                                                                  <  1:8172/56‑1 (MQ=255)
                                                                                                        gCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGa                                                                                                                                  >  2:8172/1‑56 (MQ=255)
                                                                                                           aTAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGcc                                                         >  3:671326/1‑126 (MQ=255)
                                                                                                           aTAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGcc                                                         <  4:671326/126‑1 (MQ=255)
                                                                                                                aTCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGt                                                                                             <  4:754831/85‑1 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGtt                                                                                            <  3:991910/85‑1 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             >  3:170237/1‑116 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             <  3:856355/116‑1 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             >  4:856355/1‑116 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             <  4:170237/116‑1 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAg                                  >  1:291963/1‑143 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGc                             >  1:777491/1‑148 (MQ=255)
                                                                                                                 tCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGc                             >  2:663261/1‑148 (MQ=255)
                                                                                                                  cGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGcc                            >  1:1027291/1‑148 (MQ=255)
                                                                                                                          ggCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCa                    <  1:663261/148‑1 (MQ=255)
                                                                                                                            cTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGccc                           >  4:860840/1‑139 (MQ=255)
                                                                                                                            cTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGccc                           <  3:860840/139‑1 (MQ=255)
                                                                                                                                accacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcagca                                                                                                        <  4:379562/58‑1 (MQ=255)
                                                                                                                                accacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCcagca                                                                                                        >  3:379562/1‑58 (MQ=255)
                                                                                                                                 ccacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCtt                                                              >  1:460284/1‑99 (MQ=255)
                                                                                                                                 ccacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCtt                                                              <  2:460284/99‑1 (MQ=255)
                                                                                                                                 ccacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTg                         <  2:593741/136‑1 (MQ=255)
                                                                                                                                 ccacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTg                         >  1:593741/1‑136 (MQ=255)
                                                                                                                                 ccacGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAg             >  4:664406/1‑148 (MQ=255)
                                                                                                                                   acGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             <  3:802780/98‑1 (MQ=255)
                                                                                                                                   acGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             >  4:802780/1‑98 (MQ=255)
                                                                                                                                   acGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGcc                                                         >  1:616912/1‑102 (MQ=255)
                                                                                                                                   acGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCCTTTGcc                                                         <  2:616912/102‑1 (MQ=255)
                                                                                                                                     gTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTa         <  2:777491/148‑1 (MQ=255)
                                                                                                                                     gTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTa         <  2:1027291/148‑1 (MQ=255)
                                                                                                                                     atCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTg                                       >  2:678420/2‑118 (MQ=255)
                                                                                                                                     atCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTg                                       <  1:678420/117‑1 (MQ=255)
                                                                                                                                      tCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                       >  1:266750/1‑85 (MQ=255)
                                                                                                                                      tCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTAc                                                                       <  2:266750/85‑1 (MQ=255)
                                                                                                                                      tCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTGAACGGCGGCTTTg                                       >  3:723978/1‑117 (MQ=255)
                                                                                                                                      tCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGTTTTg                                       <  4:723978/117‑1 (MQ=255)
                                                                                                                                      tCACCGGTAACGGTGTTGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCt                                                                                                   <  3:698702/57‑1 (MQ=255)
                                                                                                                                      tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCt                                                                                                   >  4:698702/1‑57 (MQ=255)
                                                                                                                                      tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             <  1:910803/95‑1 (MQ=255)
                                                                                                                                      tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCttt                                                             >  2:910803/1‑95 (MQ=255)
                                                                                                                                      tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTg                         <  2:162916/131‑1 (MQ=255)
                                                                                                                                      tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTg                         >  1:162916/1‑131 (MQ=255)
                                                                                                                                       cACCGGTAGCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCa                                                        >  2:865545/1‑99 (MQ=255)
                                                                                                                                       cACCGGTAGCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCa                                                        <  1:865545/99‑1 (MQ=255)
                                                                                                                                       cACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcgg                                               <  3:865119/108‑1 (MQ=255)
                                                                                                                                       cACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcgg                                               >  4:865119/1‑108 (MQ=255)
                                                                                                                                       cACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACg       >  1:64030/1‑148 (MQ=255)
                                                                                                                                         ccGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTg                                                                                                              <  4:297935/43‑1 (MQ=255)
                                                                                                                                         ccGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTg                                                                                                              >  3:297935/1‑43 (MQ=255)
                                                                                                                                         ccGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATc                                                                            <  4:630098/77‑1 (MQ=255)
                                                                                                                                         ccGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATc                                                                            >  3:630098/1‑77 (MQ=255)
                                                                                                                                         ccGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATc                                                                            >  2:1034658/1‑77 (MQ=255)
                                                                                                                                         ccGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATc                                                                            <  1:1034658/77‑1 (MQ=255)
                                                                                                                                            ttACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCaa  >  4:112795/1‑148 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
CTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAA  >  NC_000913/4121819‑4122106

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: