Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,093,965 C→G 5.2% D166E (GAC→GAG hisB → histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,093,9650CG 5.2% 329.9 / 6.0 94D166E (GAC→GAGhisBhistidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (46/43);  new base G (2/3);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

AAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCA  >  NC_000913/2093820‑2094113
                                                                                                                                                 |                                                                                                                                                    
aaGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                                                                                                                   >  2:460133/1‑149 (MQ=255)
aaGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCATGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                                                                                                                   >  2:328663/1‑149 (MQ=255)
aaGCGATGGATCGCACTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGt                                                                                                                                                   >  2:349446/1‑149 (MQ=255)
 aGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGTGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGtt                                                                                                                                                  <  1:850241/149‑1 (MQ=255)
 aGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGtt                                                                                                                                                  >  4:339830/1‑149 (MQ=255)
 aGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGtt                                                                                                                                                  <  1:536587/149‑1 (MQ=255)
  gCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTa                                                                                                                                                 >  3:594939/1‑149 (MQ=255)
   cGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTAc                                                                                                                                                >  3:476743/1‑149 (MQ=255)
   cGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTAc                                                                                                                                                <  4:68183/149‑1 (MQ=255)
          tCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACg                                                                                                                                         >  2:938299/1‑149 (MQ=255)
             cgcTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAGCGTAACGGTCACGtct                                                                                                                                      <  3:115248/149‑3 (MQ=255)
              gcTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGt                                                                                                                                     <  4:873623/149‑1 (MQ=255)
                    aGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                                                                                                                               <  3:429284/149‑1 (MQ=255)
                    aGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                                                                                                                               >  3:643962/1‑149 (MQ=255)
                           tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                        <  4:438617/149‑1 (MQ=255)
                           tgATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCaaa                                                                                                                        <  3:789472/149‑1 (MQ=255)
                            gATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAag                                                                                                                       >  2:856663/1‑149 (MQ=255)
                            gATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAag                                                                                                                       >  4:961870/1‑149 (MQ=255)
                                ggCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGAc                                                                                                                   >  2:756998/1‑149 (MQ=255)
                                     tCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGa                                                                                                              <  1:856663/149‑1 (MQ=255)
                                     tCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGa                                                                                                              <  3:245465/149‑1 (MQ=255)
                                       atgcgAACGACAGTCAAGTAGCGGAAAACATGGGCAGTACTGGTTAACGTTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGAGCGTTACGCACGCGTAGTGCGTAATACCCAAGAGACGCAGAtt                                                                                                            <  2:1018978/147‑1 (MQ=255)
                                               gacATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCa                                                                                                    >  3:830069/1‑149 (MQ=255)
                                                acATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAg                                                                                                   <  1:534526/149‑1 (MQ=255)
                                                  aTTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgt                                                                                                 <  2:391822/149‑1 (MQ=255)
                                                    tCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtgt                                                                                               <  2:48040/149‑1 (MQ=255)
                                                     cAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                                                                                              <  3:259957/149‑1 (MQ=255)
                                                           gCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGa                                                                                        >  3:395921/1‑149 (MQ=255)
                                                            cGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGAt                                                                                       <  1:184531/149‑1 (MQ=255)
                                                             ggAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                                                                                      <  2:844748/149‑1 (MQ=255)
                                                             ggAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                                                                                      >  2:1089369/1‑149 (MQ=255)
                                                              gAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCg                                                                                     >  2:47933/1‑149 (MQ=255)
                                                                aaaCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTg                                                                                   <  1:700888/149‑1 (MQ=255)
                                                                  aCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                 <  3:680489/149‑1 (MQ=255)
                                                                  aCATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                 <  4:594939/149‑1 (MQ=255)
                                                                   cATGGGCATTACTGGTTTACGCTACGACCGCGAAACCATGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAg                                                                                <  1:328663/149‑1 (MQ=255)
                                                                    aTGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAgg                                                                               <  1:756998/149‑1 (MQ=255)
                                                                     tGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGt                                                                              <  1:611971/149‑1 (MQ=255)
                                                                      gggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTg                                                                             >  4:920084/1‑149 (MQ=255)
                                                                       ggCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTgg                                                                            >  3:864219/1‑149 (MQ=255)
                                                                        gCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                 >  2:415194/1‑143 (MQ=255)
                                                                        gCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGaa                                                                                 <  1:415194/143‑1 (MQ=255)
                                                                           ttACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                                                                                            <  1:274905/129‑1 (MQ=255)
                                                                           ttACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                                                                                            >  2:274905/1‑129 (MQ=255)
                                                                             aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCaa                                                                      <  1:550803/149‑1 (MQ=255)
                                                                             aCTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCaa                                                                      >  2:1068314/1‑149 (MQ=255)
                                                                              cTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAg                                                                     >  3:92704/1‑149 (MQ=255)
                                                                                      cGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACGCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAacac                                                             >  1:858149/1‑149 (MQ=255)
                                                                                       gCTGCGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGTAGGTGGCAGCAAGATTAACAcc                                                            >  4:789101/1‑149 (MQ=255)
                                                                                        cTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCg                                                           <  2:5615/149‑1 (MQ=255)
                                                                                          aCGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGc                                                         <  4:422351/149‑1 (MQ=255)
                                                                                             aCCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGtt                                                      >  4:931560/1‑149 (MQ=255)
                                                                                                   aaaCCCTGAACTGGCCAATTATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttc                                                >  3:894380/1‑149 (MQ=255)
                                                                                                   aaaCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttc                                                >  2:880578/1‑149 (MQ=255)
                                                                                                     aCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttctt                                              <  1:873929/149‑1 (MQ=255)
                                                                                                        cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                           >  1:1016258/1‑149 (MQ=255)
                                                                                                        cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGa                                           <  3:138819/149‑1 (MQ=255)
                                                                                                         tGAACTGGCCAACGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                          <  3:819235/149‑1 (MQ=255)
                                                                                                            aCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCat                                       <  1:47933/149‑1 (MQ=255)
                                                                                                              tGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCatat                                     >  4:458843/1‑149 (MQ=255)
                                                                                                                gCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGc                                   <  3:829573/149‑1 (MQ=255)
                                                                                                                 ccAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCt                                  >  1:606851/1‑149 (MQ=255)
                                                                                                                  cAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttct                                               >  2:988965/1‑135 (MQ=255)
                                                                                                                  cAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGcttct                                               <  1:988965/135‑1 (MQ=255)
                                                                                                                  cAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTg                                 >  3:181189/1‑149 (MQ=255)
                                                                                                                    aTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGAACGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGa                               <  4:820006/149‑1 (MQ=255)
                                                                                                                     tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGAt                              <  1:938299/149‑1 (MQ=255)
                                                                                                                        ttGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTg                                                                                   >  3:1016706/1‑93 (MQ=255)
                                                                                                                        ttGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTg                                                                                   <  4:1016706/93‑1 (MQ=255)
                                                                                                                         tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGa                          <  1:460133/149‑1 (MQ=255)
                                                                                                                         tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGa                          <  2:63967/149‑1 (MQ=255)
                                                                                                                               gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTa                    <  1:129745/149‑1 (MQ=255)
                                                                                                                               gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTa                    >  4:946647/1‑149 (MQ=255)
                                                                                                                               gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTa                    >  3:385187/1‑149 (MQ=255)
                                                                                                                                 aaCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTAcc                  >  2:384843/1‑149 (MQ=255)
                                                                                                                                  aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                          <  3:762477/124‑1 (MQ=255)
                                                                                                                                  aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                          >  4:762477/1‑124 (MQ=255)
                                                                                                                                  aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGc                                   >  4:387626/1‑131 (MQ=255)
                                                                                                                                  aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGc                                   <  3:387626/131‑1 (MQ=255)
                                                                                                                                  aCTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTAccc                 >  1:322975/1‑149 (MQ=255)
                                                                                                                                      aCCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcgg             >  2:617001/1‑149 (MQ=255)
                                                                                                                                       ccAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcg           >  2:578465/1‑149 (MQ=255)
                                                                                                                                        cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcg           >  4:38308/1‑149 (MQ=255)
                                                                                                                                        cAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcg           >  2:234186/1‑149 (MQ=255)
                                                                                                                                         aGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCAcggcgg          <  1:998404/149‑1 (MQ=255)
                                                                                                                                          taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgc                                                                           >  3:760527/2‑83 (MQ=39)
                                                                                                                                          taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgc                                                                           <  4:760527/82‑1 (MQ=39)
                                                                                                                                          taCGTGAGCGTAACGGTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACCGCGGt         >  2:579097/2‑149 (MQ=255)
                                                                                                                                          gACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGt         >  4:1075166/1‑149 (MQ=255)
                                                                                                                                             gTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTc      >  4:1079585/1‑149 (MQ=255)
                                                                                                                                             gTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTc      >  3:792698/1‑149 (MQ=255)
                                                                                                                                                aCCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                          <  3:247917/110‑1 (MQ=255)
                                                                                                                                                aCCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGAt                                          >  4:247917/1‑110 (MQ=255)
                                                                                                                                                aCCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGc   <  4:792698/149‑1 (MQ=255)
                                                                                                                                                 ccGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCa  <  2:390509/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
AAGCGATGGATCGCGCTAACAGTTATGTGATTGGCGATCGCGCGACCGACATTCAACTGGCGGAAAACATGGGCATTACTGGTTTACGCTACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAAGATTAACACCGGCGTTGGCTTCTTTGATCATATGCTGGATCAGATCGCTACCCACGGCGGTTTCCGCA  >  NC_000913/2093820‑2094113

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