Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 180,794 T→C 6.8% intergenic (+40/‑90) dgt → / → degP deoxyguanosine triphosphate triphosphohydrolase/serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913180,7940TC 6.8% 161.1 / 3.1 61intergenic (+40/‑90)dgt/degPdeoxyguanosine triphosphate triphosphohydrolase/serine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base C (2/2);  total (33/28)
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.17e-01

CTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGTTTAGGTTTGGC  >  NC_000913/180647‑180936
                                                                                                                                                   |                                                                                                                                                
cTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTg                                                                                                                                                 <  1:7575/149‑1 (MQ=255)
        atatTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACt                                                                                                                                         >  1:722680/1‑149 (MQ=255)
         tattatTCCCGTTGCCGCCTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACtt                                                                                                                                        <  1:103769/149‑1 (MQ=255)
              ttaCCGTTGCCGCCTGCGGCAGGATTATACCAGCTGTATGACCCACCTCTACTCGTGGGATGAATACCGACGTCGGATGGCCGTAGAACAATAACCAGGCTTTTGTAAATACGAACAATAAATTTTTACCTTTTGCAGAAACCTTAGtt                                                                                                                                   <  2:292903/149‑1 (MQ=255)
                   gTTGCCGCCTGCTGCAGGATTATATCAGTGTTATGACGGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAACACATTTTTACCTTTTGCAGAAACTTTAGTTCAGac                                                                                                                              <  2:562241/149‑2 (MQ=255)
                       ccgccTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTGTTGCAGAAACTTTAGGGCGGAACTTc                                                                                                                          >  2:631218/1‑149 (MQ=255)
                        cgccTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCCTTTGTCAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTcc                                                                                                                         >  1:685673/1‑148 (MQ=255)
                          ccTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGCACAATAAATTTTTACCTTCTGGAGAAACTTTAGTTCGGAACTTCAgg                                                                                                                       <  1:448441/149‑1 (MQ=255)
                          ccTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAgg                                                                                                                       >  2:161715/1‑149 (MQ=255)
                          ccTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTGTTGCAGAAACTTTAGTTCGGAACTTCAgg                                                                                                                       >  2:161567/1‑149 (MQ=255)
                                 cAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTCGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATaaa                                                                                                               >  2:566185/1‑149 (MQ=255)
                                 cAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTCGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATaaa                                                                                                               >  2:575885/1‑149 (MQ=255)
                                    gaatataTCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCGTTTGCAGAAACTTTAGTTCGCAACTTCAGGCT‑ATAAAACg                                                                                                            >  1:722331/4‑149 (MQ=255)
                                    gATTATATCAGCGGTATGACCGAGCTATATGCGTGGGATGAAGACCGACGTCTGATGGCCGTAGAACGATAGCCAGGCTGTTGTAAAAACGAACAATACGTCTTTGCCGTTTGCAGAAAGTTTCGGTCGGACCTTCAGGTC‑ATATAACg                                                                                                            >  2:446621/1‑149 (MQ=255)
                                    gATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACg                                                                                                            >  1:89013/1‑149 (MQ=255)
                                       tataTCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAAt                                                                                                         <  2:658261/149‑1 (MQ=255)
                                        atatCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCAGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑CTAAAACGAATc                                                                                                        >  1:615453/1‑149 (MQ=255)
                                              gCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTgaaga                                                                                                  >  1:185479/1‑149 (MQ=255)
                                              gCGGTATGACCGACCTCTATGCGTGGGAAGAAGACCGATGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTgcaga                                                                                                  >  2:174425/1‑149 (MQ=255)
                                               cGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGGTGGCCGTAGAACAATAACCAGGCTTTTGTACAGCCGAACAATGAATTTTTACCTTTTGCAGAAACTTTAGTGCGGAACTTCAGGCT‑ATAAAACGAATCTgaagaa                                                                                                 >  2:684288/1‑149 (MQ=255)
                                               cGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTgaagaa                                                                                                 >  2:162038/1‑149 (MQ=255)
                                                      accgaccTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAGCTTTAGTTCGGAACTTAAGACT‑ATAAAACGAATCTGAAGAACACAGCa                                                                                          >  2:348495/1‑149 (MQ=255)
                                                        cgaccTCTATGTGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGGAGAACACAGCagt                                                                                        >  1:21570/1‑147 (MQ=255)
                                                         gaccTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCATGCTTTTGTAAAGACGAACAATAGATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAAtt                                                                                       <  2:374989/149‑1 (MQ=255)
                                                         gaccTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAAtt                                                                                       >  2:388495/1‑149 (MQ=255)
                                                            ctctATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTCCCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATATAACGAATCTGAAGAACACAGCAATTTTg                                                                                    <  1:480664/149‑1 (MQ=255)
                                                            cactATGCATGGGATGAATACCGACGTCGGATGGCCGTAGGACAAGAACCATGCGTTGGGAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAACTTTg                                                                                    >  2:469016/3‑149 (MQ=255)
                                                              ctATGCGTGGGATGAATACCGACTTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCg                                                                                  >  1:443833/1‑149 (MQ=255)
                                                                   cGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCACTTTTGCGTTATc                                                                             >  1:32276/1‑149 (MQ=255)
                                                                   cGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATc                                                                             <  2:375639/149‑1 (MQ=255)
                                                                   cGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAAGTTTGCGTTATc                                                                             >  1:32190/1‑149 (MQ=255)
                                                                    gcaggATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATcc                                                                            >  2:336497/4‑148 (MQ=255)
                                                                            taTACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCg                                                                    <  2:385820/148‑1 (MQ=255)
                                                                             tcaCCGATGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACGTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCga                                                                   <  2:631758/147‑1 (MQ=255)
                                                                             aTACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCga                                                                   <  1:283940/149‑1 (MQ=255)
                                                                               aCCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCgaga                                                                 >  1:576937/1‑149 (MQ=255)
                                                                                          aTGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACAt                                                      <  1:722062/149‑1 (MQ=255)
                                                                                          aTGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACAt                                                      <  1:333526/149‑1 (MQ=255)
                                                                                             gCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑aa                                                  >  1:437951/1‑149 (MQ=255)
                                                                                               cGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGGTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGATATCTGTTAATCGCGACGGAAATACATG‑acaa                                                >  1:735462/1‑149 (MQ=255)
                                                                                                   gAACAATCGCCAGGCTTTTGTAAACACGAACAATAAAGTGTTCCCCTTTGCAGAAAGTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAAc                                            <  1:411530/149‑1 (MQ=255)
                                                                                                     aCAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACca                                          <  1:11758/149‑1 (MQ=255)
                                                                                                            ccAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATATGAAGAACACAGCAATTTTACGTAATCTATTAATCGAAACTAAAGTGCATT‑AAAAAAAGAAAATTAGc                                   >  1:364291/1‑149 (MQ=255)
                                                                                                             cAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCa                                  >  1:344738/1‑149 (MQ=255)
                                                                                                             cAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCa                                  >  2:544066/1‑149 (MQ=255)
                                                                                                                 ccgttatAAAGACGAACAATAAATTATTACCATCTCCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCCGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTAAAATACAAA‑AAAAAAACCACATTAGCACTGa                              <  2:42978/143‑1 (MQ=255)
                                                                                                                 cTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGa                              <  1:285392/149‑1 (MQ=255)
                                                                                                                  ttttGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATGAAAAAAAACCACATTAGCACTGa                              <  1:559617/149‑1 (MQ=255)
                                                                                                                        aaaGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGCACTTCAGGCT‑ATAACACGAATCTGAAGAACACAGCAATTTTGCGTTAACTGTTAATCGAGACTGAGATACCTG‑AAAAAAACCACATTAGCACTGAGTGCCCt                       >  1:433426/1‑149 (MQ=255)
                                                                                                                          aGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTgg                     <  2:185479/149‑1 (MQ=255)
                                                                                                                           gACGAACAATAAATTTTTACCTTTCGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGc                    <  1:566185/149‑1 (MQ=255)
                                                                                                                           gACGAACAATAAATTTTTACCTTTCGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGc                    <  1:575885/149‑1 (MQ=255)
                                                                                                                               aaCTATGAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTg                <  1:298680/149‑1 (MQ=255)
                                                                                                                                 cAATAAATTTTTACCTTTAACAGAAACTTTAGTACGGAA‑ACCATGTTGCCAAAATGCATCTTAAGGACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTggg              >  1:586827/1‑147 (MQ=255)
                                                                                                                                  aagaaaTTTTTCCCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGt             <  1:282121/146‑1 (MQ=255)
                                                                                                                                  aaTAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTTAAGAACACAGCAATTTTGCGTTATATGTTAATCGAGACTGAAATACATG‑ATAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGt             <  1:517838/149‑1 (MQ=255)
                                                                                                                                      aaTTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTACTCGAGACCGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGTTTAg         <  1:550977/149‑1 (MQ=255)
                                                                                                                                      aaTTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGTTTAg         >  1:777755/1‑149 (MQ=255)
                                                                                                                                      aaTTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTAAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGTTTAg         >  1:625175/1‑149 (MQ=255)
                                                                                                                                            tACCTTTTGCAGAAACTTTAGGTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAAGCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGAGCACTGGCTCTGAGTTTAGGTTTgg   <  1:399591/149‑1 (MQ=255)
                                                                                                                                             aCCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGTTTAGGTTTGGc  <  1:162038/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CTATGGGAATATTATTACCGTTGCCGCCTGCTGCAGGATTATATCAGCGGTATGACCGACCTCTATGCGTGGGATGAATACCGACGTCTGATGGCCGTAGAACAATAACCAGGCTTTTGTAAAGACGAACAATAAATTTTTACCTTTTGCAGAAACTTTAGTTCGGAACTTCAGGCT‑ATAAAACGAATCTGAAGAACACAGCAATTTTGCGTTATCTGTTAATCGAGACTGAAATACATG‑AAAAAAACCACATTAGCACTGAGTGCACTGGCTCTGAGTTTAGGTTTGGC  >  NC_000913/180647‑180936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑